USAARL  Contract  Report  No.  CR  95-3 


Development  of  a  Standard  for  the  Health  Hazard  Assessment 
of  Mechanical  Shock  and  Repeated  Impact  in  Army  Vehicles 

Phase  3,  Pilot  Tests 
By 

Judy  Village 
George  Roddan 
Brian  Remedios 
James  Morrison 
Julia  Rylands 
Dale  Brown 
Dan  Robinson 
Barbara  Cameron 

B.C.  Research 
Vancouver,  B.C.,  Canada 


United  States  Army  Aeromedical  Research  Laboratory 
Fort  Rucker,  Alabama  36362-0577 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


1b.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 

Approved  for  piiblic  release,  distribution 
unlimited 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
USAARL  Contract  Report  No.  CR  95-3 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
U.S.  Army  Aeromedical 
Research  Laboratory 


6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

P.O.  Box  620577 

Fort  Rucker,  AL  36362-0577 


6b.  OFFICE  SYMBOL 
(If  applicable) 

MCMR-UAD 


7a.  NAME  OF  MONITORING  ORGANIZATION 
U.S.  Army  Medical  Research  and  Materiel 
Command 


7b.  ADDRESS  (City,  Slate,  and  ZIP  Code) 
Fort  Detrick 

Frederick,  MD  21702-5012 


8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 


8b.  OFFICE  SYMBOL 
(If  applicable) 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
DAMD17-91-C-1115 


8c.  ADDRESS  (City,  Stale,  and  ZIP  Code) 


PROGRAM 
ELEMENT  NO. 


62787A 


PROJECT 

TASK 

NO. 

NO. 

30162787A878 

FA 

WORK  UNIT 
ACCESSION  NO. 


1 1 .  TITLE  (Include  Security  Classificalion) 

Development  of  a  standard  for  the  health  hazard  assessment  of  mechanical  shock  and 
repeated  impact  in  Army  vehicles.  Phase  3 


12.  PERSONAL  AUTHOR(S) 

J.  Village, J.  Morrison,  D.  Robinson, G.  Roddan,J.  Rylands,B.  Cameron,  B.  Remedios,  D.  Brovm 


13a.  TYPE  OF  REPORT 
Final 


13b.  TIME  COVERED 
FROM 


14.  DATE  OF  REPORT  fYear,  Month,  Day) 
1995  May  _ 


15.  PAGE  COUNT 
267 


16.  SUPPLEMENTAL  NOTATION 

Fulfills  one  of  the  requirements  of  Contract  DAMD17-91-C-1115 


COSATT  CODES 


GROUP  I  SUB-GROUP 


18.  SUBJECT  TERMS  (Conbnue  on  teverae  If  necessary  and  identify  by  block  number) 
mechanical  shock,  repeated  impact,  vibration  exposure, 
jolt,  jolt  steuidards,  biomechanic  modeling 


1 9.  /VBSTRACT  (Conbnue  on  reverse  If  necessary  and  identify  by  block  number) 

New  tactical  ground  vehicles  developed  by  the  U.S.  Army  are  lower  in  weight  and 
capable  of  higher  speeds  than  their  predecessors.  This  combination  produces  repetitive 
mechanical  shocks  that  are  transmitted  to  the  soldier  primarily  through  the  seating 
system.  Under  certain  operating  conditions,  this  exposure  poses  health  and  safety  threats 
to  the  crew  as  well  as  performcuice  degradation  due  to  fatigue.  The  Army  Surgeon  General 
urgently  required  the  Medical  Research  6uid  Materiel  Command  to  develop  exposure  steuidards 
for  repetitive  impacts  that  are  relevant  to  the  environment  of  soldiers  operating  modem 
tactical  vehicles. 

A  five-phase  research  study  was  designed  to  develop  a  standard  for  the  health  hazard 
assessment  of  mechcuiical  shock  and  repeated  impact  in  Army  vehicles.  Phase  1  reviewed  the 
relevant  scientific,  medical,  and  military  literature.  Phase  2  analyzed  and  characterized 
the  shock  and  vibration  environment  of  Army  tactical  ground  vehicles  (TGVs) .  Phase  3,  a 
pilot  study,  determined  the  most  sensitive  human  response  measures  to  mechcUiical  shock  emd 

(continued  on  next  page) 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
^  UNCLASSIFIEDAJNLIMITED  EH  SAMEASRPT.  U  PRC  USERS 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Chief,  Science  Support  Center 


21.  ABSTRACT  SECURITY  CLASSIFICARON 
Unclassified 


22b.  TELEPHONE  (Include  Area  Code) 
(334)  255-6907 


22c.  OFFICE  SYMBOL 
MCMR-UAX-SS 


DDForm  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 
Unclassified 


19.  Abstract  (Continued) . 


repeated  impact  for  use  in  the  development  of  the  experimental  phase  and  in 
a  dose-effect  model.  Ten  subjects  participated  in  pilot  tests  using  the 
multiaxis  ride  simulator  at  Fort  Rucker.  Subjects  experienced  mechanical 
shocks  of  0.5  to  3.0  g  magnitude  with  wave  forms  of  2  to  11  Hz  in  a  series 
of  short-term  exposures.  Four  of  these  subjects  also  were  subjected  to  long 
duration  (1  and  2  hour)  exposures  to  investigate  fatigue  mechanisms.  These 
experiments  involved  exposure  to  repeated  impacts  from  0.5  to  3.0  g  at  6  Hz . 
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Introduction 


The  overall  objective  of  the  project  is  to  develop  a 

•  dose-response  model  that  will  predict,  and  ultimately 
minimize,  the  risk  of  injury  to  a  soldier  when  exposed  to 
the  repeated  impact  environment  of  tactical  ground  vehicles. 
The  project  spans  five  years  and  five  phases.  In  the  first 

*  phase  a  review  of  literature  was  conducted.  This  phase 
concluded  with  a  list  of  potential  measures  or  indices  that 
might  be  sensitive  to  vibration  and  impact  and  could  be 
measured  in  the  pilot  phase.  Phase  two  ran  concurrently 
with  phase  one.  Phase  two  involved  receiving  and  analyzing 
tapes  of  data  from  USAARL  containing  vibration  measurements 
from  tactical  ground  vehicles  (TGVs) .  A  variety  of 
characterization  methods  for  data  containing  mechanical 
shocks  and  repeated  impacts  were  developed  and  programmed. 
These  methods  are  meant  to  be  more  sensitive  to  shocks  than 
previously  available  methods.  This  allowed  "typical" 
vibration  and  impact  environments  to  be  defined  based  upon 
the  TGV  data  tapes .  These  were  used  to  plan  and  produce 
motion  signatures  to  drive  the  multi-axis  ride  simulator  for 
Phase  3  pilot  tests. 

In  preparation  for  the  Phase  3  pilot  tests,  pre-pilot 
experiments  were  undertaken  at  B.C.  Research  using  a  single¬ 
axis  shaker  and  impact  device.  The  impact  device,  which  was 
designed  to  produce  vertical  impacts,  was  pneiomatically 
powered  and  computer  controlled  to  deliver  ±1  g  impacts  at 
random  time  intervals.  Many  of  the  experimental  techniques 
were  developed  and  refined  at  this  stage  in  preparation  for 
the  pilot  phase.  The  phase  3  pilot  tests  were  conducted 
using  the  multi-axis  ride  simulator  (MARS)  at  Fort  Rucker. 

A  niomber  of  biomechanical,  physiological  and  biochemical 
variables  were  investigated  in  a  series  of  short  duration 
(5.5  minutes)  and  long  duration  (1  and  2  hour)  experiments. 
This  phase  3  report  details  the  experimental  methods,  data 
analyses,  and  results  of  the  pilot  tests. 

»  Phase  four  of  this  project  is  the  full  experimental  phase 

to  be  conducted  at  the  MARS  facility.  It  is  anticipated 
that  the  full  experiments  will  examine  a  smaller  number  of 
indices,  but  a  larger  number  of  subjects.  The  duration  of 

^  exposure  will  likely  be  longer,  and  could  involve  repeated 

exposures  over  consecutive  days.  Phase  five  is  the  analysis 
and  model  development  phase.  The  final  outcome  of  phase 
five  will  be  recommendations  for  a  health  hazard  assessment 
index  sensitive  to  the  health  effects  of  shocks  and  repeated 
impacts . 
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Project  Objectives 


The  immediate  objective  of  the  pilot  phase  is  to  evaluate 
biomechanical,  physiological  and  biochemical  human  responses 
to  repeated  impact,  determined  from  a  review  of  literature 
to  be  the  most  promising  for  the  prediction  of  injury  risk,  • 

and  to  determine  which  of  those  responses  are  the  most 
appropriate  for  use  in  the  main  experiment.  To  this  end, 
the  pilot  tests  contained  both  short-term  exposures  (5.5 
minutes)  and  long  duration  exposures  (1  and  2 -hour) . 

The  objective  of  the  short-term  exposures  was  to  evaluate 
the  immediate  human  response  to  a  range  of  impact 
situations,  including:  varying  rise- times  (frequencies)  of 
shocks;  peak  acceleration  amplitudes;  number  of  shocks  per 
unit  time;  and  combinations  of  shock  and  rms  background 
vibration.  These  exposures  (described  in  Table  3)  were 
designed  to  provide  information  about  the  transmission 
characteristics  and  human  responses  to  single  shocks  in  each 
of  the  X,  y  and  z-axes. 

The  objectives  of  the  long  duration  exposures  were  to 
examine  the  potential  fatigue  and  recovery  effects  of 
repeated  impact;  to  provide  information  about  the  human 
response  to  varying  severity  of  exposure;  and  to  investigate 
how  response  may  change  during  prolonged  exposure.  Seven 
conditions  of  2 -hour  exposures  (four  with  z-axis  impacts 
(vertical  direction) ,  one  with  z-axis  RMS  vibration,  one 
with  X,  y  and  z  impacts  and  one  with  no  vibration  or 
impacts)  and  2  conditions  of  1-hour  exposures  (with  x  and  y- 
axis  impacts  respectively)  were  presented  to  each  subject 
(details  of  the  design  are  presented  in  Table  5  and  Figure 
3)  . 


Both  the  short  and  long  duration  exposures  were  designed 
to  provide  information  for  the  planning  of  the  full 
experimental  phases  to  be  conducted  during  1994.  From 
results  of  the  pilot  tests,  appropriate  shock  signatures 
will  be  designed,  and  the  experimental  measures  that  can 
best  be  correlated  with  the  motion  environment  will  be 
selected  for  use  in  the  full  experimentation  phase. 


Background 


Many  epidemiological  studies  have  been  conducted  of  heavy 
equipment  operators  from  industries  such  as  agriculture, 
construction,  mining,  forestry,  and  the  military  (Rosegger 
and  Rosegger,  1960;  Konda  et  al.,  1985;  Beevis  and  Forshaw, 
1985;  Boshuizen,  Bongers,  and  Hulshof,  1990;  and  Milby  and 
Spear,  1974)  .  Some  studies  focus  on  subjective  symptoms  of 
health  problems,  such  as  backache.  Other  studies 


2 


investigate  objective  findings  of  disease,  such  as  back 
disorders  (intervertebral  disc  herniation  and 
spondylolisthesis) ,  as  diagnosed  through  clinical  or 
radiological  findings.  Most  disorders  associated  with 
vibration  are  not  specific  to  vibration,  but  occur  generally 
in  the  population.  They  may  be  associated  with  other 
ergonomic  or  environmental  problems.  However,  sufficient 
evidence  is  available  to  conclude  that  long-term  exposure  to 
vibration  can  be  harmful  to  the  spine  and  possibly  other 
systems  of  the  body  (gastro-intestinal  and  cardio¬ 
respiratory)  .  Several  hypotheses  have  been  developed  to 
explain  the  etiology  of  back  disorders.  One  hypothesis 
suggests  vibration  alters  nutrition  of  the  disc  (Dupuis  and 
Zerlett,  1986) .  A  second  suggests  dynamic  loading  of  the 
intervertebral  joints  causes  fatigue  damage  to  the  annulus 
of  the  intervertebral  discs  (Sandover,  1981)  . 

There  is  no  single  disorder  linked  to  whole  body 
vibration  or  impact.  Instead,  there  is  some  acceptance  that 
vibration  accelerates  the  onset  and  progression  of  currently 
recognizable  syndromes,  rather  than  causing  specific 
pathologies  (Seidel  and  Heide,  1986;  Dupuis  and  Zerlett, 

1986;  Hulshof  and  van  Zanten,  1987) .  Although  some  studies 
include  exposure  to  vibration  and  repeated  impact,  there 
have  been  few  studies  where  repeated  impact  is  a  variable 
under  consideration. 

Vibration  and  impacts  have  been  shown  to  have  acute 
effects  on  a  number  of  different  systems  in  the  body, 
including  the  cardiovascular,  respiratory,  and 
gastrointestinal.  Numerous  reviews  are  available  (Guignard, 
1972;  Guignard,  1985;  Guignard,  1974;  Weaver,  1979;  Barnes, 
1987;  Ramsay  and  Beshir,  1981).  Physiological  effects  are 
related  to  two  potential  mechanisms:  the  movement  of  organs 
and  tissues  and  a  generalized  stress  response  related  to 
intensity  and  duration  of  vibration  exposure.  Many  of  the 
responses  to  vibration  are  attributed  to  stimulation  or 
over-activation  of  the  sympathetic  nervous  system.  This  can 
result  in  increased  concentrations  of  catecholamines  and 
vasoactive  metabolites  which  in  turn  cause  a  generalized 
stress  response.  An  increase  in  heart  rate,  cardiac  output, 
respiration  rate  and  oxygen  uptake  occurs  in  response  to 
whole  body  vibration  (WBV) .  In  some  cases,  peripheral 
vasoconstriction  has  also  been  reported  (Spaul,  Spear,  and 
Greenleaf ,  1986;  Abu-Lisan,  1979) .  Acute  pathological 
effects  of  vibration  and  shock  have  included  injury  to 
viscera,  lung  and  myocardium  (Guignard,  1972),  bleeding  in 
the  gastro-intestinal  system  (Sturges  et  al.,1974),  and 
occasionally,  hemorrhage  of  kidney  and  brain  (Guignard, 

1972)  .  The  majority  of  this  work  has  been  conducted  using 
animals . 

Vibration  and  impacts  have  been  shown  to  have  a  niunber  of 
different  effects  on  the  ECG  signal.  Changes  have  been  seen 
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in  the  R-R  interval  (Ullsperger,  Seidel,  and  Menzel,  1986; 
Harada,  Kondo  and  Kinura,  1990) ,  heart  rate  variability 
(Harstela  and  Pilirainen,  1985;  Auffret,  Demange,  and 
Vettes,  1974) ,  P-R  interval  (Abu-Lisan,  1979)  and  T-wave 
amplitudes  (Roman  et  al,,1968).  It  is  also  possible  that 
mechanical  vibration  of  the  intestines  will  increase 
motility,  or  movement  of  ingested  material  without 
appropriate  breakdown  or  absorption  taking  place.  A  number 
of  epidemiology  papers  have  suggested  that  hearing  loss  due 
to  noise  is  exacerbated  by  vibration  (Rehm  and  Wieth,  1984; 
Chernyuk  and  Tashker,  1989) . 

Literature  was  reviewed  to  help  find  a  biochemical  marker 
for  general  stress,  fatigue,  and  tissue  or  organ  damage  in 
response  to  vibration  and  shock.  Animal  studies  have  shown 
that  exposure  to  vibration  resulted  in  damage  to  heart, 
lung,  brain,  kidney,  gastro-intestinal  (GI)  tract,  liver, 
skeletal  muscle,  adrenal  glands,  and  reproductive  organs. 
Some  damage  was  detectable  in  blood  and  urine,  while  others 
required  histological  examination  of  tissue.  Biochemical 
measures  in  blood  and  urine  are  routinely  used  to  evaluate 
general  stress  in  humans.  While  inflammation  is  not  specific 
to  vibration,  an  inflammatory  response  has  been  linked  to 
Raynaud ' s  phenomenon  both  of  occupational  and  non- 
occupational  origin  (Langauer-Lewowicka,  1976) .  A  marker  of 
fatigue  is  particularly  important  in  this  study,  since 
increased  fatigue  may  impair  physical  and  psychological 
performance,  and  increase  recovery  time.  Peripheral  muscle 
fatigue  has  been  related  to  changes  in  carbohydrate 
metabolism  (lactate,  glucose) ,  protein  and  energy  metabolism 
(ammonia) ,  cortisol,  and  electrolyte  balance  (K+,  Mg2+, 

Ca2+)  (Roberts  and  Smith,  1989) .  Maintenance  of  blood 
glucose  is  important  in  occupational  settings  to  prevent 
hypoglycemia  which  interferes  with  task  performance. 


Electromyography  (EMG)  can  be  used  to  monitor  various 
aspects  of  muscle  function.'  The  response  of  paraspinal 
muscles  to  whole-body  vibration  has  been  studied  using  EMG 
to  assess  localized  muscle  fatigue  (Hansson,  Magnusson,  and 
Broman,  1991;  Hosea  et  al.,  1986;  Magnusson,  Hansson,  and 
Broman,  1988;  Robertson  and  Griffin,  1989;  and  Wilder, 
Frymoyer,  and  Pope,  1983),  phase  and  timing  relationships 
between  muscle  response  and  acceleration  (Hagena  et 
al.,1986;  Robertson,  1987;  Robertson  and  Griffin,  1989),  and 
to  estimate  compressive  loading  and  torque  about  the  spine 
(Marras  and  Mirka,  1991;  Seidel,  Bluethner  and  Hinz,  1986; 
and  Ortengren,  Andersson,  and  Nachemson,  1981) .  These 
parameters  are  of  interest  because  of  their  association  with 
stabilization  of  the  spine,  and  their  subsequent  association 
with  back  pain  and  injury  to  spinal  tissues.  Muscle  fatigue 
may  diminish  the  ability  of  muscle  to  adequately  compensate 
for  perturbing  forces,  while  out-of -phase  or  untimely  muscle 
response  can  contribute  to  postural  destabilization  and 
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increase  both  torque  and  compressive  loading  of  the  spine 
(Seroussi,  Wilder,  and  Pope,  1989;  Seroussi  et  al.,1987). 

When  the  human  body  is  subject  to  vibration  or  impact,  it 
demonstrates  a  dynamic  response.  The  displacement  of 
tissues  and  the  forces  transmitted  by  them  alters  as  a 
function  of  time.  A  useful  method  to  assess  the  potentially 
harmful  effects  is  to  measure  the  relative  displacements  and 
hence  stresses  of  different  regions  of  the  body  in  response 
^  to  vibration  amplitude  and  frequency.  Transmission  of 

acceleration  can  be  expressed  in  terms  of  a  transfer 
function  that  defines  the  relative  magnitude  and  phase 
relationship  of  the  output  acceleration  in  a  particular 
region  (for  example,  the  spine)  compared  with  the  input 
acceleration  (for  example,  at  the  seat) .  Knowledge  of 
acceleration  transfer  functions  provide  insight  into 
behaviour  of  the  body  sub-systems,  and  enables  assessment  of 
input  acceleration  levels  and  frequencies  where  a  particular 
tissue  is  more  likely  to  be  damaged.  A  substantial  body  of 
knowledge  has  been  reported  concerning  the  transmission  of 
vibration.  However,  little  is  known  about  the  repeated 
impact  environment  and  the  dynamic  response  of  individual 
body  segments  to  vibration  and  impact.  Niimerous  attempts 
have  been  made  to  model  the  biodynamic  characteristics  of 
the  human  body,  from  simple  mass  spring  models  (Payne,  1991) 
to  highly  complex  representations  of  the  hxaman  body 
containing  multiple  degrees  of  freedom  (Amirouche  and  Ider, 
1988)  .  A  well  developed  model  could  prove  to  be  an  ideal 
tool  for  assessing  the  health  effects  of  impacts  and 
vibration. 

All  of  the  biomechanical,  physiological  and  biochemical 
measures  and  analytical  procedures  that  were  proposed  in  the 
pilot  testing  protocol  have  been  conducted  previously  in 
research  applications.  Many  of  the  measures,  such  as  ECG 
and  EMG,  are  common  clinical  tools  and  several  people  on  the 
research  team  have  used  them  in  other  research  studies . 

Some  of  the  measures  were  investigated  in  small-scale  pre¬ 
pilot  tests  using  a  sinusoidal  shaker  with  an  impact  device 
at  B.C.  Research  (Pre-pilot  testing  protocol  August,  1992, 
approved  September,  1992) .  Most  of  the  measures  have  been 
reported  by  others  in  vibrating  environments  (Zagorski  et 
al.,  1976;  Hansson,  Magnusson,  and  Broman  1991;  Harada, 

Kondo  and  Kinura,  1990;  Kjellberg  and  Wikstrom,  1987; 
Robertson  and  Griffin  1989;  Seidel,  Bluethner  and  Hinz, 

1986;  Spaul,  Spear  and  Greenleaf,  1986;  Ullsperger,  Seidel 
and  Menzel,  1986) .  However,  few  studies  have  involved 
shocks  such  as  would  be  encountered  in  military  vehicles. 
Fewer  studies  still  have  involved  exposure  to  vibration  with 
repeated  shocks  in  a  controlled  laboratory  environment  to 
study  the  physiological,  biochemical  and  biomechanical 
responses  of  the  body  which  are  linked  to  health  problems. 
This  study  fills  an  important  gap  in  the  literature.  A  more 
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detailed  review  of  the  literature  may  be  found  in  the  Phase 
1  report  (Village  et  al.,  1992). 


Methods 


Multi-Axis  Ride  Simulator 

The  multi-axis  ride  simulator  (MARS)  at  USAARL  (Fort 
Rucker)  was  used  to  reproduce  vibration  and  shock  signatures 
from  computer-generated  signals.  The  MARS  consists  of  a 
shake  table,  two  large  hydraulic  pumps  in  parallel,  three 
hydraulic  actuators  each  having  a  3-stage  valve  system, 
three  fail-safe  valves,  one  multi-channel  servo  controller 
(Schenck/ Pegasus  5900)  and  one  DEC-PDPll  based  multi-channel 
iterative  transfer  function  compensation  computer  (ITFC) 
with  associated  anti-aliasing  filters,  analog-to-digital 
converters  and  digital-to-analog  converters.  The  MARS  is 
capable  of  up  to  600  lbs  load,  a  frequency  response  range  of 
2-40  Hz,  up  to  ±4  g  peak  acceleration  and  up  to  ±3.5  inches 
peak  displacement  (7  inches  peak  to  peak) .  The  vibration 
frequency  and  acceleration  amplitude  are  controlled  by  a 
computer  synthesized  displacement  command  signal.  The 
output  acceleration  signal  at  the  shake  table  is  corrected 
for  the  table  transfer  function  by  means  of  a  correction 
matrix  and  an  iterative  process  in  which  the  input  and 
output  signals  are  compared  for  quality  of  fit. 

To  ensure  the  safety  of  the  subject,  several  automatic 
safety  and  shut-down  features  are  built  into  the  MARS 
facility  to  ensure  that  vibration  and  impacts  do  not  exceed 
pre-determined  levels.  Additionally,  the  MARS  operator 
continuously  monitors  the  vibration  and  shock  exposure  from 
the  control  booth.  He  is  also  in  clear  view  of  the  subject 
and  in  voice  communication  through  headsets  with  other 
experimenters  and  the  subject  monitor.  A  dedicated  subject 
monitor  is  positioned  facing  the  subject  within  3  feet  to 
exclusively  monitor  their  well  being.  Should  the  subject 
wish  to  discontinue  the  exposure  for  any  reason  they  need 
only  speak  to  the  monitor  or  wave  to  the  controller  of  the 
MARS  in  the  control  booth  to  stop  the  motion.  The  shake 
table  controls  include  a  hand-operated  stop  button  for  the 
subject  monitor,  and  touch  sensitive  switches  placed  in  the 
area  of  the  MARS  controller  and  subject  monitor. 

The  computer  control  system  includes  an  automatic  shut 
down  in  case  of  control  system  failure.  Automatic  shut  down 
can  be  triggered  by  table  velocity  and  position  feedback, ■ 
feedback  failure,  or  power  failure.  The  end  stops  of  the 
shaker  table  are  also  buffered  mechanically  to  limit  any 
overshoot  deceleration  to  a  maximum  of  6  g.  Acceleration  at 
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the  tabletop  is  displayed  on  a  screen  and  monitored  in  the 
control  booth  by  the  operator. 

The  signals  were  input  to  the  shake  table  and  iterated  to 
ensure  a  close  match  between  input  and  output.  The 
•  corrected  signal  was  stored  in  the  DEC-PDPll  computer  system 

at  the  MARS  facility.  Signals  for  the  shake  table  were 
generated  on  an  IBM  compatible  personal  computer.  The  shock 
waveforms  were  based  on  data  from  tactical  ground  vehicles 
(TGVs)  analyzed  on  a  VAX  4000  at  B.C.  Research. 

A  solid  metal  seat  was  securely  mounted  on  the  shake 
table  and  a  bean-bag  cushion  taped  on  top  of  the  seatpan. 

The  seat  did  not  have  a  backrest,  since  it  was  determined 
from  communication  with  USAARL  personnel  that  most  drivers 
and  occupants  of  TGVs  do  not  utilize  a  backrest.  The  seat 
could  be  adjusted  such  that  the  subject's  feet  could  rest 
comfortably  on  the  MARS  table  with  the  knees  and  hips  at 
approximately  90°.  The  cushion  was  designed  to  distribute 
the  subject's  weight  without  altering  the  input  acceleration 
signal.  The  cushion  was  taped  firmly  to  the  metal  seatpan 
in  order  to  minimize  lateral  shear  effects  between  the 
seatpan  and  cushion. 

Instriimentation 

The  main  data  acquisition  and  analysis  system  was  a  VAX 
4000  computer  system  (Digital  Equipment  Corporation  and 
Alpine) .  Up  to  sixteen  channels  of  data  were  collected  on 
the  VAX,  and  the  channels  were  simultaneously  recorded  on 
VHS  tape  using  a  TEAC  XR-510  FM  recorder  as  back-up.  Data 
collection  was  performed  using  the  Generalized  Data 
Acquisition/Analysis  Programs  (GEDAP)  Software  (National 
Research  Council,  Ottawa,  Canada,  1990)  modified  for  our 
purposes.  Calibrations  for  relevant  channels  were  stored  in 
a  data  control  file  (known  as  a  port  file) ,  developed  for 
each  experimental  condition  of  short-term  x,  y  and  z,  and 
long-term  x,  y,  z,  and  xyz  combined.  The  port  file 
(containing  data  control  information)  was  automatically 
linked  to  the  appropriate  experimental  data  file  during  the 
analysis  process.  Two  PC  computers  (provided  by  USAARL) 
were  also  utilized;  one  for  collection  of  Optotrak 
positional  data,  and  one  for  collection  of  Synthetic  Work 
Environment  performance  data.  A  timer  was  used  that 
produced  a  voltage  pulse  every  5  seconds  to  later 
synchronize  Optotrak  data  with  data  collected  from  the  VAX 
4000.  The  timer  was  connected  to  both  an  Optotrak  Data 
Acquisition  Unit  (ODAU)  channel  and  channel  1  of  the  VAX 
4000  during  acquisition  of  exposure  data. 

Acceleration  at  the  spine  was  measured  by  Entran 
miniature  accelerometers  (weight  0.3  gm;  range  of  ±  10  g  or 
±25  g) .  Acceleration  at  the  seat  was  measured  with  three 
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single  axis  Entran  accelerometers  {±  25  g) .  The 
accelerometers  were  positioned  in  a  triaxial  accelerometer 
block  which  was  housed  within  a  moulded  flexible  epoxy 
seatpad  constructed  according  to  the  specifications  in  SAE 
(1974) .  A  calibration  box  was  constructed  to  secure  the 
seatpad  accelerometer  during  calibrations  and  a  second 
calibration  box  was  constructed  to  hold  the  nine  single-axis 
accelerometers.  Accelerometers  were  calibrated  dynamically 
at  USAARL  by  securing  the  boxes  to  the  MARS  shaker  and 
running  a  swept-sine  wave  signal.  Calibration  factors  were 
calculated  by  comparing  accelerometer  RMS  output  values  to  a 
factory  calibrated  piezo-electric  accelerometer  (P.C.B. 

Model  301A03)  using  a  power  supply  (482A05) . 

Internal  pressure  was  measured  by  a  specially  constructed 
rectal  pressure  probe.  The  probe  was  50  cm  in  length  and 
terminated  in  an  Entran  (model  EPB-140W-5S)  miniature 
pressure  transducer  (range  ±  5  psi) .  The  transducer  and 
wiring  were  covered  with  2  layers  of  heat  shrink  tubing  20 
cm  in  length  to  provide  a  suitable  degree  of  strength  and 
flexibility.  Latex  rubber  was  injected  around  the 
transducer  walls  and  wiring  junction  to  seal  and  protect  it. 
Tygon  plastic  tubing  was  placed  over  the  distal  end  of  the 
connection  wires  and  sealed  at  the  heat  shrink  interface. 

An  amphenol  electrical  connector  was  attached  at  the  distal 
end.  Wires  within  the  tubing  were  kept  slack  so  that  the 
tubing  provided  strain  relief. 

The  transducer  output  (volts)  was  initially  calibrated 
with  a  mercury  column  manometer  using  air  pressure  and  a 
sealed  pressure-ported  vessel  in  which  the  transducer  was 
installed.  To  provide  a  simple  field  calibration,  the  pump 
valve  and  pressure  gauge  of  a  blood  pressure  cuff 
(sphygmomanometer)  was  calibrated  using  the  mercury 
manometer  within  ±  2  mmHg  and  found  to  be  accurate.  At 
USAARL,  the  transducer  was  then  calibrated  regularly  with 
the  pressure  gauge  and  a  digital  volt  meter  connected  to  the 
amplified  output  of  the  transducer.  The  probe  was  sealed  in 
a  pressure-ported  vessel,  and  the  side  tap  of  the  vessel  was 
attached  to  the  blood  pressure  pump  valve.  Pressure  in  the 
vessel  was  changed  in  steps  and  readings  from  the  pressure 
gauge  and  voltmeter  were  recorded  manually.  Data  were  later 
entered  into  a  Microsoft  EXCEL  spreadsheet  and  a  regression 
was  performed  to  yield  the  calibration  factor. 

The  accelerometers,  load  cell,  and  pressure  transducers 
were  powered  by  a  Terrascience  14  channel  Signal 
Conditioning  unit  (Terrascience  Ltd. ,  Vancouver,  Canada) . 
Gains  and  filters  were  established  for  each  channel 
according  to  Table  1  below. 
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Tzdsle  1 


Amplifier  gains  and  filters 


Amplifier 

Channel 

Signal  Description 

Gain 

Filter 

(Low  pass,  Hz) 

1 

Rectal  Probe  #1 

200 

391 

2 

Rectal  Probe  #2 

200 

391 

3 

T3  z  Accelerometer 

200 

391 

4 

L4  z  Accelerometer 

500 

391 

5 

Sz  Accelerometer 

500 

391 

6 

T1  X  Accelerometer 

200 

391 

7 

L2  X  Accelerometer 

500 

391 

8 

Sx  Accelerometer 

500 

391 

9 

T2  y  Accelerometer 

200 

391 

10 

L3  y  Accelerometer 

500 

391 

11 

Sy  Accelerometer 

|500 

391 

12 

Load  Cell 

500 

50 

13 

Spare 

2 

1 

14 

Spare  Accelerometer 

500 

50 

Note:  T3  z;  T1  x;  Sx;  etc  denotes  location  and  measurement 
axis  of  accelerometer .  Detailed  information  is  provided  in 
Table  7. 


Thoracic  (RT-CH)  and  abdominal  (RT-AB)  circumference  was 
measured  using  Respitrace  bands  fitted  around  the  chest  and 
abdomen  of  the  subject  and  connected  to  a  Respitrace  monitor 
(Ambulatory  Monitoring  Inc.,  N.Y. ) .  Changes  in 
circumference  indicated  thoracic  or  abdominal  displacement. 

ECG  was  monitored  during  exposures  by  a  Hewlett  Packard 
(HP)  3 -lead  electrocardiograph  (Monitor  model  no.  783 OA  with 
patient  lead  HP  14056B  and  HP14358-A) .  The  pre-and  post¬ 
experiment  12-lead  ECG  was  from  Marquette  Electronics  Inc. 
(MAC15)  . 

EMG  was  measured  using  the  Telemg  Electromyograph  from 
Bioengineering  Technology  Systems  (BTS) ,  Milan,  Italy.  The 
Telemg  system  consists  of  sensors  with  lOOx  pre-amplifiers 
at  the  electrodes,  a  battery  powered  fibre-optic  converter 
box,  and  amplifiers  with  band-pass  filters.  During  maximum 
voluntary  contractions,  amplification  was  set  at  200x  (lOOx 
at  the  electrode  preamplifier  and  2x  at  the  amplifiers) . 

For  submaximal  calibrations  and  vibration  and  shock 
exposures,  the  gain  was  set  at  lOOOx  (lOOx  at  electrode  and 
lOx  at  amplifiers) .  The  low  pass  filter  (anti-aliasing)  of 
the  EMG  amplifier  was  set  to  400  Hz  and  the  high  pass  to  10 
Hz.  A  notch  filter  at  60  Hz  was  not  used,  since  the  60  Hz 
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background  noise  was  negligible.  For  EMG  calibrations  a 
force  transducer  (Maywood  Instruments  Ltd.,  Basingstoke, 
England,  Model  U4000  Load  Cell)  was  used  (range  ±  lOON) . 

The  force  transducer  has  a  linear  response  and  was 
calibrated  with  known  weights  suspended  vertically. 

Positional  data  was  collected  from  the  spine  using  the 
Optotrak  System  (Northern  Digital,  Version  0.9,  June,  1992). 
The  markers  used  with  the  Optotrak  system  consist  of  infra¬ 
red  emitting  diodes  (IREDS) .  These  markers  are  attached  by 
cable  to  a  strober  unit  that  can  accommodate  six  IREDS,  and 
is  connected  to  the  Optotrak  system.  The  IREDS  are  sensed 
by  an  Optotrak/3020  which  contains  three  cameras  mounted 
within  a  1.1  m  long  bar.  The  cameras  were  positioned 
approximately  2  m  behind  the  subject.  Error  of  measurement 
for  calibrated  displacements  are  of  the  order  of  0.05  mm  for 
objects  placed  3  m  from  the  camera  position.  The  Optotrak 
system  includes  the  Optotrak/3020 ,  a  system  unit,  the 
Optotrak  Data  Acquisition  Unit  (ODAU) ,  computer  adapter 
interface  card,  and  data  acquisition  and  analysis  software. 

Instrumentation  and  Layout 

A  block  diagram  of  the  instrumentation  is  provided  in 
Figure  1 . 
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Data  Acquisition 

Output  from  the  timer,  force  transducer  (during  EMG 
calibrations),  ECG,  EMG,  spine  and  seat  accelerometers, 
internal  pressure  transducer,  and  abdominal  and  thoracic 
Respitrace  bands  were  collected  to  the  VAX  computer.  All 
channels  were  sampled  at  a  rate  of  800  Hz  to  accommodate  the 
frequency  response  of  EMG  signals.  During  short-term 
experiments,  15  channels  of  data  were  collected  in  287 
second  blocks.  This  was  the  maximum  collection  time 
possible  by  the  VAX  4000  when  sampling  15  channels  at  800 
Hz.  In  the  long-term  experiments  16  channels  of  data  were 
collected  in  270  second  blocks.  Details  of  data  acquired  to 
the  VAX  4000  during  short-  and  long-term  exposures  are 
provided  in  Table  2.  The  shock  direction (s)  for  each 
exposure  are  indicated  at  the  foot  of  Table  2 . 
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Table  2 


Channel  configuration  VAX  4000  data  acquisition 


Channel 

Short-term  Experiments 

Long-term  Experiments 

1 

Timer 

Timer 

Timer 

Timer 

Timer 

Timer 

Timer 

2 

ECG 

ECG 

ECG 

ECG 

ECG 

ECG 

ECG 

3 

EMG 

EMG 

EMG 

EMG 

EMG 

EMG 

EMG 

4 

EMG 

EMG 

EMG 

EMG 

EMG 

EMG 

EMG 

5 

EMG 

EMG 

EMG 

EMG 

EMG 

EMG 

EMG 

6 

EMG 

EMG 

EMG 

EMG 

EMG 

EMG 

EMG 

7 

IP 

IP 

IP 

IP 

IP 

IP 

IP 

8 

T3  z 

T3  z 

T3  z 

T3  z 

T3  z 

T3  z 

T3  z 

9 

L4  z 

L4  z 

L4  z 

L4  z 

L4  z 

L4  z 

L4  z 

10 

Sz 

Sz 

Sz 

Sz 

Sz 

Sz 

Sz 

11 

T1  X 

T1  X 

T2  y 

T1  X 

T1  X 

T1  X 

T2  y 

12 

L2  X 

L2  X 

L3  y 

L2  X 

L2  X 

L2  X 

L3  y 

13 

Sx 

Sx 

Sy 

Sx 

Sx 

Sx 

Sy 

14 

RT  AB 

RT  AB 

RT  AB 

RT  AB 

RT  AB 

T2  y 

RT-AB 

15 

RT  CH 

RT  CH 

RT  CH 

RT  CH 

RT  CH 

L3  y 

RT-CH 

16 

— 

— 

EMG 

EMG 

Sy 

EMG 

Shock 

X 

z 

y 

X 

z 

xyz 

y 

Direction 

Note:  T3  z;  L2  x;  Sz;  etc.  denotes  location  and  measurement 
axis  of  accelerometers.  Detailed  information  is  provided  in 
Tabl e  7 . 

During  most  short-term  experiments,  Optotrak  positional 
data  were  collected  at  a  rate  of  100  Hz  from  5  markers 
attached  to  the  spine.  During  the  long-term  experiments,  an 
additional  marker  was  placed  on  the  seat  and  the  6  markers 
were  collected  at  80  Hz.  At  the  end  of  each  day  the 
Optotrak  data  files  were  downloaded  to  the  VAX  computer 
through  an  Ethernet  connection.  All  data  collected  by  the 
VAX  were  stored  to  digital  data  tapes  at  the  end  of  each 
day. 


Synthetic  Work  Environment  data  was  collected  by  a 
dedicated  PC  (as  explained  under  Experimental  Methods) . 
Subjective  responses  were  recorded  manually  (see 
Experimental  Methods  and  Appendix  B  for  details) . 

Resting  data  was  collected  for  150  seconds  prior  to  all 
vibration  exposures.  The  first  40  seconds  of  each  short 
duration  exposure  was  a  warm-up  period.  Data  collection  was 
timed  and  collection  was  started  to  obtain  the  last  287 
seconds  of  the  5.5  minute  exposure. 

During  the  1-hour  and  2-hour  experiments,  data  was 
collected  for  150  seconds  pre-exposure.  Once  the  motion 
started,  data  was  collected  for  the  first  4.5  minutes  of 
every  second  5.5  minute  signature.  Thus,  during  the  2  hour 
and  1-hour  experiments,  there  were  a  total  of  12  and  6  data 
collection  epochs,  respectively. 

Experimental  Procedures 

Selection  of  Subjects 

Twelve  male  subjects  were  recruited  from  U.S.  Army 
personnel  assigned  to  Fort  Rucker.  Subjects  were  recruited 
through  a  notice  published  in  the  USAARL  newsletter  and 
personal  communication.  Two  of  the  twelve  subjects  were 
eliminated,  based  on  selection  criteria,  leaving  ten  to 
participate.  Selection  criteria  are  detailed  in  the 
Experimental  Protocol  developed  for  the  Human  Use  Review 
Committee.  Subjects  were  selected  to  be  between  20  and  40 
years  of  age  and  within  one  standard  deviation  of  the  mean 
for  height  and  weight  based  on  standard  military  data.  All 
were  required  to  have  experience  with  motion,  either  TGVs  or 
air  transport.  Ten  subjects  participated  in  the  short 
duration  exposures  and  four  of  the  ten  subjects  continued  to 
participate  in  the  2-hour  and  1-hour  experiments. 

Procedures  prior  to  experiments 

Prior  to  participating  in  the  experimental  protocol,  all 
subjects  were  required  to: 


1.  Read  and  sign  the  "Volunteer  Agreement  Affidavit", 
"Volunteer  registry  data  sheet"  and  "Information 
for  Subjects"  which  describes  the  experimental 
protocol  and  all  associated  risks  or  benefits. 

2 .  Participate  in  a  medical  examination  conducted  by 
USAARL  physicians,  complete  a  medical  history 
questionnaire  and  satisfy  all  exclusion  criteria. 
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3.  Read  the  document  "Instructions  to  Subjects"  and 
comply  with  physical  activity,  caffeine  and 
smoking  restrictions  prior  to  experimentation. 

4 .  Attend  a  briefing  at  Fort  Rucker  to  hear  an 
explanation  of  the  protocol  and  instructions . 

5a.  Perform  maximal  and  submaximal  exertions  with  the 
paraspinal  muscles  for  calibration  of  EMG 
(procedures  described  under  EMG  methods) . 

5b.  Perform  a  sustained  maximal  isometric  contraction 
of  the  paraspinal  muscles  for  60  seconds. 

6.  Participate  in  a  trial  exposure  of  16.5  minutes  on 
the  MARS  (2  g  shocks  for  5. 5  minutes  in  each  of  x, 
y  and  z  directions)  to  familiarize  them  with  the 
motion  and  instrumentation.  Subjects  were 
instrumented  but  no  data  was  collected. 

7.  Perform  a  minimum  of  10  learning  and 
familiarization  sessions  with  the  synthetic  work 
environment . 

Design  of  short-term  experiments 

The  short-term  experiments  were  designed  to  determine  the 
human  response  to  impacts  of  varying  amplitude,  rise  time, 
shock  number  and  background  rms/shock  relationship.  A 
series  of  seven  5.5  minute  exposures  in  each  of  the  x,  y  and 
z-axes  was  designed  to  investigate  these  parameters  (see 
Table  3) .  In  each  experiment  (consisting  of  seven 
exposures),  subjects  were  exposed  to  shocks  in  one  axis. 
Subjects  returned  on  two  more  occasions  to  complete  the 
experimental  design  in  all  three  axes. 

Shock  profiles  for  the  7  short-term  experiments  are  shown 
graphically  in  the  Phase  2  report.  The  Phase  2  report  also 
details  the  methods  of  development  of  the  signals. 

Prior  to  exposure,  a  12 -lead  ECG  was  performed  on  each 
subject  while  in  a  prone  position.  Subjects  were  then 
instrumented  for  internal  pressure,  ECG  and  EMG,  as 
described  in  the  experimental  methods  section.  The  subject 
then  mounted  the  MARS  and  the  internal  pressure,  ECG  and  EMG 
leads  were  connected  to  the  appropriate  amplifiers. 
Calibration  of  EMG  signals  were  performed  using  the 
procedures  described  in  experimental  methods .  Following  EMG 
calibration,  Respi trace  bands  were  placed  around  the  chest 
and  abdomen  and  calibrated.  Accelerometers  and  IREDs  were 
attached  to  the  skin  over  the  spinous  processes  as  described 
in  experimental  methods .  To  obtain  information  on  the 
transfer  function  (natural  frequency  and  damping 
coefficient)  between  the  spinous  process  and  the 
accelerometer,  the  skin  adjacent  to  the  accelerometer  was 
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displaced  and  released,  and  the  resultant  acceleration  time 
history  recorded. 


Table  3 

Design  of  short-term  exposures 


Experiment  Nxamber 
and  Variable 

Background 
rms  Level 

Shock 
Frequency 
(rise  time) 

Shock 

Amplitude 

Duration 

1 :  Rise  Time 

0.05  g 

2  ,  4 ,  6 ,  8 , 

0.5  g 

5 . 5  minutes 

3  shocks /min 

broad  band 

11 

2 :  Rise  Time 

2,  4,  6,  8, 

1  g 

5 . 5  minutes 

3  shocks /min 

11 

3 :  Rise  Time 

4,  5,  6,  8, 

2  g 

5 . 5  minutes 

3  shocks /min 

11 

4 :  Rise  Time 

4,  5,  6,  8, 

3  g 

5 . 5  minutes 

3  shocks /min 

11 

5 :  Shock  Number 

0.05  g 

2  shocks /min 

broad  band 

6  Hz 

2  g 

5 . 5  minutes 

10  shocks /min 

6  Hz 

2  g 

30  shocks /min 

6  Hz 

2  g 

6 :  Backgro\ind 

0.1  g 

6  Hz 

2  g 

5 . 5  minutes 

rms/ shock 

0.17  g 

6  Hz 

2  g 

relationship 

0.25  g 

6  Hz 

2  g 

2  shocks /min. 

0.4  q 

6  Hz 

2  g 

7 :  Response  to 

0.3  g 

sinusoidal 

N/A 

79  seconds 

single  amplitude 

at  2  Hz  to 

swept  sine  wave 

17  Hz  at  1/4 

Hz  per 

second 

■■■■■■ 

Subjects  were  instructed  to  face  forward  and  sit  in  a 
comfortable  upright  posture.  They  were  then  exposed  to  a 
standardized  30  second  signature  consisting  of  1  g  shocks  at 
3  shocks /minute  presented  randomly  in  the  x,  y  and  z 
directions.  Subjects  answered  a  subjective  questionnaire 
pertaining  to  discomfort  caused  by  the  vibration  and  shocks 
(Appendix  B) .  At  the  end  of  the  30  seconds,  the 
instrumentation  and  data  recordings  were  checked  to  ensure 
all  channels  were  functional  and  free  of  artifact.  Resting 
data  was  then  collected  for  150  seconds. 

# 

Seven  short-term  exposures  were  separated  by  150  seconds 
of  data  collection  at  rest.  The  exposure  sequences  were 
ordered  in  a  balanced  design  in  terms  of  the  seven  exposures 
and  the  three  directions  of  exposure.  At  minute  1.0  and  5.0  ^ 

of  exposure,  subjects  answered  the  subjective  questionnaire. 

Subjects  had  approximately  5  minutes  of  recovery  between  the 
exposures  (including  the  period  of  resting  data  collection) . 

At  the  completion  of  the  short-term  experiments, 
instruments  were  detached  and  the  subject  dismounted  the 
table.  When  sensors  had  been  removed  and  the  skin  cleaned. 
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a  12 -lead  ECG  was  again  performed.  The  same  subject 
returned  on  two  more  occasions  to  complete  the  exposure 
sequence  in  the  remaining  two  axes.  A  summary  of  the 
experimental  measures  collected  is  presented  in  Table  4 . 


TzQsle  4 


Summarv  of  depend 

Lent  variables  for  each  experiment 

Short  Duration 
Experiments 

1-hour  Experiments 

2 -hour  Experiments 

ECG 

ECG 

ECG 

EMG 

EMG 

EMG 

Spine  Acceleration 

Spine  Acceleration 

Spine  Acceleration 

Seatpad  Acceleration 

Seatpad  Acceleration 

Seatpad  Acceleration 

Internal  Pressure 

Internal  Pressure 

Internal  Pressure 

Abdominal  &  thoracic 
displacement 

Abdominal  &  thoracic 
displacement 

Abdominal  &  thoracic 
displacement 

Optotrak 
(spine  motion) 

Optotrak 
(spine  motion) 

Optotrak 
(spine  motion) 

Blood,  urine  and 
feces 

Synthetic  work 
environment 

Synthetic  work 
environment 

Subjective  response 

Subjective  response 

Subjective  response 

Design  of  long-term  (2  hour)  experiments 


The  long-term  experiments  were  designed  to  evaluate 
fatigue  and  the  related  severity  of  different  combinations 
»  of  vibration  and  shock.  Four  subjects  experienced  seven 

different  exposure  conditions  of  2  hours  duration.  Details 
of  each  exposure  are  provided  in  Table  5.  As  the  objective 
was  to  evaluate  the  severity  of  different  shock  amplitudes 
and  shocks /min,  the  frequency  of  the  shock  waveform  was  set 
at  6  Hz  in  all  cases. 

The  first  condition  was  a  control  condition  with  no 
vibration  or  shock  exposure.  The  second  condition  contained 
background  rms  vibration,  but  no  shocks.  Conditions  3  to  6 
contained  various  combinations  of  shock  amplitudes  and 
shocks/min,  in  the  z-axis,  while  condition  7  contained  x,  y 
and  z  shocks . 
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Design  of  long-term  (  2  hour)  exposures 


Number 

z-axis 

rms 

x-axis 

inns 

y-axis 

rms 

Shock 

Amplitude 

Shock 

Number 

1. 

Control 

none 

none 

none 

none 

none 

2. 

Background  rms 

0.16  g 

0.05  g 

0.05  g 

none 

none 

3. 

shocks  +  rms 

0.05  g 

0.05  g 

0.05  g 

1  g  z-axis 

32/min. 

4. 

shocks  +  rms 

0.05  g 

0.05  g 

0.05  g 

2  g  z-axis 

2 /min. 

5. 

shocks  +  rms 

0.05  g 

0.05  g 

0.05  g 

3  g  z-axis 

0.4/ 

min. 

6. 

shocks  +  rms 

0.05  g 

0.05  g 

0.05  g 

2  g  z-axis 

32/min . 

7  . 

shocks  +  rms 

0.05  g 

0.05  g 

0.05  g 

1  g  shocks 
in  X,  y, 
and  z-axes 

32/min. 

random 

Graphical  representations  of  the  7  signals  for  the  long¬ 
term  experiments  are  found  in  the  Phase  2  report,  along  with 
details  of  methods  for  computer  generation  of  the  signals. 
Experiments  2,  3,  4,  5  and  7  were  designed  to  have  equal 
vibration  dose  values  (VDVs)  according  to  the  British 
Standards  (BS  6087;  1987) .  To  achieve  equal  VDV  values, 
higher  amplitudes  of  shocks  are  presented  less  frequently 
(e.g.  3  g  every  2.5  minutes)  than  lower  levels  of  shocks 
{e.g.  1  g  at  3 2 /minute) .  Experiment  6  containing  2  g  shocks 
at  32 /minute  is  designed  to  provide  two-times  the  VDV  of  the 
other  experiments .  The  time  sequence  of  the  experiment  is 
shown  in  Figure  3.  Long-term  experiments  differed  only  with 
respect  to  collection  of  recovery  data  as  indicated  in 
Figure  3 . 

As  in  the  short  duration  experiments,  a  resting  12 -lead 
ECG  was  performed  on  each  subject  prior  to  exposure.  In 
four  of  the  2 -hour  experiments,  (exposures  1,2,6  and  7; 

Table  5)  biological  specimens  (blood,  urine,  and  feces)  were 
collected  from  each  subject  before,  during  and  up  to  36- 
hours  after  each  experiment.  Subjects  also  returned  dietary 
and  bowel  movement  recall  questionnaires  prior  to 
experimentation.  An  IV  Specialist  inserted  a  heparin  lok 
and  approximately  20-30  minutes  later  pre-exposure  blood 
samples  were  drawn.  Subjects  were  then  instrumented  for 
internal  pressure,  ECG  and  EMG.  Subjects  mounted  the  MARS 
and  performed  EMG  calibration  procedures  as  described  in 
experimental  methods .  Respitrace  bands  were  positioned 
around  the  chest  and  abdomen  and  calibrated.  Accelerometers 
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and  Optotrak  markers  were  attached  to  the  skin  over  the 
spinous  processes. 


Exposure 


P  =  Synthetic  work  battery  test  performed 
B  =  Blood  sample  taken 

*  Recovery  period  is  either  0,  150  or  600 
seconds  in  length 


Figure  3.  Time  History  for  long  duration  (2  hour) 
experiments 


Subjects  were  exposed  to  the  same  standardized  30  second 
signature  as  in  the  short-term  experiments  and  answered  the 
questionnaire  on  subjective  discomfort  (see  Appendix  B) . 
Resting  data  were  then  collected  for  150  seconds. 

The  2-hour  exposure  was  composed  of  5.0  minute  signatures 
of  shock  and  vibration.  The  MARS  briefly  ramped  down  and  up 
again  between  each  signature,  with  a  ramp  time  of 
approximately  10  seconds  In  selected  exposures,  the  MARS 
was  briefly  held  stationary  during  the  ramp  down  after  1- 
hour  while  mid-point  blood  samples  were  drawn  through  the 
heparin  lok.  Subjects  performed  the  Synthetic  Work 
Environment  tasks  and  watched  a  video  on  a  monitor 
alternately  for  approximately  20  minute  periods  throughout 
the  2 -hour  exposure.  The  subjective  questionnaire  was 
answered  at  10  minute  intervals  throughout  the  exposure. 

The  questions  were  presented  during  the  first  minute  of 
alternate  5  minute  signatures.  When  performing  the  Synthetic 
Work  Environment,  subjects  also  responded  to  the  subjective 
questionnaire.  An  additional  subjective  rating  was  obtained 
of  interference  with  performance  due  to  the  exposure. 

At  the  end  of  the  exposure,  while  the  subject  was  still 
on  the  table,  end-point  blood  samples  were  drawn.  In 
approximately  33%  of  conditions,  subjects  rested  while  2.5 
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minutes  of  recovery  data  (ECG)  were  collected,  followed  by 
EMG  fatigue  trials  (described  in  experimental  methods) .  For 
the  remaining  conditions,  two  subjects  rested  on  the  MARS 
for  10  minutes  during  which  ECG  data  was  collected.  The 
other  two  subjects  performed  EMG  calibration  immediately 
following  exposure  and  no  resting  data  was  collected. 

Following  these  collections,  subjects  dismounted  the 
MARS,  sensors  were  removed  and  the  skin  was  cleaned.  A 
resting  ECG  (12-lead)  was  performed.  In  the  experiments 
where  biological  specimens  were  collected,  subjects  provided 
post-exposure  urine  and  feces  samples  to  the  technicians . 
Usually  the  urine  was  provided  before  the  subject  left  for 
the  day,  but  the  feces  sample  was  returned  at  a  later  time. 
Subjects  also  provided  a  blood  sample  1-hour  post-exposure, 
and  at  this  time  the  heparin  lok  was  removed.  Subjects 
returned  to  the  MARS  facility  approximately  12,  24  and  36- 
hours  post-exposure  for  subsequent  blood  samples  (14.5  ml.) 
to  be  taken  by  syringe. 

Each  of  the  four  subjects  undertook  all  of  the  7  2 -hour 
exposure  conditions.  Presentation  of  signatures  was  based 
on  a  balanced  order.  A  summary  of  the  experimental  measures 
collected  is  presented  in  Table  4. 

Design  of  one-hour  experiments 

The  same  4  subjects  also  participated  in  one-hour 
experiments .  These  experiments  were  designed  to  assess  the 
effects  of  shocks  in  the  y  and  x  directions.  Details  of  the 
background  vibration  and  shock  levels  for  the  1-hour 
exposures  are  outlined  in  Table  6. 

Subjects  performed  each  one-hour  exposure  on  different 
days,  and  the  order  of  presentation  was  balanced  among  the 
four  subjects.  The  same  procedures  were  used  for  the  1-hour 
experiments  as  with  the  2 -hour;  however,  there  was  no 
collection  of  blood,  urine  and  feces  (Table  4) . 
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Telble  6 

Design  of  l-hour  exposures 


Number 

z-axis 
rms  1 

x-axis 

rms 

y-axis 

rms 

Shock 

Amplitude 

Shock 

Number 

8 .  shocks  &  rms 

0.05  g 

0.05  g 

0.05  g 

1  g  in 
y-axis 

3 2 /min 

9 .  shocks  &  rms 

0.05  g 

0.05  g 

0.05  g 

1  g  in 
x-axis 

32 /min 

*  Experimental  Methods 

Acceleration 

Acceleration  was  measured  at  the  seatpad,  and  at  the  skin 
surface  over  the  spinous  processes  of  the  lumbar  and 
thoracic  spine  at  the  vertebrae  indicated  in  figures  4a  and 
4b.  Miniature  accelerometers  were  attached  to  the  skin  by- 
small  squares  (<  1  cm^)  of  two-sided  carpet  tape  and  small 
strips  of  tape  were  placed  over  the  accelerometers  for 
further  security.  Due  to  limitations  of  data  acquisition, 
only  two  acceleration  axes  were  collected,  depending  upon 
the  direction  of  the  motion,  as  described  in  Table  7.  In 
the  experiment  involving  x,  y  and  z  shocks,  all  three 
acceleration  axes  were  collected,  but  no  Respitrace  or 
abdominal  EMG  data  were  obtained. 

Skin  transfer  functions  were  obtained  for  the  y  and  z 
accelerometers  while  the  subject  was  sitting  on  the  MARS 
table.  The  y-axis  accelerometer  was  perturbed  by  placing  a 
finger  on  the  skin  beside  the  accelerometer,  pulling  the 
skin  sideways,  and  then  releasing.  The  z-axis  accelerometer 
was  similarly  perturbed  by  pulling  the  skin  downwards  and 
releasing.  During  a  20  second  data  collection,  several 
perturbations  of  both  the  thoracic  and  the  lumbar 
accelerometers  were  performed. 

Triaxial  acceleration  was  also  measured  at  the  seat  by 
three  accelerometers  fitted  into  a  seatpad  that  was 
positioned  on  the  seat  between  the  subject  and  the  cushion, 
and  firmly  secured  by  tape.  As  seat  accelerations  were 
»  measured  at  the  seatpad,  any  modification  of  the  MARS 

acceleration  waveform  due  to  either  vibration  of  the  metal 
seat  structure  or  damping  effects  of  the  bean  bag  cushion, 
would  not  affect  the  accuracy  of  the  results  reported. 

Internal  pressure 

Internal  pressure  was  monitored  using  a  rectal  probe 
containing  a  miniature  pressure  transducer.  The  probe  was 
inserted  by  the  subject  to  a  depth  of  15  centimeters  beyond 
the  anal  sphincter,  as  marked  by  tape  on  the  probe.  After 
the  experiment,  the  probe  was  removed,  washed  with  detergent 
and  rinsed  in  clean  water  by  the  subject.  The  probe  was 
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Table  7 

Direction  of  shocks  at  the  seat  and  acceleration  axes 
_ collected. _ 


Anatomical 

Location 

Measure 

Orientation  of 
Accelerometer 

X  Shocks 

y  Shocks 

z  Shocks 

xyz  Shocks 

T1 

X 

forward 

yes 

no 

yes 

yes 

T2 

y 

left 

no 

yes 

no 

yes 

T3 

z 

up 

yes 

yes 

yes 

yes 

L2 

X 

forward 

yes 

no 

yes 

yes 

L3 

y 

left 

no 

yes 

no 

yes 

L4 

z 

up 

yes 

yes 

yes 

yes 

Abdominal  and  thoracic  displacement 

Abdominal  and  thoracic  circijmf erential  displacement  was 
monitored  during  vibration  exposure  using  an  inductive  coil 
sewn  into  belts  worn  at  the  chest  and  abdomen  (Respitrace, 
Ambulatory  Monitoring  Inc.).  Due  to  shape  and  length 
dependence  of  the  bands,  the  Respitrace  was  calibrated  with 
each  subject  for  each  trial.  Subjects  were  asked  to  provide 
a  maximal  inhalation  and  exhalation,  as  well  as  a  maximal 
and  minimal  abdominal  distention  to  provide  a  two  point 
calibration  of  the  chest  and  abdominal  belts,  respectively. 
During  calibration,  body  circiimference  at  the  level  of  the 
belts  was  measured  at  minimum  and  maximum  excursions  using  a 
tape  measure.  The  maximum  and  minimum  voltage  output  was 
read  simultaneously  from  the  Respitrace  digital  display. 
Sensitivity  of  the  system  in  mm/V  was  calculated  from  the 
data  and  this  calibration  factor  entered  into  the  GEDAP  data 
control  file. 

Electrocardiography 

Electrocardiography  (ECG)  involved  taking  a  12-lead  EGG 
prior  to  and  following  vibration  exposure  and  a  single  lead 
ECG  (3  electrodes)  during  vibration  exposure.  Twelve  lead 
ECG  records  were  obtained  with  the  subject  in  a  prone 
position.  For  the  single  lead  ECG,  the  skin  was  shaved 
where  necessary,  cleaned  with  an  alcohol  solution  and 
disposable  electrodes  (Graphic  Controls,  Ag/AgCl  electrodes 
F23T-30)  were  applied  to  the  left  and  right  thorax.  The 
reference  electrode  was  placed  on  the  left  scapula. 
Positioning  of  the  electrodes  is  shown  in  Figure  4 .  ECG  was 
monitored  at  rest  prior  to  vibration  exposure,  throughout 
the  vibration  exposure,  and  during  recovery  following  the 
long-term  exposures. 

Electromyography 

Electromyography  (EMG)  was  performed  using  disposable  0.5 
cm  diameter  Ag/AgCl  surface  ECG  electrodes  (F23-T,  Graphic 
Controls) .  The  skin  was  prepared  by  shaving  hair,  where 


necessary,  and  rubbing  the  skin  with  an  alcohol  pad. 
Electrodes  were  placed  one  centimeter  apart  bilaterally  over 
the  belly  of  the  erector  spinae  muscles  at  the  level  of  L3 
and  T9  (see  Figure  4) .  During  long-term  experiments  and 
short-term  experiments  in  z-  and  x-axes,  abdominal  EMG  was 
measured  from  the  second  band  of  rectus  abdominus  on  the 
right  side  approximately  5  cm  below  the  rib  cage.  The  EMG 
reference  electrode  was  placed  on  the  lower  ribs  of  the 
subject  on  the  right  side  of  the  back  and  slightly  posterior 
of  the  saggital  midline.  EMG  was  concurrently  recorded  from 
four  sites  during  the  short-term  experiments  and  from  five 
sites  during  the  long-term  experiments  as  discussed  in  the 
experimental  procedures . 

Subjects  participated  in  calibration  of  the  erector 
spinae  EMG  and  muscle  fatigue  procedures  by  performing  back 
extensions  against  a  force  transducer.  Subjects  were 
secured  to  the  seat  by  a  seatbelt  around  the  hips  and  cargo 
straps  around  each  leg  at  the  groin  to  prevent  forward 
movements  on  the  seat.  A  force  transducer  was  connected  to 
a  vertical  bar  at  the  level  of  the  sternum  and  the  bar  was 
secured  to  the  MARS  table  approximately  3  feet  in  front  of 
the  subject.  The  vertical  distance  between  the  subject's  L3 
vertebra  and  sterniim  was  measured  and  recorded.  A  strap 
connected  the  load  cell  to  a  harness  worn  over  the  shoulders 
of  the  subject.  Calibration  contractions  required  the 
subject  to  attempt  to  extend  the  trunk  at  the  waist  against 
the  resistance  of  the  chest  harness  and  hip  belts.  Subjects 
were  asked  to  perform  three  brief  maximal  isometric 
voluntary  contractions  and  a  series  of  five  controlled 
submaximal  contractions  (0%,  5%,  10%,  20%  and  30%  MVC) .  A 
digital  voltmeter  was  connected  to  the  load  cell  and  held  in 
front  of  the  subject  to  provide  feedback  regarding  target 
force  levels . 


Muscle  fatigue  was  evaluated  by  a  30 -second  sustained 
isometric  contraction  at  50%  MVC,  performed  against  the 
resistance  of  the  shoulder  harness  and  load  cell. 

Positional  data 

Positional  data  from  the  spine  were  measured  using  the 
Optotrak  System  (Northern  Digital,  Version  0.9,  June  1992). 
Five  infrared  emitting  diodes  (IREDS)  were  taped  with  two- 
sided  tape  to  the  skin  of  the  subject  over  the  vertebral 
processes  (C7,  T4,  T8,  LI  and  L5)  as  shown  in  Figure  4.  The 
Optotrak  is  self  calibrating,  therefore  no  calibration 
routine  was  required.  During  long-term  exposures,  a  sixth 
marker  was  placed  on  the  seat  cushion. 
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S\ibj  active  response 

Subjective  responses  were  elicited  through  aural 
administration  of  a  seven-point  scale,  similar  to  those  used 
in  other  vibration  studies.  Subjects  had  been  familiarized 
with  the  scale  in  advance.  During  the  vibration  and  shock 
*  exposures,  they  were  asked  to  identify  their  perceived  level 

of  discomfort  on  the  7 -point  scale.  Separate  response 
ratings  were  obtained  for  vibration,  shocks,  tiredness,  and 
interference  with  the  performance  task. 

Synthetic  Work  Environment 

A  cognitive  performance  test  battery  was  performed  during 
the  long-term  (1  and  2-hour)  experiments.  The  battery 
consists  of  four  tasks  performed  simultaneously  on  a  single 
large  screen  monitor:  Sternberg  memory  task;  a  Stability 
tracking  task;  a  4 -column  addition  task;  and  an  auditory 
monitoring  task.  The  subject  was  asked  to  give  each  task 
equal  priority.  Subjects  used  a  trackball  and  integral 
buttons  to  complete  the  task.  To  provide  stabilization  of 
the  subject's  arm  and  the  trackball  during  vibration,  an  arm 
support  was  constructed.  The  trackball  and  buttons  were 
attached  to  the  support  using  Velcro,  and  the  subject's  arm 
was  strapped  into  the  support.  The  arm  support  was  then 
attached  by  Velcro  to  a  flat  board  strapped  to  the  subject’s 
leg.  The  Synthetic  Work  Environment  required  approximately 
15  minutes  to  complete. 

Biological  specimens 

Blood,  urine,  and  fecal  samples  were  collected  and 
analyzed  in  4  of  the  two  hour  exposure  conditions: 
experiment  1  (control  condition) ,  experiment  2  (random 
vibration  alone) ,  experiment  6  (2  g  shocks  in  z-axis  at 
3 2 /minute) ,  and  experiment  7  (1  g  shocks  in  x,  y  and  z-axes 
at  32 /min)  (Table  5) . 

Urine  and  fecal  samples  were  collected  by  each  subject 
pre-exposure  and  post-exposure  in  the  four  experiments . 
Subjects  were  given  sterile  containers  for  collection  of 

*  their  own  urine,  and  cards  for  collection  of  fecal  samples. 
Fecal  cards  were  brought  to  the  lab  by  the  subjects,  and 
urine  samples  were  collected  at  the  MARS  facility. 

*  Blood  was  sampled  at  seven  time  points  for  each 
experiment:  pre-exposure,  at  1-hour  (mid)  exposure,  peak- 
exposure  (end  of  experiment),  1-hour,  12-hours,  24-hours, 
and  36-hours  post-exposure.  Due  to  logistics,  some  samples 
were  not  collected  at  the  exact  time  intervals  indicated. 
Experiments  which  were  conducted  from  14:00  to  16:00  had  the 
1-hour  post-exposure  blood  sample  collected  at  17:00,  and 
the  next  sample  was  taken  13  hours  later  at  06:00  hours. 
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All  blood  samples  were  drawn  by  a  trained  IV  Specialist 
(910  or  Phlebotomist  (92B) .  An  indwelling  catheter  was 
inserted  in  an  antecubital  vein  and  kept  patent  with  a 
heparin  lok  solution.  The  first  four  blood  samples  were 
collected  from  the  catheter.  Blood  samples  were  drawn  from 
the  catheter  into  a  syringe  and  were  transferred  immediately 
to  a  series  of  vacutainer  tubes  for  analysis  of  specific 
metabolites.  Subsequent  blood  samples  were  collected  by 
standard  venipuncture  into  a  series  of  vacutainers .  A 
maximum  of  27.5  ml.  of  blood  was  taken  at  pre-exposure,  at 
1-hour  (mid)  exposure,  and  post-exposure.  A  maximum  of  14.5 
ml.  of  blood  was  taken  at  1-hour,  12-hours,  24-hours,  and 
36-hours  post-exposure.  The  maximum  blood  draw  was, 
therefore,  140.5  ml.  in  the  38  hour  period  of  each 
experiment.  A  minimum  of  three  days  separated  experiments 
where  biological  specimens  were  collected.  A  maximum  total 
of  562  ml.  of  blood  was  collected  from  each  subject  over  a 
three  week  period. 

Blood  samples  were  processed  at  the  MARS  facility  by  a 
team  of  trained  technicians.  Serxam  separator  tubes  (SST) 
were  allowed  to  clot  for  30  minutes,  centrifuged  for  15 
minutes,  and  the  serum  was  pipetted  into  a  storage 
container.  The  serum  was  refrigerated  until  analysis. 

Tubes  with  7.5%  EDTA  for  hematology  were  inverted  gently 
then  refrigerated.  Three  other  blood  samples  collected  in 
tubes  containing  specific  preservatives  (potassium  oxalate 
and  sodium  fluoride  for  lactate  analysis;  sodium  heparin  for 
ammonia  and  catecholamine  analysis;  buffered  sodium  citrate 
for  von  Willebrand's  factor  analysis)  were  inverted  gently, 
centrifuged  for  15  minutes  and  the  supernatant  pipetted  into 
storage  vials.  These  samples  were  frozen  on  dry  ice  until 
analysis . 

Blood  and  urine  samples  were  transported  each  day  by 
courier  to  Allied  Clinical  Laboratory  in  Columbus,  Georgia 
for  analysis.  Fecal  cards  were  analyzed  at  Lyster  Hospital 
at  Fort  Rucker  two  days  per  week.  Results  from  both  Lyster 
Hospital  and  Allied  Clinical  Laboratory  were  returned  to 
researchers  on  a  regular  basis. 

Blood  samples  were  analyzed  for  the  metabolites  listed  in 
Table  8.  Only  the  metabolites  indicated  with  an  asterisk 
(*)  were  analyzed  in  the  samples  collected  at  1-hour,  12- 
hours,  2 4 -hours  and  3 6 -hours  post-vibration.  Plasma  volxome 
was  calculated  using  the  method  of  Dill  and  Costill  (1974) 
to  determine  the  percent  change  in  plasma  voliome  during 
experiments .  Metabolites  measured  in  urine  samples  are 
listed  in  Table  9 . 

Fecal  samples  were  analyzed  for  the  presence  of  fecal 
hemoglobin,  which  could  indicate  injury  of  the 
gastrointestinal  system.  Subjects  were  also  asked  to 
complete  a  24  hour  dietary  recall  and  a  bowel  movement 
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questionnaire  for  the  day  preceding  and  following  the 
experiment.  This  record  was  collected  to  assist  in  the 
explanation  of  potential  irregularities  in  biochemical 
results . 

A  summary  of  the  method  used  in  each  biochemical 
analysis,  and  the  coefficient  of  variation  of  the  method 
reported  by  Allied  Clinical  Laboratories  is  shown  in 
Appendix  C.  In  addition  to  the  subject's  samples,  extra 
♦  blood  and  urine  samples  were  sent  to  Allied  Clinical 

Laboratories  to  verify  the  quality  control  of  their  methods 
on  two  occasions.  The  coefficient  of  variation  determined 
without  the  knowledge  of  Allied  Clinical  Laboratories  is 
also  shown  in  Appendix  C.  A  coefficient  of  variation  could 
not  be  calculated  for  qualitative  measures . 

Other  measures 

Gastrointestinal  Measures 

During  the  pre-pilot  experimentation  at  B.C.  Research, 
the  feasibility  of  acquiring  Electrogastrogram  (EGG)  data  as 
a  measure  of  stomach  motility  and  transit  time  was  examined. 
It  was  originally  planned  to  use  the  BTS  surface  EMG  system 
for  this  measure.  However,  the  system  has  a  high  pass 
filter  of  1  Hz  which  could  not  be  removed.  EGG  frequencies 
of  interest  are  well  below  1  Hz.  A  Physiometer  Model  400 
EMG  system  was  used  in  place  of  the  BTS.  Significant 
problems  were  experienced  with  respect  to  pre-conditioning 
and  filtering  of  the  data,  reliability  of  the  equipment,  and 
electrode  placement  with  short  pre-amplifier  leads .  Data 
collected  were  not  reliable  and  it.  appeared  that  the  EGG 
signal  itself  was  being  masked  by  body  sway.  On  further 
consultation  with  a  researcher  experienced  in  EGG 
measurement,  it  was  decided  that  collecting  EGG  during  the 
vibration  and  shock  period  was  not  viable  given  the  effect 
of  the  movement  on  the  signal.  This  measure  was,  therefore, 
abandoned . 

Other  methods  to  determine  gastrointestinal  transit  time 
k  were  assessed  for  feasibility  and  cost.  The  hydrogen  breath 

method  appeared  to  offer  significant  advantages  over  methods 
such  as  Surface  Vibration  Analysis  (Campbell,  Storey  et  al., 
1989)  or  Radionucleotide  imaging.  However,  we  were  unable 
^  to  rent  or  borrow  a  hydrogen  analyzer  and  purchase  of  this 

equipment  was  not  included  in  the  budget.  This  measure  was, 
therefore,  abandoned. 
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Table  8 


Metedaolites  measured  from  blood  samples 


Metabolite  Name 

Relevance  to  Study 

Alkaline  Phosphatase 

Bone  remodeling 

Ammonia 

Fatigue 

Bilirubin 

RBC  damage 

Blood  Urea  Nitrogen 

Kidney/ liver  function 

Calcium 

Electrolyte  shift 

Catecholamines 

Stress/adrenal  function 

Cortisol* 

Stress /Adrenal  Function 

CPK* 

Skeletal  muscle  damage 

CPK- ( iso-enzyme  profile)* 

Skeletal  vs  cardiac  muscle 
damage 

Creatinine 

Renal  function 

Free  Hemoglobin 

RBC  damage 

Glucose 

Hypoglycemia 

Haptoglobin 

RBC  damage 

Hematocrit* 

Fluid  shift 

Hemoglobin* 

Fluid  shift 

Lactate 

Fatigue 

LDH* 

Skeletal  muscle  damage 

LDH- (iso-enzyme  profile)* 

Skeletal  vs  cardiac  muscle 
damage 

Magnesium 

Electrolyte  shift 

Myoglobin 

Skeletal  muscle  damage 

Potassium 

Electrolyte  shift 

Sodium 

Electrolyte  shift 

Total  protein 

tluid/protein  shift 

Total  vs  Conjugated  Bilirubin 

RBC  vs.  liver  dysfunction 

Uric  Acid 

O  O 

1  or  2  Hyperuricemia 

von  Willebrand's  Factor 
Antigen* 

Capillary  endothelial  damage 

White  blood  cell  profile* 

Inf lammat ion/ infect ion 

1 
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Table  9 


Metabolites  measured  from  urine  s^^wples 


Metabolite  Name 

Relevance  to  Study 

Bacteria 

Infection 

Casts 

Urinary  tract  injury 

Crystals 

Kidney  stone  formation 

Hemoglobin 

RBC  damage 

Hydr oxypr o line 

Joint  damage 

Mucous 

Infection 

Myoglobin 

Skeletal  muscle  damage 

Protein 

Kidney  dysfunction 

RBCs 

Urinary  tract  bleeding 

Specific  Gravity 

Urinary  concentration 

WBCs 

Inflammation 

Vasoconstriction 

During  the  pre-pilot  phase,  a  nuinber  of  clinical 
technicians  and  thermal  physiologists  were  contacted  about 
the  best  way  to  assess  vasoconstriction.  The  following 
equipment  and  methods  were  examined  or  investigated; 
impedance  plethysmography,  photo-electric  plethysmography, 
laser-doppler  flow  meter,  ultrasonic  doppler  flow  detector, 
and  heat  flow  sensors  with  integral  thermister.  The 
majority  of  the  methods  could  not  be  used  to  reliably 
measure  blood  flow  while  the  subject  was  moving.  Large 
artifacts  appeared  in  the  data.  Since  vasoconstriction  was 
likely  to  disappear  on  cessation  of  the  stressor,  a  resting 
and  recovery  measurement  approach  would  not  have  been 
appropriate.  The  remaining  methods  did  not  directly  measure 
blood  flow  and  results  had  to  be  extrapolated.  Given  the 
lack  of  tight  environmental  control  at  the  MARS  facility  and 
the  problems  of  acquiring  artifact-free  data,  this  measure 
was  abandoned . 

Temporary  Threshold  Shift 

The  literature  survey  revealed  that  exposure  to  vibration 
as  well  as  noise  may  cause  greater  hearing  loss  (permanent 
and  temporary)  than  that  caused  by  noise  alone.  One  method 
of  evaluating  this  effect  is  by  measuring  Temporary 
Threshold  Shift  (TTS) .  The  inclusion  of  noise  as  a  stressor 
in  the  experiment  and  the  involvement  of  USAARL  staff  and 
facilities  to  measure  the  TTS  was  discussed  with  the 
Scientific  Authority  and  staff  at  USAARL.  Although  the 
topic  is  of  inherent  interest,  the  staff  responsible  for 
noise  studies  did  not  have  the  resources  to  devote  to  this 


work.  As  this  work  was  not  in  the  original  proposal,  and 
therefore  not  budgeted,  the  measure  was  discontinued. 


Data  Analysis 


Data  acquired  during  the  pilot  experiments  at  USAARL  were 
downloaded  from  the  VAX  4000  to  digital  tapes  each  day. 
Subsequent  data  analysis  of  the  pilot  experiments  took  place 
at  B.C.  Research.  As  the  available  disc  space  on  the  VAX 
4000  was  insufficient  to  accommodate  all  data  files,  the 
data  were  selectively  reloaded  to  the  VAX  for  analysis.  The 
data  were  organized  within  two  separate  directories  for 
short  duration  experiments  (5.5  min  exposures)  and  long 
duration  experiments  (1  and  2  hour  exposures)  with  each 
directory  stored  on  separate  discs  of  approximately  1 
gigabyte  capacity.  Short  duration  data  were  organized 
within  3  subdirectories  according  to  the  direction  of  the 
shocks  (short_x;  short_y;  and  short_z) .  The  long  duration 
data  were  organized  within  nine  subdirectories  according  to 
exposure  condition  (Ll  to  L9) . 

Data  analysis  programs  were  written  to  investigate  the 
effects  of  individual  shocks  on  ECG,  EMG,  spinal 
accelerations,  internal  pressure,  and  chest  and  abdominal 
displacement.  Programs  were  also  written  to  investigate 
cumulative  effects  of  repetitive  shocks  on  ECG  (T  wave 
amplitude,  R-R  interval,  frequency  spectrum  of  R-R  interval, 
heart  rate)  and  EMG  activity  (mean  power  frequency  and  lEMG 
time  series) .  It  was  hypothesized  that  fatigue  effects  due 
to  repetitive  shocks  would  cause  a  change  in  the  parameters 
of  these  signals.  Indices  of  fatigue  and  tissue  damage  were 
also  sought  in  the  biochemical  data. 

Transmission 


Acceleration  data 

A  number  of  seat  and  spinal  acceleration  records  were 
inspected  using  the  VAX  4000  GEDAP  software.  Individual 
shock  events  in  the  x,  y  and  z-axes  were  selected  manually 
from  the  graphical  display,  expanded  in  the  time  domain 
using  a  ZOOM  feature,  then  plotted  to  hard  copy.  The  seat 
acceleration  waveforms  contained  higher  frequency  components 
superimposed  on  the  original  shock  signal  computed  as  input 
to  the  MARS  controller.  A  number  of  software  filters  were 
evaluated  to  improve  signal  quality.  A  low  pass  filter  with 
a  cut-off  frequency  of  60  Hz  removed  most  of  the  high 
frequency  components  without  attenuating  the  peak  of  the 
underlying  shock  waveform.  Higher  frequency  components  of 
acceleration  generally  were  not  observed  in  the  spinal 
acceleration  data,  with  the  exception  of  some  z-axis 
negative  2  g  and  negative  3  g  shocks.  However,  to  maintain 


consistency  in  the  treatment  of  acceleration  data,  and  avoid 
distortion  of  transmission  ratios,  the  spinal  acceleration 
records  were  also  low-pass  filtered  at  60  Hz.  The 
acceleration  data  were  also  high-pass  filtered  at  0.5  Hz  to 
remove  any  baseline  bias  or  zero  drift  from  the  data. 

Details  of  the  choice  of  filter  cut  off  and  the  effect  of 
filtering  on  the  underlying  signal  are  presented  in  the 
results . 

Vertebra-skin  transfer  function 


Prior  to  analysis  of  shock  transmission,  the  bone-to-skin 
transfer  function  of  accelerometers  measuring  spinal 
response  in  the  z  and  y  directions  was  calculated  for  each 
subject.  The  acceleration  response  of  the  skin  to  a  single 
perturbation  (as  described  in  methods)  was  viewed  on  the 
display  terminal  of  the  VAX  4000.  The  GEDAP  software  was 
used  to  digitize  the  magnitude  and  time  of  adjacent 
acceleration  peaks  within  the  signal.  By  considering  the 
accelerometer  and  subcutaneous  tissue  to  act  as  a  simple 
Kelvin  element,  the  damping  and  natural  frequency  of  the 
accelerometer-tissue  subsystem  was  calculated  from  the 
logarithmic  decrement  in  amplitude,  and  the  period  of  the 
waveform,  as  described  by  Hinz  et  al.,  1988. 


and  GJ_  = 


where,  6=  logarithmic  decrement 


1=  period  (s) 

fraction  of  critical  damping 
G5^  =  angular  natural  frequency  (rads/s) 

The  transfer  function  between  the  spinous  process  and 
accelerometer  is  characterized  by  the  amplitude  ratio 


H(03^  ,0  = 


1  +  (2C  TO  /  (g^) 


1(1-  (CJ  /  +(2^(3/  GJ^)^ 


and  the  phase  angle. 


0  (C3„  ,t^)  =  TAN' 


-1 


2C(ro  / 


(1-  (ta  /  +(2Cn3  / 


where,  CJ  =  angular  frequency 

())  =  phase  angle. 
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The  transfer  function  of  each  accelerometer  for  each 
subject  was  based  on  the  average  values  of  ^  and  obtained 
from  a  number  of  perturbations. 

To  estimate  the  acceleration  response  of  the  vertebra 
underlying  the  accelerometer,  the  spinal  acceleration  data 
was  converted  from  the  time  domain  to  the  frequency  domain. 
The  frequency  spectrum  was  multiplied  by  the  inverse  of  the 
transfer  function,  and  the  data  then  reconstructed  in  the 
time  domain.  This  mathematical  treatment  of  the  data 
provided  an  estimate  of  the  input  acceleration  signal  at  the 
spinous  process  necessary  to  produce  the  output  acceleration 
signal  measured  at  the  skin  surface. 

A  computer  program  was  written  to  apply  an  inverse 
transfer  function  to  each  y-and  z-axis  acceleration  record 
measured  at  the  spine.  The  program  was  tested  on  several 
sample  records  of  spinal  shocks,  and  the  output  data 
visually  inspected.  The  acceleration  data  were  filtered  at 
60  Hz  prior  to  application  of  the  inverse  transfer  function. 
It  was  found  that  this  cut-off  frequency  smoothed  the  low 
amplitude,  higher  frequency  components  within  the  signal 
without  disturbing  the  underlying  shape  of  the  response. 
Details  of  skin  transfer  functions,  and  their  effect  on  the 
accelerations  measured  at  the  spine  in  response  to  shocks 
are  presented  in  Results  and  Discussion  sections. 

Internal  pressure 

Visual  inspection  of  the  internal  pressure  data  showed  a 
characteristic  pressure  response  to  shock  inputs  at  the 
seat,  superimposed  on  slower  fluctuations  of  internal 
pressure  with  time.  To  analyze  the  shock  response,  the 
pressure  data  were  high  pass  filtered  at  0.5  Hz  to  remove 
these  fluctuations  and  provide  a  stable  baseline.  The  data 
were  also  low  pass  filtered  at  60  Hz  to  conform  with  the 
seat  acceleration  data.  Inspection  of  typical  shock 
responses  showed  that  the  low  pass  filter  generally  did  not 
affect  the  waveform,  with  frequency  components  of  the 
internal  pressure  being  well  below  this  value. 

Abdominal  and  thoracic  displacement  (Respi trace) 

Inspection  of  Respitrace  records  showed  a  shock  response 
superimposed  on  the  respiratory  pattern.  The  data  were  high 
pass  filtered  at  1  Hz  to  remove  the  oscillations  of 
respiration,  and  low  pass  filtered  at  60  Hz  to  conform  with 
the  treatment  of  acceleration  data.  The  Respitrace  data  did 
not  display  the  higher  frequencies  noted  in  the  acceleration 
records  and  hence  the  low  pass  filter  did  not  affect  the 
shape  of  the  Respitrace  waveform. 
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Shock  peak  detection  and  transmission  ratio 

To  determine  the  transmission  of  acceleration  from  the 
seat  to  the  spine,  a  data  analysis  program  was  written  to 
identify  the  discrete  shock  events  within  the  short-term 
acceleration  data.  The  program  allowed  the  user  to  input 
the  number  of  shocks  to  be  identified  at  the  seat,  the 
spinal  acceleration  (input  data  file)  to  be  compared  to  the 
accelerations  at  the  seat,  and  the  time  period  (window)  of 
each  shock  event  to  be  analyzed.  The  program  then 
identified  the  highest  acceleration  values  occurring  within 
the  seat  acceleration  data  file.  To  prevent  two  maxima 
being  identified  within  a  single  shock  (due  to  higher 
frequencies  superimposed  on  the  shock  waveform) ,  an 
exclusion  window  was  set  around  the  shock  peak  once  the 
(highest)  maxima  had  been  identified.  As  the  timing  of  the 
shocks  within  the  MARS  control  signature  was  known,  the 
timing  of  each  shock  peak  detected  at  the  seat  was  provided 
as  output  from  the  program  to  confirm  that  the  correct  peaks 
had  been  identified. 

Once  the  positive  peak  value  of  seat  acceleration  was 
identified,  the  program  identified  the  acceleration  minima 
(negative  peak)  occurring  within  the  shock  window.  The 
program  then  analyzed  the  spinal  response  data  to  identify 
the  corresponding  positive  acceleration  peak  and  negative 
acceleration  peak  occurring  within  each  shock  window,  and 
the  time  delay  between  the  seat  and  spinal  acceleration. 

Analysis  and  visual  inspection  of  a  number  of  data  files 
established  that  a  shock  window  of  200  msecs  was  adequate  to 
provide  identification  of  the  positive  and  negative 
acceleration  peaks  in  response  to  the  shock  input  at  the 
seat.  As  the  0.5  g  and  1.0  g  data  contained  2  Hz  shocks  it 
was  found  that  a  shock  window  of  400  msec  was  required  to 
accommodate  the  slower  waveform.  The  output  from  the 
program  provided  the  time  of  each  shock  event;  initial  peak 
and  peak-to-peak  acceleration  at  the  seat;  initial  peak  and 
peak-to-peak  acceleration  at  the  spine;  seat  to  spine 
transmission  ratio  for  the  initial  and  peak-to-peak 
accelerations;  and  time  delay  between  the  initial  seat  and 
spinal  acceleration  peaks.  Data  output  was  presented  in  two 
files;  a  GEDAP  table  and  an  ASCII  file  that  could  readily  be 
exported  to  Microsoft  EXCEL. 

The  same  peak  detection  program  was  used  to  analyze  shock 
transmission  for  internal  pressure  and  Respitrace  data.  In 
these  cases,  transmission  was  expressed  as  the  ratio  of  the 
peak  output  to  acceleration  input  (in  units  of 
inmHg.m-i.s2  and  mm.m“i.s2  respectively).  A  window  of  250 
msec  was  used  in  the  case  of  the  internal  pressure,  and  400 
msec  in  the  case  of  the  Respitrace  data. 
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Data  processing 

A  batch  program  was  written  to  process  the  transmission 
data  of  all  ten  subjects  for  the  short-term  exposures  to  0.5 
g,  1.0  g,  2.0  g  and  3.0  g  shocks.  The  batch  program  first 
demultiplexed  the  data  files  of  one  subject  at  one  exposure 
condition.  The  acceleration,  internal  pressure  and 
Respitrace  data  files  were  then  treated  by  GEDAP  software  to 
remove  the  mean.  This  eliminated  any  offset  error  in  the 
transducer  outputs.  The  data  files  were  filtered  as 
described  in  the  previous  section  and  the  spinal 
accelerations  treated  with  the  appropriate  inverse  transfer 
function.  The  peak  detection  program  was  then  used  to 
compare  the  spinal  accelerations,  internal  pressures  and 
Respitrace  displacements  to  the  input  seat  accelerations . 
Shock  transmission  ratios  were  calculated  for  each  shock 
event.  A  printout  of  tabulated  results  was  obtained  from 
the  peak  detection  program  for  each  response  measure. 

This  process  was  repeated  for  all  four  shock  conditions 
(amplitudes) ,  in  all  three  directions  of  acceleration  and 
for  all  ten  subjects.  In  EXCEL,  tables  were  created  to 
combine  transmission  data  and  delays  (between  peak  shock 
input  and  peak  response)  for  all  10  subjects  for  1,  2,  and  3 
g  shocks  for  z-axis  shocks:  lumbar  to  seat  (L4  z/Sz) ; 
thoracic  to  seat  (T3  z/Sz) ;  lumbar  to  seat  (L2  x/Sz) ; 
thoracic  to  seat  (T1  x/Sz) ;  internal  pressure  to  seat 
(IP/Sz) ;  abdominal  Respitrace  to  seat  z;  and  chest 
Respitrace  to  seat  z.  These  21  tables  were  repeated  for  x 
and  y  directions.  Also  in  EXCEL,  plots  were  developed  to 
show  each  of  the  10  subjects’  individual  response  to  the 
average  of  2  shocks  at  each  frequency  and  amplitude.  The 
mean  of  the  10  subjects  was  calculated  and  superimposed  on 
the  plots  as  points  at  the  appropriate  frequencies . 

Electrocardiography 

Analysis  of  the  ECG  was  aimed  at  answering  two  separate 
questions:  1)  what  instantaneous  effect  do  shocks  have  on 
the  ECG  R-R  interval;  and  2)  does  the  long-term  exposure  to 
shocks  cause  fatigue? 

Instantaneous  effect  of  shocks 

It  was  hypothesized  that  any  effect  of  shock  on  the  R-R 
interval  may  vary  depending  on  where  in  the  ECG  cycle  the 
shock  occurs.  For  example,  the  mass  of  the  heart  may  be  a 
significant  variable  with  respect  to  the  shock  response.  To 
test  this  hypothesis,  the  ECG  waveform  was  divided  into 
three  sections,  based  on  the  filling  and  emptying  cycle. 
These  were  categorized  as  full  (PQRS) ,  emptying  (T-wave)  and 
filling  (T-P)  interval  as  shown  in  Figure  5. 
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Figure  5 .  Demarcation  of  ECG  waveform  complexes  for  R-R 
analysis  and  the  associated  left  ventricular 
volxjme 


An  analysis  program  was  then  written  to  categorize  the 
ECG  signals  by  shock  position.  The  program  first  locates  the 
largest  peaks  in  the  ECG  waveform  (usually  the  R-wave  peak) 
and  then  travels  forwards  and  backwards  in  time  to  locate 
the  adjacent  Q  &  S-wave  depressions.  Since  the  P  &  T-wave 
peaks  are  susceptible  to  masking  by  shock  artifacts,  P  and 
T-wave  recognition  by  simple  peak  and  slope  detection 
strategies  proved  to  be  unworkable.  Instead,  a  simple  set  of 
ratios  was  used  to  approximate  the  location  of  the  beginning 
of  the  P-wave  and  the  ending  of  the  T-wave  based  on  the 
previous  R-R  interval .  The  interval  between  the  R  wave 
preceding  and  succeeding  the  shock  was  then  assessed, 
together  with  the  following  14  R-R  intervals.  Data  from  the 
long-term  experiment  4  (2  g  shocks,  2 /minute)  were  analyzed 
for  one  subject  from  the  following  periods  of  the 
experiment : 
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Rest 


Epochs  1-3  of  Exposure  (minutes  0-5 ; 10-15 ; 20-25 ) 

Epochs  10-12  of  Exposure  (minutes  90-95 ; 100-105 ; 110-115 ) 

Long-term  experiment  4  was  chosen  for  analysis  as  it  '*■ 

offered  a  balance  between  number  of  shocks  and  shock 
amplitude.  There  was  enough  time  after  the  shock  to  analyze 
the  resulting  R-R  intervals  before  the  next  shock  occurred. 

Fatigue  effects 

The  fatigue  effect  of  long-term  exposure  to  shocks  was 
addressed  by  analyzing  T-wave  amplitude-s,  R-R  interval, 
heart  rate,  and  changes  in  spectral  components  of  ECG 
waveform.  The  analysis  of  T-wave  amplitudes  was  initially 
undertaken  manually  where  the  peak  of  the  T-wave,  R-wave  and 
base  of  the  T-wave  were  displayed  on  the  computer  screen  and 
hand-digitized.  Corrections  were  made  for  changes  in 
baseline  and  for  changes  in  overall  signal  strength  (in  case 
impedance  had  varied) .  The  T-wave  amplitudes  were  therefore 
expressed  as  (T-wave  -  Baseline) / (R-wave  -  Baseline).  Part 
way  through  the  manual  analysis,  a  program  was  written  which 
automatically  detected  the  T-wave  using  the  same  principles 
as  in  the  waveform  analysis  described  above.  The  output  of 
the  program  was  transferred  into  EXCEL.  The  data  were  then 
graphed  for  visual  verification.  Means  and  standard 
deviations  were  calculated  of  the  corrected  T-wave 
amplitudes  for  the  following  segments  of  the  experiment: 

Rest 

Epoch  1 

Epoch  6 

Epoch  12 

Recovery  (where  available) 


The  programs  used  to  analyze  R-R  interval,  heart  rate  and 
spectral  components  (Variance  of  Baroreceptor  or  Traube- 
Hering-Mayer  wave  0.06-0.1  Hz;  Variance  of  Sinus  Arrythmia 
0.12-0.4  Hz)  firstly  identify  the  R-waves  using  a  threshold 
approach.  The  R-R  peak  intervals  are  then  calculated  and 
plotted  against  time.  Since  the  shocks  produce  movement 
artifacts  that  can  be  identified  as  R  peaks,  interval 
rejection  limits  were  built  into  the  program.  These  were 
initially  set  at  less  than  0.666  sec  and  greater  than  1.166 
sec.  R-R  intervals  outside  these  limits  were  rejected. 
Subject  1  (long-term  experiments)  had  a  very  high  heart  rate 
during  the  control  run  and  the  rejection  limits  were  changed 
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to  0.545  sec  and  1.5  sec  to  accommodate  this.  Plots  of  R-R 
intervals  were  visually  examined  and  the  output  was  passed 
to  EXCEL  where  means  and  standard  deviations  were 
calculated. 

To  analyze  the  spectral  components  of  the  signal,  the 
aperiodic  trends  are  removed  from  the  R-R  interval  plot. 
This  is  achieved  by  passing  the  data  points  through  a  21- 
point  moving  window  filter  (cubic  polynomial) ,  then 
subtracting  the  derived  aperiodic  trend  from  the  original 
signal.  The  spectral  analysis  was  performed  using  the  VSD 
GEDAP  routine.  Graphs  of  the  spectral  components  were 
visually  examined.  The  variance  in  each  of  the  two  bands 
under  consideration  (Baroreceptor  0.06-0.1  Hz;  Sinus 
Arrythmia  0.12-0.4  Hz)  was  calculated  by  the  program. 

R-R  interval,  heart  rate,  and  spectral  variance  were 
analyzed  over  the  full  270  second  data  collection  period 
during  the  following  segments  of  the  experiment: 

Rest 

Epoch  1-3 

Epoch  6 

Epoch  10-12 

Recovery  (where  available) 

Analysis  was  performed  on  3  of  the  4  subjects  for 
Experiment  L6  the  2  g  z-axis  exposure  at  32  shocks  per 
minute.  The  4th  subject's  ECG  data  could  not  be  analyzed 
due  to  baseline  shifts  and  movement  artifacts  within  the 
signal.  Data  from  all  4  subjects  were  analyzed  for 
Experiment  LI  (the  control)  and  Experiment  L5  (3  g  shocks 
1/2.5  minutes) . 

Given  the  small  number  of  subjects  and  the  stage  the 
analysis  has  reached  to  date,  no  statistical  analysis  has 
been  undertaken. 

-Electromyography 

Electromyographic  (EMG)  data  were  analyzed  to  determine 
if  paraspinal  muscle  fatigue  was  produced  during  one  hour 
exposures  to  x-axis  and  y-axis  impacts  and  during  two  hour 
exposures  to  z-axis  impacts.  EMG  signals  were  obtained 
during  30  second  static  back  extensions  of  approximately  50% 
maximal  effort,  performed  immediately  before  and  after  one 
and  two  hour  exposures.  Both  mean  frequency  (ME)  and 
integrated  EMG  (lEMG)  time  series  were  generated  from  the  30 
second  signals.  Signal  processing  to  produce  a  MF  time 
series  from  the  raw  EMG  data  involved  the  following  steps: 


subtraction  of  the  signal  mean  to  remove  DC  offset; 
extraction  of  thirty  consecutive  one  second  windows; 
calculation  of  the  variance  spectral  density  of  each  window; 
determination  of  the  zero  and  first  moments  of  each 
spectrum,  and  finally  calculation  of  MF  as  the  ratio  of  the 
first  to  the  zero  moment  for  each  window.  Each  window  was 
assigned  a  time  value  centered  on  the  window,  such  that  the 
first  window  corresponds  with  0.5  seconds,  the  second  with 
1.5  seconds,  and  so  on. 

The  mean  power  frequency  is  expected  to  decay 
exponentially  with  time  during  a  sustained  submaximal 
contraction  (Hagg  and  Suurkula,  1991,  Deeb,  Drury,  and 
Pendergast,  1992,  and  Kadefors  et  al . ,  1978).  The  initial 
mean  frequency  at  time  zero  (MPq)  and  the  time  constant  (T) 
of  the  decay  in  MF  with  time  were  determined  from  a  least 
squares  linear  regression  of  the  natural  logarithm  (In)  of 
MF  against  time. 

MF  =  MFq  e  -t/t 

Signal  processing  to  produce  the  lEMG  time  series  from 
raw  EMG  data  involved  the  following  steps :  subtraction  of 
the  signal  mean  to  remove  DC  offset;  full  wave  rectification 
of  the  signal;  and  integration  of  sequential  0.5  second 
windows.  A  least  squares  linear  regression  was  then  applied 
to  the  lEMG  time  series  to  determine  a  slope  (m)  and  an 
initial  lEMG  value  (Iq)  • 

lEMG  =  m  t  +  Iq 

The  pre-  and  post-exposure  mean  and  standard  deviation  of 
MFq,  X,  m,  and  Iq  were  compared  within  conditions  to 
determine  whether  there  was  evidence  of  paraspinal  muscle 
fatigue,  indicated  by  either  a  decrease  in  MFq  or  x,  and  an 
increase  in  either  m  or  Ig. 

The  EMG  burst  magnitude  in  response  to  impacts  of 
differing  frequency,  amplitude  and  direction  was  quantified 
from  the  short-term  exposures.  The  muscle  response  to 
individual  impacts  was  quantified  as  the  difference  in 
magnitude  of  lEMG  measured  within  a  125  msec  window  placed 
one  second  before  an  impact  {lEMG  rest) ,  compared  to  the 
lEMG  measured  within  a  window  of  the  same  size  occurring  at 
a  specified  time  delay  after  an  impact  (lEMG  burst) .  This 
analysis  was  performed  on  both  2  g  and  3  g  shocks  for  all 
three  axes.  A  window  duration  of  125  msecs  was  chosen  after 
visual  examination  of  typical  burst  responses  to  impact  in 
each  of  the  three  axes.  The  seat  acceleration  signal  for_ 
the  dominant  impact  axis  was  used  to  determine  the  specific 
time  of  impact  acceleration  peaks .  The  location  of 
acceleration  peaks  was  used  to  position  lEMG  windows  on  the 


EMG  signals.  The  lEMG  rest  window  was  placed  one  second 
before  the  time  of  the  impact,  and  the  lEMG  burst  window  was 
placed  at  six  different  positions  around  the  time  of  peak 
impact.  The  lEMG  burst  windows  were  located  125,  150,  and 
175  milliseconds  before  the  impact,  and  25,  50,  and  75 
milliseconds  after  the  impact.  Window  placement  prior  to 
the  peak  impact  acceleration  was  intended  to  identify  any 
muscle  response  to  the  motion  of  the  seat  that  preceded  the 
peak  acceleration. 

Processing  of  the  EMG  signal  involved  removing  the  mean 
to  eliminate  DC  offset,  rectification,  signal  averaging 
using  a  5  millisecond  moving  average,  and  finally 
integration  of  the  specified  rest  and  burst  windows .  The 
greatest  difference  between  lEMG  rest  and  the  series  of  lEMG 
burst  values  was  accepted  as  the  optimal  value  for  further 
analysis.  Means  and  standard  deviations  were  then  compared 
to  identify  trends  in  the  lEMG  response  across  impact 
frequency,  amplitude,  and  direction. 

A  qualitative  assessment  of  the  muscle  response  to 
impacts  of  differing  frequency,  amplitude  and  direction  was 
performed  on  selected  plots  of  EMG  and  acceleration  data. 
This  assessment  involved:  description  of  EMG  burst  patterns; 
the  relationship  of  the  EMG  response  to  displacement, 
velocity,  and  acceleration;  and  the  relationship  of  EMG  to 
the  internal  pressure  response. 

Biochemistry 

Biochemical  measures  were  obtained  from  the  four  subjects 
in  four  of  the  seven  experimental  conditions.  All  the  data 
were  time  dependent.  Some  dependent  variables  were  not 
measured  in  all  blood  samples,  either  because  of  difficulty 
in  obtaining  blood,  or  because  of  obvious  hemolysis  of  the 
blood  sample  which  affected  the  reliability  of  the  measure. 

The  data  are  reported  as  the  group  mean  and  standard 
deviation  of  the  mean.  Each  dependent  variable  was  analyzed 
by  a  one  way  repeated  measures  multiple  analysis  of  variance 
(MANOVA)  with  either  time  or  condition  as  the  within 
subjects  factor.  The  small  sample  size  prevented 
examination  of  the  interaction  of  time  and  condition  by  two 
way  repeated  measures  MANOVA.  Significant  differences  were 
accepted  at  an  alpha  value  of  less  than  0.05. 

Some  individual  data  were  plotted  for  variables  in  which 
a  trend  was  observed  in  the  group  mean  data.  These  plots 
were  examined  for  trends  which  may  have  been  obscured  by  the 
group  data. 

For  dependent  biochemical  variables  which  showed  no 
change  during  vibration  exposure,  but  did  show  a  change 
several  hours  post -exposure,  the  data  for  pre-,  mid-,  and 
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post-exposure  were  averaged  and  compared  with  data  averaged 
at  12,  24  and  36-hours  post-exposure.  These  data  were  also 
examined  by  a  one  way  repeated  measures  multiple  analysis  of 
variance  (MANOVA)  with  either  time  or  condition  as  the 
within  subjects  factor. 


Results 


Shock  directions  at  the  seat 

Examination  of  acceleration  profiles  indicated  that  the 
directions  of  accelerations  input  at  the  seat  were  +x,  +y 
and  -z.  These  were  the  opposite- -tO'  those  which  had 
originally  been  assiamed.  The  directions  of  seat 
accelerations  were  confirmed  independently  by  analysis  of 
the  Optotrak  displacement  data. 

The  uncertainty  regarding  the  direction  of  seat 
accelerations  originated  from  the  mounting  of  the  MARS 
accelerometers,  and  the  direction  of  table  displacement. 

The  accelerometers  providing  output  from  the  MARS  are 
mounted  on  the  underside  of  the  table  in  a  triaxial 
arrangement  and  do  not  conform  with  the  normal  biodynamic 
axes.  It  was  unclear  whether  a  correction  factor  was 
provided  in  the  MARS  software,  and  as  the  table  is  normally 
used  to  generate  a  random  vibration  signal,  it  had  not 
previously  been  of  importance. 

The  control  signal  generated  at  B.C.  Research  and 
provided  to  the  MARS  was  a  damped  sinusoid  having  a  positive 
peak  acceleration  followed  by  a  (lower)  negative  peak 
acceleration  (see  Phase  2  report) .  It  was  assumed  that  the 
output  response  of  the  table  should  be  a  positive 
displacement,  and  this  observation  was  used  to  confirm  the 
direction  of  acceleration.  It  was  subsequently  determined 
(from  analysis  of  the  Optotrak  displacement  and  acceleration 
data)  that  to  compensate  for  any  cumulative  displacement  of 
the  table  with  time,  the  MARS  software  generates  a  modified 
waveform.  To  maintain  a  neutral  reference  position  while 
responding  to  random  accelerations,  the  MARS  software 
applies  a  correction  factor  to  the  input  control  signal  to 
maintain  a  net  zero  displacement  in  the  time  domain.  Hence 
when  presented  with  a  shock  consisting  of  a  single  positive 
acceleration,  the  MARS  software  superimposes  a  low  amplitude 
negative  acceleration  to  balance  the  displacement  waveform. 
The  resultant  control  signal  delivered  to  the  table  appears 
as  a  transient  negative  displacement  followed  by  a  return  to 
a  neutral  position.  The  peak  (positive)  acceleration  occurs 
at  the  peak  (negative)  displacement,  rather  than  at  the 
leading  edge  of  a  (positive)  displacement  as  would  be 
expected.  Although  the  output  acceleration  waveform  at  the 
table  closely  approximates  the  input  acceleration  signal  to 


the  controller,  the  displacement  waveform  is  very  different 
from  that  expected  in  that  the  table  undergoes  a  negative 
displacement  in  response  to  a  positive  acceleration  signal. 
This  artifact  in  the  displacement  data  due  to  positional 
compensation  resulted  in  a  misinterpretation  of  the 
*»  direction  of  acceleration  during  the  programming  of  the  MARS 

control  signals. 

The  reasons  that  the  true  direction  of  acceleration  were 
^  not  identified  during  data  collection  are  threefold.  First, 

as  explained  above,  a  positive  transient  displacement 
observed  at  the  seat  corresponded  to  a  negative  peak 
acceleration;  second,  the  accelerometers  were  calibrated 
using  the  RMS  value  of  a  swept  sine  wave,  and  hence  the 
calibration  procedure  was  insensitive  to  direction;  and 
third,  the  test  signal  displayed  prior  to  each  data 
collection  consisted  of  random  ±  Ig  shocks . 

As  the  long  term  experiments  contained  both  positive  and 
negative  shocks  these  experiments  were  unaffected  by  this 
error.  The  resultant  shocks  delivered  during  the  short  term 
exposures  were  in  the  +x,  +y  and  -z  directions.  As  the  body 
is  symmetrical  about  the  saggital  (x-z)  plane,  the  response 
to  ±  g  shocks  in  the  y-axis  should  be  identical  but  opposite 
in  sign.  The  response  to  shocks  in  the  x-and  z-axes, 
however,  will  differ  according  to  the  direction  of  the 
shocks .  In  the  pilot  study,  due  to  time  constraints  and 
experimental  design,  it  was  only  possible  to  collect  shock 
data  in  one  direction  for  each  axis.  It  is  intended  that 
both  positive  and  negative  shocks  will  be  analyzed  in  these 
axes  in  the  full  experimental  protocol  (phase  4) . 

Analysis  of  individual  shocks 

The  input  control  signal  developed  at  B.C.  Research 
consisted  of  a  background  random  vibration  having  a 
frequency  spectriom  of  2  to  40  Hz.  This  spectrum  was  chosen 
to  correspond  with  the  band-pass  filter  routine  contained  in 
the  MARS  software,  and  automatically  applied  to  control 
signals.  The  individual  shocks  were  inserted  into  the 
background  vibration  record,  but  to  avoid  variance  in  the 
absolute  value  of  peak  acceleration,  the  background 
vibration  was  not  superimposed  on  the  shock  profile. 

Inspection  of  the  output  data  indicated  that  the  results 
obtained  from  the  peak  detection  program  were  unreliable  for 
0.5  g  shocks.  This  was  due  to  difficulties  in  the  software 
correctly  identifying  the  shock  peaks  within  the  background 
vibration  signal.  Although  the  background  vibration  level 
of  0.5  m.s  rms  represented  only  one  tenth  of  the  peak 
acceleration  input  of  the  shock  waveform,  the  signal  was  ISO 
weighted,  and  contained  crest  factors  of  3  to  5.  Hence,  as 
the  unweighted  signal  was  approximately  1.0  m.s  rms,  the 
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peak  amplitude  of  some  frequency  components  approached  that 

of  the  0.5  g  shock  waveform.  Although  it  was  possible  to 

identify  the  shocks  visually  within  the  acceleration  record, 

the  peak  detection  program  was  not  sufficiently 

sophisticated  to  provide  consistent  identification  of 

shocks.  A  similar  problem  was  experienced  in  identifying 

the  shock  response  of  spinal  acceleration,  internal  pressure 

and  Respitrace.  This  was  frequently  due  to  the  absence  of  a 

significant  response,  particularly  in  the  internal  pressure 

and  Respitrace  records.  For  these  reasons,  the  results 

obtained  for  the  0.5  g  exposures  are  incomplete  and  have  not  > 

been  reported. 

Analysis  of  results  showed  that  the  shock  acceleration 
waveforms  measured  by  the  seatpad  accelerometer  contained 
high  frequency  resonances  superimposed  on  the  sinusoidal 
shocks  generated  as  input  to  the  MARS.  These  frequency 
components  were  most  likely  attributable  to  the  resonances 
of  the  seat-subject  subsystem  mounted  on  the  shake  table. 

These  resonances  were  greatest  in  the  4  Hz  shocks  at  3  g  in 
the  X  and  y  directions,  and  decreased  in  severity  with 
increasing  shock  frequency  and  decreasing  shock  amplitude. 

A  sample  of  the  unfiltered  acceleration  data  measured  at  the 
seat  (Sx) ,  and  spine  (L2  x  and  T1  x)  in  response  to  a  3  g,  4 
Hz  shock  in  the  x-axis  is  shown  in  Figure  1 .  It  can  be  seen 
that  the  higher  frequencies  contained  in  the  seat 
acceleration  data  tend  to  magnify  the  peak  of  the 
fundamental  shock  waveform.  These  high  frequency  resonances 
at  the  seat  do  not  transmit  to  the  accelerometers  attached 
to  the  spine.  As  the  objective  of  the  experiments  was  to 
evaluate  the  effect  of  the  nominal  shock  frequency  on  seat 
to  spine  transmission,  the  data  was  filtered  at  various  cut¬ 
off  frequencies  to  remove  the  unwanted  higher  frequencies 
without  attenuating  the  underlying  shock  waveform.  It  was 
found  that  a  filter  cut-off  frequency  of  60  Hz  removed  many 
of  the  higher  frequencies  from  the  4  Hz  shock  at  the  seat . 

Tests  using  a  lower  cut-off  frequency  of  40  Hz  tended  to 
attenuate  the  3  g,  11  Hz  shocks  (>  10%  attenuation) .  Hence, 
all  acceleration  data  were  filtered  at  60  Hz  prior  to 
analysis.  The  effects  of  the  60  Hz  low  pass  filter  on  the  4 
Hz  and  11  Hz  3g,  x-axis  shocks  are  shown  in  Figures  2  and  3. 

The  low-pass  filter  cut-off  did  not  affect  the  spinal 
acceleration  data  in  response  to  the  x-axis  or  y-axis 
shocks.  Examples  of  unfiltered  and  filtered  accelerations 
measured  at  the  lumbar  spine  (L2  x)  are  shown  in  Figure  4. 

The  high  frequency  resonances  superimposed  on  the  z-axis  ^ 

shocks  were  less  severe  than  in  the  x-axis  and  y-axis 
shocks.  However,  high  frequency  accelerations  occurred  at 
the  seat  following  the  4  Hz  shocks  at  -3  g  (see  Figure  5) . 

This  was  attributed  to  the  subject  leaving  the  seat  during 
the  downwards  acceleration  of  the  table,  and  then  striking 
the  seat  as  the  table  returned  to  the  neutral  position.  The 
frequency  of  these  oscillations  were  generally  >  100  Hz,  and 
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were  removed  from  the  raw  data  by  the  60  Hz  low  pass  filter 
as  shown  in  Figure  5 .  This  high  frequency  resonance  at  the 
seat  was  only  evident  in  the  lower  frequency  (4,  5  and  6 
Hz)  and  higher  amplitude  (-3  g  and  -2  g)  shocks,  and 
diminished  with  increasing  frequency  and  decreasing 
amplitude.  However,  unlike  the  x-axis  and  y-axis  shocks, 
results  show  that  these  higher  frequencies  were  transmitted 
from  the  seatpad  to  the  spinal  accelerometers  in  some 
subjects.  An  example  of  this  effect  in  the  acceleration 
measured  at  the  lumbar  spine  {L4  z)  is  shown  in  Figure  6. 

The  results  also  show  that  the  peak  acceleration  at  the 
spine  in  response  to  shocks  in  the  -z  direction  occurred 
when  the  subject  hit  the  seat  and  not  as  a  direct  response 
to  the  seat  acceleration  waveform.  This  effect  is  more 
pronounced  in  the  higher  amplitude  and  lower  frequency 
shocks .  As  the  maximum  acceleration  values  were  used  in 
calculating  the  spine/seat  acceleration  ratio,  the  results 
are  expressed  as  a  shock  response  ratio  rather  than  as  a 
transmission  ratio.  The  underlying  period  of  the  shock 
waveform  at  the  spine  was  also  much  shorter  than  that  of  the 
shock  input  at  the  seat,  indicating  a  higher  frequency  event 
when  the  subject  hit  the  seat.  In  some  shock  conditions, 
the  peak  acceleration  response  of  the  spine  included 
frequency  components  in  the  range  of  30  to  80  Hz.  This 
occurred  when  responding  to  the  higher  magnitude,  lower 
frequency  shocks.  In  a  few  cases  this  resulted  in 
attenuation  of  acceleration  peaks  due  to  low  pass  filtering. 

The  internal  pressure  data  showed  a  clear  response  to 
shocks  in  the  x,  y  and  z-axes.  The  waveform  of  internal 
pressure  showed  no  evidence  of  the  high  frequency  components 
(>100  Hz)  seen  in  the  seatpad  acceleration  records.  The 
most  severe  response  was  to  negative  z-axis  shocks. 

Internal  pressure  showed  a  very  similar  pattern  to  that  of 
spinal  acceleration,  with  a  rapid  increase  in  pressure  as 
the  subject  hit  the  seat  after  the  shock.  Although  this 
waveform  also  contained  high  frequency  components  {>20  Hz) , 
no  evidence  was  found  of  peak  attenuation  due  to  low  pass 
filtering  at  60  Hz  with  the  exception  of  -3  g  z-axis  shocks, 
where  minor  attenuation  occurred  in  a  few  subjects  at  the 
lowest  frequency  (see  Figure  7) , 

Vertebra-skin  transfer  function 

The  natural  frequencies  (GJn)  damping  ratios  (^) 

calculated  for  the  vertebra-skin  transfer  functions  were 
generally  in  the  range  200  <  On  C  0.4, 

respectively.  No  consistent  differences  were  found  between 
the  data  in  the  y  and  z  directions,  or  at  the  lumbar  and 
thoracic  levels .  The  natural  frequency  and  damping  ratio  of 
the  accelerometer-tissue  sub-system  of  each  subject  at  the 
lumbar  and  thoracic  levels  are  presented  in  Table  1.  The 


43 


frequencies  of  accelerations  measured  at  the  spine  in 
response  to  the  x-axis  and  y-axis  shocks  were  lower  than 
those  calculated  for  the  vertebra-skin  subsystem  (03^)  . 

Hence  the  application  of  the  inverse  transfer  function  had 
very  little  effect  on  the  spinal  acceleration  data.  In  some 
instances,  the  unfiltered  spinal  accelerations  measured  in 
response  to  the  z-axis  response  contained  frequency 

components  greater  than  G3n-  Application  of  the  inverse 
transfer  function  caused  these  components  to  be  amplified 
and  severely  distorted  the  transformed  acceleration  record. 
This  effect  was  considered  to  be  an  artifact  due  to  the 
asstimption  that  the  tissue  behaves  as  a  simple  Kelvin 
element.  To  remove  this  effect  the  acceleration  data  were 
low  pass  filtered  at  60  Hz  prior  to  the  application  of  the 
inverse  transfer  function.  This  treatment  of  the  data  gave 
a  satisfactory  result  in  most  cases.  An  example  of  spinal 
acceleration  measured  at  the  l\jmbar  level  (L4  z)  unfiltered, 
filtered  at  60  Hz,  and  subsequently  corrected  by  application 
of  the  inverse  transfer  function  is  shown  in  Figure  8. 

Spinal  acceleration 

Spine  (T1  and  L2)  x-axis  acceleration  response  to  positive 
x-axis  shocks  at  the  seat 

Figure  9  (Appendix  D)  shows  a  typical  plot  of  a  4  Hz,  +3 
g  shock  in  the  x-axis  at  the  seat.  Superimposed  is  the 
response  from  the  Ixambar  vertebrae  (L2),  followed  in  time  by 
the  thoracic  vertebrae  (Tl) .  Figure  10  shows  the  same 
responses  for  a  +3  g,  11  Hz  shock  in  the  x-axis.  A 
consistent  feature  of  the  spinal  response  to  shocks  in  the 
x-axis  was  that  the  acceleration  output  at  the  thoracic 
level  (Tl)  was  the  reverse  of  the  acceleration  input  at  the 
seat.  Hence,  transmission  ratios  Tl  x/Sx  are  negative 
quantities.  This  result  suggests  that  the  sudden  forwards 
acceleration  of  the  seat  induces  a  backwards  rotation  of  the 
upper  body  in  the  sagittal  (x  -  z)  plane.  Anatomically, 
this  motion  can  be  achieved  by  extension  of  the  hip  joint 
and/or  extension  of  the  spine. 

The  average  delay  between  the  seat  impact  and  lumbar 
response  for  the  3  g  shocks  ranged  from  0.04  sec.  to  0.06 
sec.  The  average  delay  between  seat  impact  and  thoracic 
response  ranged  from  0.06  to  0.07  sec.  Delays  were  slightly 
longer  for  Tl  response  compared  with  L2  response.  Although 
the  actual  differences  between  frequencies  were  very  small, 
the  4  Hz  shock  produced  the  shortest  delay  at  L2,  but  the 
longest  delay  at  Tl .  Delays  for  2  g  shocks  were  0.04  to 
0.06  for  L2  and  0.05  to  0.08  for  Tl.  For  1  g  shocks,  delays 
were  0.05  to  0.09  for  L2  and  0.07  to  .10  for  Tl .  Delays  for 
all  x-axis  shocks  across  the  10  subjects  are  presented  with 
transmission  data  in  Tables  2  to  7. 
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The  transmission  ratio,  a  non-dimensional  unit,  was 
calculated  as  the  peak  of  the  lumbar  or  thoracic 
acceleration  response  (m.s  -2),  divided  by  the  peak 
acceleration  input  measured  at  the  seat  (m.s  "2).  X-axis 
transmission  ratios  for  all  10  subjects  for  3  g  shocks  from 

♦  seat  to  lumbar  vertebrae  are  shown  in  Table  2.  Each  subject 
responded  to  2  shocks  at  the  shock  frequencies  of  4,  5,  6,  8 
and  11  Hz  presented  in  a  random  order  and  timing.  Table  2 
also  provides  the  average  transmission  ratio  across  all  10 

*  subjects  for  each  shock  frequency  and  the  average  delays  in 
a  summary  table.  Responses  of  each  individual  subject  are 
plotted  as  line  graphs  in  Figure  11a.  Superimposed  on  the 
10  lines  are  data  points  indicating  the  mean  response  at 
each  frequency  for  the  10  subjects  (20  shocks) 

Average  seat  to  Ixambar  (L2)  transmission  ratios  ranged 
from  0.2  to  0.49  in  magnitude,  although  individual  subject's 
responses  ranged  up  to  0.725.  Figure  11a  shows  that 
although  the  range  of  response  between  subjects  is  large, 
the  trends  in  response  between  subjects  are  quite  similar. 
Most  subjects  (8  of  10)  demonstrated  a  peak  response  at  5 
Hz,  Three  subjects  had  a  second  peak  at  8  Hz,  and  for  all 
subjects  the  response  was  lowest  at  the  11  Hz  shock. 

X-axis  seat  to  thoracic  (Tl)  transmission  data  are  shown 
for  10  subjects  in  Table  3  and  plotted  in  Figure  12a. 
Transmission  ratios  were  lower  (0,13  to  0.38)  than  at  the 
lumbar  spine.  The  pattern  of  responses  is  similar  to  the 
Iximbar  response.  Eight  of  10  subjects  demonstrated  a  peak 
response  at  5  Hz  and  3  showed  a  tendency  toward  a  second 
peak  response  at  8  Hz.  All  subjects,  once  again,  had  a 
lower  response  to  11  Hz  shocks. 

X-axis  lumbar  (L2)  response  (L2  x/Sx)  for  2  g  shocks  is 
shown  in  Table  4  and  Figure  lib.  At  2  g  shocks,  the 
tendency  toward  a  5  Hz  peak  is  less  evident,  with  5  of  10 
subjects  responding  with  higher  transmission  ratios  at  5  Hz 
compared  with  4  Hz .  The  average  response  of  the  group 
differed  very  little  between  4,  5  and  6  Hz.  The  second  peak 
at  8  Hz  found  with  the  3  g  shocks  is  not  apparent  with  2  g 
>  shocks.  All  10  subjects  displayed  a  reduction  in  response 

from  6  Hz  to  11  Hz.  The  magnitude  of  the  transmission  ratio 
varied  from  0.19  to  0.39.  This  was  consistently  lower  at 
any  particular  frequency  than  the  transmission  at  3  g. 

*  X-axis  thoracic  (Tl)  acceleration  transmission  for  2  g 

shocks  (Tl  x/Sx)  is  shown  in  Table  5  and  Figure  12b. 

Average  levels  of  transmission  range  from  0.11  to  0.32, 
slightly  lower  than  levels  measured  at  3  g.  There  were 
equal  or  higher  levels  of  Tl  transmission  at  4  Hz  compared 
with  5  Hz  for  8  of  the  10  subjects  at  2  g,  suggesting  the 
dominant  frequency  for  2  g  shocks  is  4  Hz  or  lower. 

However,  at  l\ambar  transmission,  the  4  Hz  response  was 
dominant  for  fewer  subjects  (6  of  the  10) ,  resulting  in  a 
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flatter  mean  response.  The  dominant  peak  at  5  Hz  in  the  3  g 
transmissions  is  not  evident  for  many  subjects  at  2  g.  The 
overall  response  appeared  much  flatter,  possibly  with  a 
dominant  peak  at  4  Hz  rather  than  5  Hz. 

X-axis  lumbar  {L2)  responses  for  1  g  shocks  (L2  x/Sx)  are 
shown  in  Table  6  and  Figure  11c.  Each  subject  responded  to 
two  shocks  at  frequencies  of  2,  4,  6,  8,  and  11  Hz.  The 
response  to  1  g,  2  Hz  shocks  was  higher  than  an  any  other 
frequency  and  shock  amplitude  (0.52).  Transmission  of  4  Hz 
shocks  was  similar  at  1  g  compared  with  2  g,  but  the 
remaining  responses  at  1  g  were  slightly  less  than  that  at  2 
g.  Responses  at  1  g  ranged  from  0.17  to  0.52,  increasing  in 
magnitude  with  decreasing  shock  frequency.  It  appears  that 
the  lower  frequency  shocks,  associated  with’  the  greatest 
displacement,  also  result  in  the  greatest  transmission  of 
acceleration.  With  the  exception  of  two  subjects.  Figure 
11c  shows  a  similar  trend  of  individual  responses  with  each 
subject's  response  decreasing  as  frequency  increased  from  2 
Hz  to  11  Hz.  Since  there  is  no  data  lower  than  2  Hz,  it  is 
impossible  to  determine  whether  the  dominant  frequency  is  2 
Hz,  or  at  an  even  lower  frequency. 

X-axis  thoracic  (T1  x/Sx)  response  at  1  g  is  similar  to 
lumbar  response  (Table  7,  Figure  12c) .  Transmission  ratios 
ranged  from  0.10  at  11  Hz  up  to  0.50  at  2  Hz.  Again,  the  2 
Hz  response  was  dominant  and  the  4  Hz  response  similar  in 
magnitude  to  2  g.  At  the  remaining  frequencies,  the 
responses  were  slightly  lower  than  responses  at  2  g,  and 
lower  also  compared  with  lumbar  transmissions.  There  was  a 
trend  for  all  individuals  of  decreasing  transmission  from  2 
Hz  to  11  Hz.  However,  a  few  subjects  showed  a  slight 
flattening  or  increased  response  between  6  and  8  Hz. 

In  general,  the  dominant  shock  frequencies  at  3  g  are  not 
the  same  for  lower  (2  g  and  1  g)  amplitudes.  The  pattern  of 
response  changes  from  a  5  Hz  peak  at  3  g  to  a  4  Hz  at  2  g 
and  2  Hz  at  1  g.  For  1  and  2  g  shocks,  the  dominant 
frequency  may  be  lower  still,  but  no  data  could  be  collected 
to  test  responses  at  ranges  lower  than  those  presented. 

The  acceleration  response  of  the  body  at  the  spine  to  an  * 

input  at  the  seat  appears  to  be  non-linear  for  various 
amplitudes  of  shocks.  The  highest  average  transmission 
ratios  were  found  for  3  g  shocks  and  lowest  for  1  g  shocks. 

In  addition,  a  second  resonance  frequency  began  to  appear  in  • 

the  response  to  3  g  shocks,  which  was  not  evident  in  the 
response  to  1  and  2  g  shocks.  It  cannot  be  determined  from 
this  study  whether  these  trends  would  continue  at  higher  and 
lower  amplitudes. 
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Spine  (T3  and  L4)  z-axis  acceleration  response  to  positive 
x-axis  shocks  at  the  seat 

Figure  13  shows  the  z-axis  response  at  the  thoracic  {T3) 
and  lumbar  (L4)  vertebrae  to  a  positive  3  g,  4  Hz,  x-axis 
shock  at  the  seat  (T3  z/Sx) .  A  positive  x-axis  shock 
produced  a  negative  z-axis  acceleration  at  both  levels  of 
the  spine.  The  average  delays  between  peak  seat 
acceleration  response  and  peak  L4  acceleration  in  the  z-axis 
are  from  0.036  to  0.047  sec  and  at  T3  are  from  0.06  to  0.076 
sec.  These  delays  are  similar  to  those  reported  for  seat  to 
lumbar  and  thoracic  transmission  in  the  x-axis. 

Transmission  of  x-axis  3  g  shocks  from  the  seat  to  the  z- 
axis  of  the  Itimbar  spine  (L4  z/Sx)  for  all  10  subjects  are 
shown  in  Table  8  and  Figure  14a.  The  average  transmission 
of  an  x-axis  input  is  higher  in  the  lumbar  z-axis  (0.29  to 
0.67)  than  in  the  lumbar  x-axis  (0.22  to  0.49)  for  all 
frequencies  by  18  to  27%.  A  5  Hz  dominant  trend,  as  found 
in  the  x-axis  spinal  response,  is  seen  in  the  z-axis  spinal 
response  in  eight  of  the  ten  subjects.  All  but  one  of  the 
ten  subjects  show  a  lower  response  at  6  Hz.  Six  subjects 
responded  to  8  Hz  with  a  higher  transmission  than  6  Hz, 
indicating  a  second  peak  in  transmission.  However,  the 
differences  in  response  between  the  6  and  8  Hz  frequencies 
was  usually  very  small.  All  responses  decrease  towards  11 
Hz. 

Z-axis  thoracic  (T3)  responses  (0.23  to  0.50)  were,  on 
average,  24  to  43%  higher  than  the  x-axis  thoracic  responses 
(0.13  to  0.38)  to  the  same  x-axis  shock  inputs  (Table  9  and 
Figure  15a) .  The  T3  response  pattern  in  the  z-axis  was 
similar  to  the  z-axis  response  at  L4  with  a  dominant 
transmission  at  5  Hz  for  8  subjects.  Six  of  the  10  subject 
had  a  second  slight  peak  at  8  Hz. 

The  2  g,  z-axis  Itimbar  and  thoracic  responses  to  x-axis 
input  (Tables  10  and  11,  Figures  14b  and  15b)  were  similar 
to  one  another.  There  was  a  flat  response  between  4,  5,  6 
and  8  Hz,  and  a  reduced  transmission  ratio  at  11  Hz.  The 
magnitude  of  response  for  both  thoracic  and  lumbar  vertebrae 
at  2  g  (0.20  to  0,50)  was  slightly  less  than  at  3  g  (0.23  to 
0.67)  . 

Z-axis  lumbar  and  thoracic  responses  to  x-axis  1  g  inputs 
are  shown  in  Tables  12  and  13,  and  Figures  14c  and  15c. 
Again,  L4  z-axis  responses  (0.23  to  0.61)  are  higher  than  T3 
z-axis  (0,19  to  0.57)  and  show  a  fairly  consistent  trend 
from  a  peak  response  at  2  Hz  downward  toward  11  Hz.  The 
transmission  ratios  at  lumbar  and  thoracic  vertebrae  for  1  g 
shocks  were  usually  lower  than  for  2  g  and  3  g  shocks.  As 
with  other  amplitudes,  the  z-axis  thoracic  response  was 
higher  than  the  x-axis  thoracic  response  by  13  to  48%  and 


the  z-axis  lumbar  response  was  higher  than  the  x-axis  by  14 
to  28%,  depending  upon  frequency. 

Spine  (T2  and  L3 )  y-axis  acceleration  response  to  positive 
y-axis  shocks  at  the  seat 

Figure  16  illustrates  a  y-axis  shock  (3  g,  4  Hz)  at  the 
seat.  Superimposed  on  the  seat  shock  is  the  liimbar  (L3) 
response  and  the  thoracic  (T2)  response  to  the  shock. 

Figure  17  shows  the  seat,  Ixambar  and  thoracic  accelerations 
for  a  3  g,  11  Hz  shock  in  the  y-axis.  The  y-axis  shocks  and 
responses  were  similar  in  pattern  to  the  x-axis  data  shown 
in  Figures  9  and  10.  The  initial  acceleration  response  at 
the  thoracic  level  (T2)  was  in  the  opposite  direction  from 
both  the  lumbar  (L3)  response,  and  from  the  shock  input  at 
the  seat.  The  response  at  T2  tended  to  be  biphasic, 
indicating  a  whipping  action  in  response  to  a  y-axis  shock. 
In  some  cases,  the  positive  acceleration  peak  exceeded  the 
initial  negative  peak.  However,  to  remain  consistent,  the 
transmission  ratios  were  in  all  cases  reported  as  the 
negative  (peak)  acceleration  response  to  seat  input 
acceleration . 

Transmission  of  3  g  shocks  from  seat  to  lumbar  spine  in 
the  y-axis  (L3  y/Sy)  are  shown  in  Table  14  and  Figure  18a 
and  seat  to  thoracic  spine  in  the  y-axis  (T2  y/Sy)  are  shown 
in  Table  15  and  Figure  19a.  Delays  between  the  seat  input 
and  the  lumbar  response  averaged  from  0.06  to  0.07  seconds, 
depending  upon  the  shock  frequency.  The  delay  between  seat 
input  and  thoracic  response  was  similarly  between  0.06  and 
0.07  seconds.  Levels  of  transmission  from  seat  to  lumbar 
spine  were  higher  in  the  y-axis  than  that  found  in  the  x- 
axis,  but  levels  of  x-and  y-axis  transmission  from  seat  to 
thoracic  spine  were  similar.  Y-axis  transmission  of  3  g 
shocks  from  seat  to  lumbar  spine  ranged  from  0.23  to  0.64 
(compared  with  0.20  to  0.49  for  the  x-axis),  and  seat  to 
thoracic  spine  ranged  from  0.1  to  0.32. 

The  5  Hz  dominance  in  the  3  g  x-axis  data  was  not  found 
in  the  3  g  y-axis  data.  As  seen  in  Figures  18a  and  19a,  the 
dominant  frequency  response  for  both  lumbar  and  thoracic 
vertebrae  was  at  4  Hz,  with  lower  levels  of  transmission 
occurring  for  subsequently  higher  shock  frequencies  for  all 
but  one  subject. 

Y-axis  transmission  of  2  g  shocks  from  seat  to  lumbar  and 
seat  to  thoracic  vertebrae  are  shown  in  Tables  16  and  17, 
and  Figures  18b  and  19b.  Levels  of  transmission  of  2  g 
shock  were  very  similar  to  3  g  shocks  (0.23  to  0.64  for  seat 
to  Ixombar  and  0.1  to  0.32  for  seat  to  thoracic)  .  The 
pattern  of  y-axis  response  at  2  g  was  similar  to  the  3  g 
response,  with  the  highest  spinal  response  occurring  at  4 
Hz,  and  decreasing  at  higher  shock  frequencies. 
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Y-axis  transmission  of  1  g  shocks  from  seat  to  lumbar  and 
seat  to  thoracic  vertebrae  are  shown  in  Tables  18  and  19, 
and  Figures  18c  and  19c.  Transmission  from  seat  to  lumbar 
vertebrae  was  similar  at  3  g  and  2  g  for  4  and  11  Hz  shocks. 
The  highest  transmission  ratio  (0.7)  was  found  for  the  2  Hz 
shock  at  1  g.  Seat  to  thoracic  transmission  at  the  higher 
frequencies  was  slightly  lower  than  the  2  g  and  3  g  shocks . 
The  pattern  of  response  was  similar  to  the  2  g  and  3  g 
response  with  the  highest  transmission  occurring  for  the 
lowest  shock  frequency  (2  Hz)  . 

Spine  (T3  and  L4)  z-axis  acceleration  response  to  positive 
y-axis  shocks  at  the  seat 

Transmission  from  y-axis  shocks  at  the  seat  to  the  z-axis 
at  the  lumbar  and  thoracic  spine  for  3  g,  2  g  and  1  g  shocks 
{L4  z/Sy  and  T3  z/Sy)  are  shown  in  Tables  20  to  25  and 
Figures  20  to  21.  The  z-axis  spinal  response  to  y-axis 
shocks  was  very  small,  usually  less  than  0.2.  The  response 
across  frequencies  resembled  the  y-axis  response  at  the 
spine.  Often  the  peak  transmission  was  at  4  Hz,  and  the 
remaining  responses  are  somewhat  lower.  As  the  magnitudes 
were  low  relative  to  the  variance  between  individuals,  the 
mean  responses  often  reseinbled  a  flat  line.  No  clear 
conclusions  can  be  drawn  about  the  dominant  frequency  in  the 
z-axis  response  to  y-axis  shocks. 

Spinal  (T3  and  L4)  z-axis  acceleration  response  to  negative 
z-axis  shocks  at  the  seat 

The  nature  of  the  spinal  response  to  negative  z-axis 
shocks  differed  from  the  response  to  x-and  y-axis  shocks. 

The  major  acceleration  at  the  spine  took  place  not  as  a 
direct  response  to  the  downwards  acceleration  of  the  seat, 
but  rather  when  the  subject  contacted  the  seat  following  its 
return  towards  a  neutral  position.  As  the  magnitude  of  this 
positive  acceleration  peak  was  substantially  greater  than 
the  initial  (negative)  acceleration  at  the  spine,  it  was 
considered  appropriate  to  use  this  value  in  calculating  a 
response  ratio.  Responses  are,  therefore,  not  direct 
transmission  ratios  from  seat  to  spine,  but  response  ratios 
to  the  seat  acceleration  input.  Figure  22  shows  the  z-axis 
seat  acceleration  with  the  Ixambar  {L4)  and  thoracic  (T3) 
response  superimposed  for  a  negative  3  g  shock  at  4  Hz 
frequency.  An  11  Hz,  3  g  shock  and  spinal  responses  are 
shown  in  Figure  23. 

The  most  significant  features  of  the  response  to  shocks 
in  the  negative  z-axis  were  that  they  are  of  a  much  higher 
magnitude  than  the  responses  to  x-and  y-axis  shocks,  and 
that  the  shock  waveform  had  a  much  higher  frequency.  The 
most  severe  accelerations  occurred  in  response  to  the  lower 
frequency  shocks  at  the  seat.  The  spinal  acceleration 
waveform  in  response  to  these  shocks  generally  showed  one  or 


more  oscillations  of  rapidly  decaying  magnitude,  following 
the  initial  peak.  Analysis  of  several  wave  forms  suggested 
a  frequency  of  30  to  80  Hz  in  the  initial  peak  response, 
with  the  subsequent  lower  amplitude  oscillations  having 
progressively  lower  frequencies.  However,  as  responses 
differed  among  subjects  and  across  shock  frequencies,  it  is 
not  possible  to  make  a  general  statement  regarding  the 
waveform  of  spinal  acceleration. 

The  average  delays  in  spinal  response,  both  lumbar  and 
thoracic,  for  3  g  shocks  were  0.1  to  0.17  seconds.  For  the 

2  g  shocks,  lumbar  and  thoracic  delays  were  between  0.09  and 
0.16  seconds.  Delays  for  spinal  responses  to  1  g  shocks  are 
from  0.07  to  0.17  seconds.  Delays  were  slightly  longer  for 

3  g  shocks,  followed  by  2  g,  then  1  g  shocks. 

Tables  26  and  27  and  Figures  24a  and  25a  show  spinal 
responses  to  3  g  shocks  in  the  z-axis  (L4  z/Sz  and  T3  z/Sz) . 
The  trend  in  responses  at  L4  across  frequencies  was 
consistent  with  other  conditions.  The  highest  response 
occurred  for  4  Hz  shocks  and  the  lowest  for  11  Hz  shocks. 
Lumbar  responses  averaged  0.55  at  11  Hz,  and  increased  to 
2.3  at  4  Hz.  Responses  at  T3  were  similar,  however  3 
subjects  had  a  slightly  higher  response  to  the  5  Hz  shock 
than  the  4  Hz.  Average  magnitudes  of  response  at  T3  were 
0.56  at  11  Hz  and  1.96  at  4  Hz. 

Tables  28  and  29  and  Figures  24b  and  25b  show  spinal 
responses  to  2  g  shocks  in  the  z-axis.  Although  the  trend 
in  average  2  g  responses  resembled  those  at  3  g,  there  is 
considerably  more  variation  among  subjects.  Three  subjects 
had  higher  l\imbar  responses  to  5  Hz  shocks  than  4  Hz  and  one 
subject  had  a  peak  at  6  Hz  and  an  increase  in  response  from 
8  Hz  to  11  Hz.  Response  magnitudes  (0.66  to  2.3)  resembled 
those  at  3  g,  although  they  are  slightly  lower  at  most 
frequencies.  In  the  thoracic  response,  2  subjects  had 
higher  responses  at  5  Hz  than  4  Hz,  while  2  other  subjects 
had  maximal  responses  to  6  Hz  shocks.  Average  responses  at 
T3  ranged  from  0.5  at  11  Hz  to  2.5  at  4  Hz.  Responses  to  4, 
5  and  6  Hz  shocks  at  2  g  were  higher  than  for  the  same 
frequencies  at  3  g. 

Responses  to  1  g  shocks  in  the  z-axis  are  shown  in  Tables 
30  and  32  and  Figures  24c  and  25c.  The  L4  response  was 
lower  in  magnitude  than  2  or  3  g  responses  (0.6  to  1.15)  and 
fairly  flat  across  2,  4  and  6  Hz.  The  mean  peak  response 
was  at  4  Hz  and  7  subjects  had  a  greater  response  to  4  Hz 
than  to  2  Hz  shocks .  The  same  trend  was  found  in  the  T3 
responses,  although  the  peak  at  4  Hz  is  much  more  apparent.  . 
Nine  of  the  10  subjects  showed  a  peak  response  to  4  Hz 
shocks  The  differences  between  2,  4  and  6  Hz  shocks  were 
larger  than  at  L4.  Average  responses  at  T3  ranged  from  0.37 
at  11  Hz  to  1.2  at  4  Hz. 
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Spinal  (T1  and  L2)  x-axis  acceleration  to  negative  z-axis 
shocks  at  the  seat 

X-axis  responses  to  the  seat  z-axis  shocks  illustrate  a 
whiplash-like  effect  occurring  at  the  spine  (Figure  26) . 
Accelerations  at  both  lumbar  and  thoracic  levels  showed  both 

♦  positive  and  negative  phases  of  acceleration  in  most 
subjects.  The  thoracic  spine  showed  an  initial  positive  x 
acceleration  in  response  to  the  z-axis  shock,  suggesting 
flexion  of  the  spine.  This  was  followed  by  a  high  frequency 

^  biphasic  acceleration  "spike"  as  the  subject  struck  the 

seat,  and  then  a  more  prolonged  negative  acceleration, 
suggesting  a  correction  of  the  initial  spinal  flexion.  The 
lumbar  spine  showed  a  response  mainly  to  the  subject 
striking  the  seat  as  it  returned  to  a  neutral  position. 

This  resulted  in  a  negative  x-axis  acceleration  of  the 
Ivimbar  spine,  followed  by  a  positive  acceleration.  This 
motion  again  suggests  a  flexion,  followed  by  an  extension  of 
the  spine.  Responses  at  L2  and  Tl  to  3  g  shocks  (L2  x/Sz 
and  Tl  x/Sz)  are  shown  in  Tables  32  and  33  and  Figures  27a 
and  28a.  Liimbar  responses  ranged  from  0.5  at  11  Hz  to  2.5 
at  4  Hz.  The  trend  for  most  subjects  was  a  decreasing 
response  from  4  to  8  Hz  and  a  flattening  between  8  and. 11 
Hz.  Responses  at  Tl  (0.26  to  1.2)  were  half  the  magnitude 
of  L2 .  Although  the  average  response  pattern  was  similar  to 
that  at  L2,  there  was  considerably  more  variation  between 
subjects  with  one  showing  a  peak  at  5  Hz  and  two  showing  a 
peak  at  6  Hz . 

X-axis  responses  to  2  g  shocks  are  shown  in  Tables  34  and 
35  and  Figures  27b  and  28b.  The  overall  trend  in  responses 
to  2  g  shocks  was  similar  to  3  g,  with  dominant  responses 
generally  occurring  at  4  Hz.  There  were  differences, 
however,  in  both  magnitude  and  trend  of  response  between 
individual  subjects.  Average  responses  at  L2  ranged  from 
0.35  to  1.29  and  at  Tl  from  0.29  to  1.07. 

The  1  g  x-axis  responses  to  z-axis  shocks  are  shown  in 
Tables  36  and  37  and  Figures  27c  and  28c.  The  lumbar 
response  at  1  g  was  fairly  flat  between  2,  4,  6  and  8  Hz, 
dropping  off  at  11  Hz.  The  average  peak  occurred  at  6  and  8 

•  Hz,  with  5  subjects  showing  their  highest  response  at  6  Hz, 
and  the  other  five  at  8  Hz.  Magnitudes  of  response  were 
low,  however,  between  0.38  and  0.46.  At  the  thoracic  spine 
(Tl)  the  responses  were  greater  than  at  L2  for  all  of  the 
frequencies,  especially  at  the  lowest  frequencies.  The 

*  average  lumbar  response  to  a  1  g,  2  Hz  shock  was  0.4,  while 
the  thoracic  response  to  the  same  shocks  averaged  0.74. 
Thoracic  shocks  ranged  from  0.35  to  0.74.  The  trend  across 
frequencies  shows  clearer  dominance  at  the  2  Hz  frequency, 
with  a  decreasing  trend  across  frequencies  for  most 
subjects . 
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Internal  pressure 


Internal  pressure  response  to  positive  x-axis  shocks  at  the 
seat 

Figure  29  shows  a  typical  plot  of  the  internal  pressure 
response  to  positive  x-axis  shocks.  Data  for  x-axis 
internal  pressure  is  not  complete  for  all  subjects  due  to 
problems  with  the  pressure  transducers.  These  included 
slippage  of  the  pressure  probe,  and  not  having  a  probe  ^ 

available  due  to  damage  of  the  transducer  diaphragm.  Two 
subjects  were  omitted  from  the  3  g  and  2  g  data,  and  3 
subjects  from  the  1  g  data.  The  x-axis  pressure  data  are 
presented  in  Tables  38  to  40  and  Figure  30.  The  mean  delay 
between  peak  seat  acceleration  and  peak  internal  pressure 
response  ranged  from  0.04  to  0.08  seconds.  In  some  cases 
the  delay  was  negative,  indicating  that  the  subject  was 
responding  slightly  before  the  peak  acceleration  had 
occurred.  This  could  be  due  to  a  subject's  anticipation  of 
the  shock  due  to  the  initial  displacement  of  the  table  prior 
to  the  shock's  occurrence. 

The  3  g  pressure  response  closely  resembled  the  spinal 
response  to  3  g  impacts.  Six  of  eight  subjects  had  higher 
pressure  at  5  Hz,  compared  with  4  Hz,  giving  the  peak 
response  for  the  mean  of  the  subjects  at  5  Hz.  Three  of  the 
subjects  had  a  second  peak  at  8  Hz.  Ratios  of  pressure  to 
seat  acceleration  ranged  from  0.75  at  11  Hz  to  1.48  at  5  Hz. 

Levels  of  pressure  and  response  ratios  were  smaller  in  the 
X-axis  compared  to  the  z-axis. 

The  pressure  response  to  2  g  impacts  in  the  x-axis  also 
resembled  the  spinal  responses.  The  response  ratio  ranged 
from  0.9  to  2.1,  and  the  peak  mean  response  was  observed  at, 

6  Hz  (Table  39) . 

Internal  pressure  response  to  positive  y-axis  shocks  at  the 
seat 

Figure  31  shows  a  typical  plot  of  the  internal  pressure 
response  to  y-axis  shocks .  Data  for  the  y-axis  internal  ^ 

pressure  responses  were  less  problematic  than  the  x-axis. 

The  data  are  complete  except  for  one  subject  in  the  2  g  set. 

Delays  in  response  were  sometimes  negative,  but  less  often 

than  for  the  x-axis .  It  was  more  common  to  observe  a  double 

peak  in  internal  pressure  occurring  after  a  y-axis  seat  * 

impact  than  after  an  x-axis  shock. 

Data  for  pressure  responses  to  3 ,  2  and  1  g  y-axis  shocks 
are  found  in  Tables  41  to  43  and  Figure  32.  Levels  of 
pressure  and  response  ratios  were  lower  in  the  y-axis  than 
the  x-or  z-axis.  For  3  g  shocks,  ratios  ranged  from  0.29  at 
11  Hz,  to  0.86  at  4  Hz.  Individual  responses  varied  widely, 
however.  For  example,  responses  at  4  Hz  ranged  from  0.16  to 
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1.44  among  subjects.  The  trend  in  responses  to  3  g  shocks 
was  similar  to  acceleration  responses  at  the  spine.  The  4 
Hz  shock  produced  the  highest  average  level  of  pressure,  and 
was  dominant  for  7  of  10  subjects.  A  second  peak  at  8  Hz 
was  found  for  2  subjects. 

In  response  to  2  g  y-axis  shocks,  the  highest  response 
occurred  at  4  Hz.  However,  the  4  Hz  dominance  was  not  as 
prominent  for  internal  pressure  as  it  was  for  the  spinal 
acceleration  responses  to  2  g,  y-axis  shocks.  There  was 
little  difference  between  5  Hz  and  6  Hz  responses.  Ratios 
for  2  g  shocks  ranged  from  0.35  to  0.88. 

Pressure  responses  to  1  g  shocks  resembled  the  thoracic 
spinal  responses  most-  closely.  Pressures  were  highest  for  4 
Hz  shocks,  however  there  was  a  second  peak  at  8  Hz  for  6  of 
the  10  subjects,  and  in  the  mean  response.  Response  ratios 
ranged  from  0.4  to  1.3. 

Internal  pressure  response  to  negative  z-axis  shocks  at  the 
seat 

Figure  33  shows  a  typical  internal  pressure  response  to  a 
negative  z-axis  shock  (2  g  at  4  Hz) .  Average  delays  between 
the  seat  shock  and  internal  pressure  response  for  3  g  shocks 
ranged  from  0.16  seconds  at  4  Hz  to  0.12  seconds  at  11  Hz. 
For  2  g  shocks,  average  delays  in  pressure  response  were 
0.13  to  0.16  seconds.  For  1  g  shocks  the  delays  were 
similar,  from  0.10  to  0.19  seconds. 

The  magnitude  of  internal  pressure  response  for  a  3  g 
impact  in  the  z-axis  varied  from  approximately  30  mmHg  to 
230  mmHg.  By  contrast,  when  subjects  were  asked  to  cough, 
they  generated  internal  pressures  of  approximately  130  mmHg, 
laughing  produced  pressures  of  80-90  mmHg  and  bearing  down 
produced  instantaneous  pressures  of  225  mmHg  and  sustained 
pressures  of  130  mmHg.  Some  internal  pressure  responses  to 
3  g  impacts  were  nearing  levels  that  can  only  be  voluntarily 
induced  by  forceful  instantaneous  maneuvers. 

Tables  44  to  46  and  Figure  34  show  the  internal  pressure 
responses,  presented  as  ratios  of  pressure  in  mmHg  to  seat 
acceleration  in  m.s  ,  for  3,  2  and  1  g  negative  z-axis 
shocks,  respectively.  The  internal  pressure  responses 
closely  resembled  the  spinal  acceleration  responses .  At  3 
g,  the  dominant  response  was  to  4  Hz  shocks.  The  range 
among  subjects  was  large,  with  individual  ratios  varying 
from  2  to  12.4  at  4  Hz .  Average  ratios  for  3  g  shocks  were 
1.5  at  11  Hz,  increasing  to  7.6  at  4  Hz.  There  was 
considerable  variation  among  subjects  at  the  2  g  amplitude 
of  shock.  However,  the  average  responses  showed  a  similar 
trend.  Average  ratios  ranged  from  1.48  (11  Hz)  to  5.77  (4 
Hz) .  Ratios  were  lower  in  response  to  1  g  shocks  (average 
1.46  at  11  Hz  to  3.7  at  2  Hz).  The  response  at  2  Hz  (3.7) 


was  largest,  but  was  similar  to  the  4  Hz  response  (3.57). 
While  spinal  acceleration  showed  a  trend  toward  4  Hz 
dominance  for  1  g  shocks,  the  internal  pressure  response  was 
largest  at  2  Hz. 

Abdominal  and  chest  displacement 

Respitrace  response  to  positive  x-axis  shocks  at  the  seat 

Data  from  two  subjects  were  eliminated  in  the  x-axis  as 
no  peaks  were  detected.  A  typical  Respitrace  response  to  an 
x-axis  shock  is  shown  in  Figure  35.  Tables  47  to  49  and 
Figure  36  show  the  abdominal  responses  to  3,  2  and  1  g 
shocks  in  the  positive  x-axis.  At  3  g,  abdominal  response 
ratios  were  0.29  to  0.45  with  the  peak  occurring  at  5  Hz, 
and  a  flattening  at  8  Hz.  At  2  g  the  peak  occurred  at  4  Hz, 
and  at  1  g  at  2  Hz .  As  with  z-axis  responses,  the  responses 
were  larger  when  the  magnitude  of  the  shocks  was  smaller. 

Chest  response  ratios  for  x-axis  shocks  are  shown  for  7 
subjects  in  Tables  50  to  52  and  Figure  37.  As  with 
abdominal  displacement,  at  3  g  for  the  chest  the  dominant 
frequency  was  at  5  Hz  with  ratios  from  0.16  to  0.39.  At  2  g 
shock  inputs,  the  dominant  response  occurred  at  4  Hz  with 
ratios  from  0.19  to  0.45.  Similar  to  abdominal,  the  chest 
response  to  1  g  x-axis  shocks  was  dominant  at  2  Hz  and 
ratios  ranged  from  0.20  to  0.78. 

Respitrace  response  to  positive  y-axis  shocks  at  the  seat 

A  typical  Respitrace  response  to  a  y-axis  shock  is  shown 
in  Figure  38.  The  response  at  the  three  magnitudes  in  the 
y-axis  was  similar  with  a  dominant  frequency  at  4  Hz  and 
decreasing  response  with  increasing  frequencies.  Tables  53 
to  55  and  Figure  39  show  the  abdominal  response  ratios  to  1, 
2  and  3  g  shocks.  The  magnitude  of  response  was  again 
higher  for  1  g  shocks  (0.11  to  0.52)  than  2  g  shocks  (0.09 
to  0.31)  and  3  g  shocks  (0.07  to  0.28). 

Chest  response  ratios  are  shown  in  Tables  56  to  58  and 
Figure  40.  Chest  displacements,  and  response  ratios  were 
small.  Dominant  frequencies  were  again  at  4  Hz,  although 
for  2  g  shocks,  some  subjects  had  peak  responses  at  5  and  6 
Hz.  Response  ratios  were  0.05  to  0.20  for  3  g  shocks,  0.07 
to  0.28  for  2  g  shocks  and  0.06  to  0.26  for  1  g  shocks. 

Respitrace  response  to  negative  z-axis  shocks  at  the  seat 

Peaks  of  abdominal  and  chest  displacement  which  occurred 
in  conjunction  with  shocks  at  the  seat  were  as  high  as  30  mm 
for  individual  subjects.  Figure  41  shows  a  typical  tracing 
of  a  Respitrace  response  to  a  z-axis  shock.  Displacement  of 
the  abdomen  was  greater  than  the  chest.  In  many  cases, 
displacement  at  the  chest  was  negligible.  For  some 


subjects,  especially  in  the  z-axis,  no  peak  displacement  was 
detected.  It  could  either  be  that  there  was  no  peak  in 
response  to  the  shock  input,  or  the  peak  was  below  the  zero 
baseline  and  therefore  not  detected  by  the  analysis  program 
as  a  positive  peak. 

Tables  59  to  61  and  Figure  42  show  abdominal  responses, 
presented  as  ratios  of  displacement  in  mm  to  seat 
accleration  in  m.s"^  for  1,  2  and  3  g  shocks  in  the  negative 
z-axis.  Data  were  eliminated  for  4  of  the  10  subjects. 
Response  ratios  for  3  g  shocks  ranged  from  0.31  to  0.69  with 
the  peak  frequency  occurring  at  5  Hz ,  For  2  g  shocks  the 
peak  frequency  was  4  Hz  with  a  fairly  flat  response  at 
higher  frequencies.  The  range  at  2  g  was  from  0.26  to  0.74. 
A  4  Hz  peak  was  also  evident  for  1  g  shocks,  with  a  peak  at 
4  Hz .  The  pattern  for  2  g  and  1  g  shocks  resembled  the 
spinal  acceleration  responses.  Response  ratios  for  1  g 
shocks  were  generally  higher  than  those  for  higher 
magnitudes  of  shocks  (0.33  to  0.72).  Thus,  the  pattern  of 
response  magnitude  is  opposite  that  found  for  other 
measures . 

Tables  62  to  64  and  Figure  43  show  the  chest  response 
ratios  for  the  three  shock  magnitudes.  At  3  g  shock  inputs, 
the  dominant  frequency  was  4  Hz  (compared  with  5  Hz  for  the 
abdomen)  and  a  slight  peak  was  noticed  at  11  Hz.  Response 
ratios  ranged  from  0.19  to  0.43.  At  2  g  and  3  g,  the  peak 
response  was  similar  to  the  abdomen,  occurring  at  4  Hz.  The 
ratios  were  0.18  to  0.4  and  0.23  to  0.49  for  2  g  and  1  g, 
respectively. 

Electrocardiography 

Instantaneous  response  to  shocks 

Analysis  of  the  instantaneous  effect  of  3  g  shocks  on  R-R 
interval  is  not  complete.  However,  results  to  date  show  a 
strong  periodic  trend  throughout  the  data  at  a  frequency  of 
approximately  0.1  Hz.  This  rhythm  (thought  to  be  a 
baroreceptor  response)  may  have  overridden  other  changes  in 
R-R  interval.  Data  from  one  subject  suggests  that  a  shock 
induced  a  decrease  in  R-R  interval  approximately  6-7  beats 
after  the  shock.  Further  filtering  and  analysis  of  control 
data  needs  to  be  performed  to  determine  the  full  effect  of 
shocks  on  these  parameters. 

Fatigue  effects 

Results  of  the  ECG  signal  analysis  during  the  long  term 
experiments  are  sxammarized  in  Tables  65  to  67  and  Figures  44 
to  47 .  The  following  parameters  were  investigated  for 
evidence  of  fatigue  effects: 


Heart  rate 


T-wave  amplitude 
R-R  interval 

Spectral  variance  in  the  Baroreceptor  (Traube- 
Hering-Mayer )  band  (0.06-0.1  Hz) 

Spectral  variance  in  the  sinus  arrhythmia  band 
(0.12-0.4  Hz) 

Table  65  shows  the  results  of  the  2  hour  control  experiment 
for  all  four  subjects.  Table  66  shows  the  corresponding 
results  from  the  2  hour  exposure  to  3  g  z-axis  shocks  (1/2.5 
min)  and  Table  67  shows  the  results  from  the  2  hour  exposure 
to  2  g  z-axis  shocks  (32/min)  for  three  subjects. 

No  strong  fatigue  effects  were  apparent  during  shock 
exposure  or  when  compared  to  the  control  condition.  There 
was  no  consistent  difference  between  resting  and  recovery  T- 
wave  amplitudes  in  any  condition.  Similarly,  there  was  no 
discernible  consistent  trend  in  R-R  interval  between  or 
within  conditions.  The  R-R  interval  itself  did  vary, 
especially  around  the  shocks.  Some  of  this  change  is 
thought  to  be  movement  artifact  induced  by  the  shock  itself. 

However,  there  were  a  number  of  effects  worth  noting, 
especially  if  individual  subject  responses  are  examined. 
Subject  1  had  an  unusual  response  to  the  control  condition, 
especially  with  heart  rate,  R-R  interval  and  T-wave 
amplitude.  This  is  likely  due  to  two  factors.  Firstly, 
subject  1  was  the  first  one  to  undertake  the  long-duration 
experiments  and  therefore  may  have  experienced  some  anxiety. 
Secondly,  the  subject  had  an  infection  with  concomitant 
effects  on  his  physiological  response  (as  discussed  in  the 
Biochemistry  Results  Section) .  In  general,  the  3  g  z-axis 
shocks  (1/2.5  min)  produced  the  highest  levels  of  each 
measure  in  subject  1,  especially  from  the  middle  to  end 
portion  of  the  exposure.  There  was  only  data  for  both  the 
control  and  3  g  z-axis  (1/2.5  min)  condition  for  subject  2. 
Therefore  a  comparison  of  the  two  shock  exposures  could  not 
be  made.  The  results  of  Subject  3  showed  evidence  that  2  g 
z-axis  shocks  (32/min)  produced  higher  heart  rate  and  T-wave 
amplitudes  than  the  3  g  condition  as  well  as  a  higher 
response  than  the  control  condition.  Spectral  variance 
results  in  both  bands  were  variable  but  showed  a  trend 
towards  decreased  variance  in  Epoch  3  and  11  when  the 
performance  tests  were  being  undertaken.  Subject  4  shows  a 
similar  trend  to  Subject  3,  with  the  2  g  z-axis  shock 
condition  inducing  higher  heart  rate  and  T-wave  amplitude. 
The  spectral  analysis  results  show  no  particular  pattern. 

In  general,  2  g  z-axis  shocks  (32/min)  produced  a  higher 
heart  rate  than  the  control  or  3  g  z-axis  shocks  (1/2.5 
min) .  There  was  no  consistent  increase  in  heart  rate  over 
the  resting  level  when  subjects  were  exposed  to  the  shocks, 
with  the  exception  of  Subject  3  (2  g  z-axis  32  shocks/min) . 


Heart  rate  may  have  been  elevated  during  the  resting  phase 
due  to  an  anxiety  response  in  anticipation  of  the  shock 
exposure.  Had  each  subject's  response  been  corrected  for 
resting  level  then  the  results  may  have  been  clearer. 

The  spectral  variance  of  both  the  baroreceptor  and  sinus 
arrhythmia  bands  tended  to  decrease  during  the  epochs  where 
the  subjects  were  performing  the  Synthetic  Work  Battery, 
especially  in  the  shock  conditions.  There  were  a  number  of 
exceptions  to  this  pattern. 

Electromyography 

Muscle  response  to  shocks 

The  lumbar  EMG  burst  responses  to  both  positive  and 
negative  6  Hz  1  g  impacts  in  the  x,  y,  and  z  axes  were 
examined  to  determine  if  there  was  a  difference  in  either 
burst  onset,  burst  duration,  or  in  the  pattern  of  bursting 
versus  silent  periods.  Figures  48  to  53  provide  typical  L3 
EMG  responses  to  both  positive  and  negative  6  Hz  1  g  impacts 
in  each  axis. 

The  EMG  burst  response  to  an  impact  was  subject 
dependent,  with  some  subjects  demonstrating  a  single  clear 
EMG  burst  and  others  demonstrating  two  to  three  clear  bursts 
of  increased  muscle  activity  separated  by  silent  periods  of 
decreased  muscle  activity.  For  those  subjects  that 
demonstrated  multiple  bursts,  there  was  no  apparent 
difference  in  the  likelihood  of  observing  a  multiple  burst 
response  between  axes . 

The  orientation  of  the  impact  influenced  the  resultant 
pattern  of  silent  periods  and  bursts .  A  pattern  of  a  burst 
preceding  a  silent  period  (impact-burst-silent-burst-silent) 
was  observed  for  impacts  of  positive  orientation  in  the  x- 
axis  and  for  impacts  of  negative  orientation  in  the  z-axis. 
The  opposite  pattern  of  a  silent  period  preceding  a  burst 
(impact-silent-burst-silent-burst)  was  observed  for  the 
negative  x-axis  and  positive  z-axis  orientations  of  impact 
acceleration.  In  the  y-axis,  the  pattern  of  muscle  response 
was  dependent  on  the  relation  between  impact  direction  and 
the  side  of  the  body  from  which  EMG  was  monitored.  Negative 
y-axis  acceleration  (to  the  right)  resulted  in  a  burst- 
silent  pattern  for  the  right  lumbar  EMG  and  a  silent-burst 
pattern  for  the  left  Ixambar  EMG.  Positive  y-axis 
acceleration  (to  the  left)  resulted  in  the  opposite 
response,  with  a  burst-silent  pattern  on  the  left  and  a 
silent-burst  pattern  on  the  right.  In  the  silent-burst 
pattern  of  response  to  y-axis  impacts,  there  was  often  a 
very  brief  burst  that  preceded  the  first  silent  period.  In 
some  instances,  the  first  EMG  burst  or  silent  period 
appeared  to  be  initiated  prior  to  the  peak  impact 
acceleration  and  at  times  prior  to  the  onset  of  the  impact 


acceleration,  indicating  an  anticipatory  response  or  a 
response  to  some  other  stimuli. 

There  was  no  obvious  difference  in  the  duration  of  EMG 
bursts  between  axes  or  orientations  of  impacts.  Burst 
durations  ranged  from  40  msec  to  220  msec,  however  they  were 
more  typically  in  the  range  of  90  to  130  msec.  When  * 

multiple  bursts  were  evident,  the  second  burst  tended  to  be 
of  similar  magnitude  but  longer  duration  than  the  first 
burst.  The  duration  of  silent  periods  ranged  from  50  msec 
to  250  msec,  with  no  apparent  differences  between  axes  or  * 

impact  orientations . 

Figure  54  shows  the  average  magnitude  of  liimbar  muscle 
response  to  negative  3  g  z-axis,  and  positive  3  g  x  and  y 
axis  impact  accelerations  at  frequencies  of  4,  5,  6,  8,  and 
11  Hz  expressed  as  percent  of  maximal  voluntary  contraction 
(MVC) .  The  values  in  Figure  54  are  averages  of  two  impact 
responses  for  each  of  ten  subjects  at  each  impact  frequency. 

The  mean  response  to  x-axis  impacts  was  greater  than  the 
response  to  y-axis  impacts  at  all  frequencies  tested.  The 
mean  response  to  z-axis  impacts  at  4  to  6  Hz  was  greater 
than  that  of  either  x-axis  or  y-axis  responses,  but  the 
magnitude  of  the  z-axis  response  progressively  declines 
above  4  Hz .  There  was  evidence  of  a  resonant  peak  at  8  Hz 
in  the  x-axis,  and  to  a  lesser  extent  in  the  y-axis  muscle 
responses.  The  mean,  maxim;am,  and  standard  deviation  for 
the  EMG  response  to  impacts  in  each  axis  are  summarized  in 
Table  68.  The  magnitude  of  muscle  response  in  all  three 
axes  is  typically  less  than  8  percent  MVC.  Occasional 
bursts  of  30  to  65  percent  MVC  were  observed  in  response  to 
z-axis  impacts,  and  15  to  30  percent  MVC  in  x  and  y  axes. 

The  muscle  activity  and  internal  pressure  during 
coughing,  laughing,  and  bearing  down  maximally  were  measured 
to  provide  a  realistic  reference  to  compare  with  the  EMG  and 
internal  pressure  responses  to  impacts .  These  are 
illustrated  in  Figures  55  and  56.  Bearing  down  maximally 
produced  an  internal  pressure  of  approximately  250  mmHg. 

The  internal  pressure  seen  during  laughter  is  approximately 
thirty  percent  of  the  maximal  pressure  produced  during  the 
bearing  down  maneuver,  yet  the  magnitude  of  the  EMG  response  * 

is  not  appreciably  greater  for  bearing  down  than  laughing. 

Similarly,  coughing  produced  approximately  sixty  percent  of 

maximal  pressure  and  about  twice  the  pressure  of  laughter, 

yet  the  magnitudes  of  rectus  abdominus  and  erector  spinae  ^ 

EMG  were  not  appreciably  different  for  either  condition. 

Figure  57  illustrates  a  typical  response  to  a  positive  3  g 
z-axis  impact  for  both  lumbar  and  abdominal  muscle,  internal 
pressure,  and  liambar  spine  (L4)  acceleration.  The  magnitude 
of  a  typical  muscle  response  to  a  positive  3  g  z-axis  impact 
of  6  Hz  was  approximately  50  percent  of  the  magnitude 
measured  during  either  laughing  or  coughing,  yet  the  peak 
internal  pressure  was  similar  to  that  seen  during  laughter. 
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There  was  an  obvious  difference,  however,  in  the  duration  of 
EMG  burst  during  each  condition.  Bearing  down  produced  a 
sustained  contraction,  coughing  produced  a  relatively  long 
burst,  and  laughing  produced  a  very  brief  burst.  The 
duration  of  EMG  burst  during  laughter  was  similar  to  the 
muscle  response  to  a  3  g  shock. 

Muscle  fatigue 

The  mean  frequency  at  time  zero  (MFq) ,  determined  by 
extrapolation  from  a  time  series  of  EMG  mean  frequency 
during  sustained  isometric  back  extensions,  showed  no 
consistent  trend  between  pre-exposure  and  post-exposure 
values  for  any  of  the  conditions  studied.  The  range  of  MFq 
was  44  Hz  to  82  Hz,  with  a  mean  value  of  61.5  Hz  and  a 
standard  deviation  of  6.4  Hz.  For  the  most  severe  exposure 
condition  (2  g  z-axis  impacts  at  a  rate  of  32  per  minute)  an 
increase  in  MFq  greater  than  one  standard  deviation  was 
found  in  twenty-five  percent  of  sustained  contractions, 
while  twelve  percent  showed  a  significant  decrease  in  MFq 
(four  subjects,  measurement  of  four  muscles) .  The  remaining 
fatigue  trials  for  this  exposure  condition  showed  no  notable 
difference  between  MFq  before  and  after  exposure  to  the 
experimental  condition. 

The  mean  frequency  time  series  during  each  sustained 
contraction  rarely  showed  the  expected  exponential  decay 
reported  to  occur  during  a  sustained  submaximal  contraction 
in  the  literature  (Hagg  and  Suurkula  1991,  Deeb,  Drury  and 
Pendergast  1992,  Kadefors  et  al.,1978).  Instead,  the  mean 
frequency  showed  a  tendency  to  increase  with  time.  In  cases 
that  demonstrated  a  decrease  in  mean  frequency  with  time, 
the  time  constant  ranged  from  ten  to  sixteen  minutes.  Given 
the  low  number  of  cases  that  provided  meaningful  time 
constants,  and  the  brevity  of  test  duration  (30  seconds) 
relative  to  those  time  constants,  no  further  analysis  of 
time  constants  was  pursued. 

The  initial  lEMG  estimated  from  the  sustained  contraction 
lEMG  time  series  showed  no  difference  between  pre-exposure 
and  post-exposure  values  for  paraspinal  muscles  at  either  L3 
or  T9 .  The  mean  initial  lEMG  for  pre-exposure  trials  was 
0.6V  with  a  standard  deviation  of  0.3  V.  In  the  post¬ 
exposure  trials,  the  mean  initial  lEMG  values  for  the 
control  condition,  the  most  severe  condition  (two  hour,  2  g 
z-axis  impacts  at  32/min),  one  hour  2  g  x-axis,  and  one  hour 
2  g  y-axis  conditions  were  estimated  to  be  0.5  V,  0.4  V,  0.6 
V,  and  0.8V  respectively,  with  standard  deviations  of  0.3 
V,  0.4  V,  0.1  V,  and  0.5  V.  There  was  no  notable  difference 
between  subjects  with  respect  to  the  change  in  initial  lEMG. 

The  slope  of  the  lEMG  time  series  showed  no  consistent 
change  when  comparing  the  pre-exposure  and  post-exposure 
sustained  contraction  of  paraspinal  muscles  at  L3  and  T9 . 


The  mean  slope  for  the  pre-exposure  trials  was  1.5  mV/sec 
with  a  standard  deviation  of  4.9  mV/sec.  The  mean  of  the 
post-exposure  slope  for  the  control  condition  was  -0.5 
mV/sec  with  a  standard  deviation  of  6.1  mV/sec,  and  for  the 
most  severe  condition  (2  g  z-axis  impacts  at  32/min)  was  2.6 
mV/sec  with  a  standard  deviation  of  3.8  mV/sec.  One  hour 
exposures  to  2  g  x-axis  impacts  and  2  g  y-axis  impacts 
resulted  in  post-exposure  mean  slopes  of  1.9  mV/sec  and  1.6 
mV/sec  with  standard  deviations  of  4.9  mV/sec  and  3.8  mV/sec 
respectively.  No  single  subject  demonstrated  a  pre-to-post- 
exposure  increase  in  slope  for  all  four  paraspinal  muscle 
sites  monitored. 

The  lEMG  burst  response  to  impacts  was  compared  for  two 
subjects  during  the  first  and  last  epochs  of  the  most  severe 
condition  (2  g  z-axis  impacts  at  32/min) .  The  impact 
response  lEMG  was  expressed  as  muscle  activity  above 
background  activity  to  delineate  changes  in  the  background 
activity  versus  changes  in  the  burst  response.  The  average 
impact  lEMG  response  and  corresponding  background  lEMG  was 
measured  at  spinal  level  L3  and  is  summarized  in  Table  69 
for  subjects  3  and  4.  For  subject  3,  the  last  sampling 
interval  demonstrated  a  marked  increase  in  background  muscle 
activity  and  a  25  percent  increase  in  the  magnitude  of  the 
lEMG  response  to  2  g  z-axis  impacts.  For  subject  4,  there 
was  a  45  percent  increase  in  the  lEMG  response  to  impacts, 
however  no  change  in  the  background  lEMG  activity.  There 
was  a  change  in  the  timing  of  the  impact  response  between 
the  first  and  last  epochs,  however  the  two  subjects  examined 
demonstrated  contrasting  trends.  Subject  3  provided  a 
maximal  lEMG  response  when  the  burst  window  was  located  75 
msec  after  impact  in  the  first  epoch  and  50  msec  after 
impact  in  the  last  epoch,  indicating  that  the  muscle 
response  was  occurring  closer  to  the  impact  in  the  last 
epoch.  Burst  window  locations  of  25  msec  after  impacts  in 
the  first  epoch  and  50  msec  after  impacts  in  the  last  epoch 
provided  the  maximal  lEMG  response  for  subject  4,  indicating 
that  the  muscle  response  was  occurring  further  from  the 
impact  in  the  last  interval. 

Qptotrak  displacement  and  acceleration 

Preliminary  trials  with  the  Optotrak  indicated  that  it 
was  possible  to  collect  three  dimensional  data  from  6 
markers  at  a  sampling  rate  of  80  Hz.  Subsequent  analysis  of 
shock  records  revealed  occasional  missing  data  (drop  outs) 
in  the  Optotrak  data  files.  These  missing  data  tended  to 
coincide  with  the  shock  input  at  the  seat.  A  number  of 
Optotrak  shock  records  which  were  unaffected  by  missing  data 
were  analyzed,  and  sample  data  are  presented  in  Figures  58 
to  67  . 

The  displacement  data  measured  at  the  seat  in  the  x  and  y 
and  z  axes  are  shown  in  Figures  58,  59  and  60  respectively. 


The  displacements  represent  3  g  shocks  of  6  Hz,  5  Hz  and  4 
Hz  frequency  in  the  positive  x,  positive  y  and  negative  z 
directions.  The  corresponding  accelerations  were  calculated 
from  Optotrak  data  and  superimposed  on  each  displacement 
record.  Results  show  that  the  shocks  input  at  the  seat  are 
preceded  by  a . large  displacement  in  the  opposite  direction 
to  the  peak  acceleration  of  the  shock.  For  example,  in 
Figures  58  and  59,  a  negative  displacement  of  the  seat 
occurs  at  a  low  acceleration  magnitude.  The  shock  is  then 
generated  by  a  sudden  reversal  velocity,  producing  a  large 
positive  acceleration,  and  a  positive  displacement  of  the 
seat.  A  smaller  negative  acceleration  peak  is  generated  as 
the  seat  velocity  slows  and  it  returns  to  a  neutral 
position. 

The  seat  accelerations  calculated  from  Optotrak  data  are 
compared  with  the  corresponding  accelerations  measured  by 
the  seatpad  accelerometers  in  Figures  61,  62  and  63.  The 
Optotrak  acceleration  records  do  not  reproduce  the  higher 
frequency  oscillations  present  in  the  acceleration  data. 
However,  the  fundamental  waveform  and  magnitude  of  shocks 
produced  by  the  two  systems  are  very  similar  in  the  x  and  y 
axes.  This  is  shown  in  Figure  64,  in  which  the  calculated 
Optotrak  accelerations  and  seatpad  accelerometer  data  are 
superimposed.  Although  the  shock  waveforms  are  also  similar 
in  the  z-axis  (Figure  65) ,  the  Optotrak  record  does  not 
reproduce  the  higher  frequency  oscillations  shown  by  the 
seatpad  accelerometer.  These  oscillations  result  from  the 
subject  hitting  the  seat  after  a  negative  z-axis  shock,  and 
hence  are  not  part  of  the  shock  waveform  generated  by  the 
MARS.  As  the  unfiltered  accelerometer  data  showed  major 
frequency  components  of  100  to  200  Hz  (see  Figure  5) ,  the 
discrepancy  between  the  Optotrak  and  accelerometer  records 
results  from  their  different  sampling  rates  and  low  pass 
filter  characteristics. 

The  accelerations  of  the  liambar  spine  (L5)  ,  calculated 
from  Optotrak  data  in  response  to  a  positive  3  g  y-axis 
shock  and  a  negative  3  g  z-axis  shock  are  shown  in  Figures 
66  and  67.  The  corresponding  accelerometer  data  measured  at 
L3  (y-axis)  and  L4  (z-axis)  are  also  shown.  Accelerations 
of  the  lumbar  spine  obtained  by  Optotrak  and  accelerometer 
are  very  similar  in  magnitude  and  shape,  with  the  exception 
that  higher  frequency  components  are  recorded  in  the 
accelerometer  data.  The  accelerometer  recorded  a  slightly 
higher  peak  acceleration  at  the  spine  (L4) ,  than  the 
Optotrak  (L5) ,  in  response  to  a  3  g  z-axis  shock.  However, 
this  trend  was  not  consistent  across  axes  and  shock 
frequencies . 

Biochemistry 

A  graph  of  the  mean  of  each  biochemical  measure  (Figures 
68  to  97)  shows  the  trends  of  change  in  response  to  four 


long  term  exposure  experiments  in  which  biochemical  measures 

were  obtained  (Ll,  control;  L2,  rms  vibration  with  no 

shocks;  L6,  2  g  z-axis  shocks  at  32-''min;  L7 ,  Ig  shocks  in  x, 

y,  and  z  at  32 /min) .  In  some  variables  (ammonia,  creatine 

phosphokinase,  potassiiom  and  white  blood  cell  count)  where 

the  group  mean  data  appeared  to  mask  changes  observed  in  an  ^ 

individual  subject,  scatter  plots  of  individual  data  were 

also  plotted  (Figures  98  to  101) .  The  frequency 

distribution  of  urinary  data  coded  subjectively  is  tabled  in 

Table  7  0 . 

k 

The  coefficient  of  variation  (CV)  of  each  analysis 
reported  by  Allied  Clinical  Laboratories  (ACL) ,  and  tested 
by  BCR  is  reported  in  Appendix  C.  For  most  measures,  the  CV 
determined  by  BCR  was  similar  to,  or  less  than  that  reported 
by  ACL.  However,  some  of  the  variables  had  a  much  greater 
CV  than  reported  by  ACL,  making  these  data  difficult  to 
interpret.  Unacceptable  variation  was  noted  on  at  least  one 
occasion  in  the  measurement  of  the  following  variables: 
ammonia,  bilirubin,  catecholamines,  creatinine,  free 
hemoglobin,  lactate  dehydrogenase  iso-enzymes,  potassiiom  and 
von  Willebrand's  factor  antigen.  In  addition,  ACL  reported 
that  the  blood  was  hemolyzed  in  many  of  the  samples  analyzed 
for  free  hemoglobin,  rendering  this  analysis  useless. 

Detailed  records  made  by  USAARL  technicians  at  the  time  of 
the  experiments  indicated  that  only  a  few  samples  showed 
mild  hemolysis  prior  to  being  shipped  for  analysis. 

Catecholamine  analysis  provided  good  data  only  for 
norepinephrine.  Frequently  the  sample  had  to  be  diluted  for 
epinephrine  analysis,  and  qualitative  rather  than 
quantitative  data  were  reported. 

The  data  were  grouped  under  categories  to  isolate  more 
specifically  the  type  and  location  of  biochemical  stress  or 
organ  dysfunction.  The  categories  were:  stress,  muscle 
damage,  electrolyte  shift,  fluid  shift,  capillary  or  red 
blood  cell  damage,  hypoglycemia,  fatigue,  inflammation, 
joint  damage  or  bone  remodeling,  kidney  dysfunction,  and 
intestinal  bleeding. 

Stress  ^ 

Blood  cortisol  and  norepinephrine  were  used  as  a  measure 
of  stress.  Neither  of  these  variables  exceeded  normal 
values’ in  any  experimental  condition,  and  there  were  no 
differences  observed  between  experimental  conditions  4 

compared  to  control . 

Muscle  damage 

Elevated  serum  creatine  phosphokinase  (CPK)  is  a  marker 
for  extreme  muscle  exertion  or  damage .  In  the  control 
experiment  (Ll)  and  exposure  to  Ig  shocks  in  x,  y,  and  z  at 
32/min  (L7) ,  no  change  was  noted  in  the  group  mean  value  of 
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CPK.  Exposure  to  rms  vibration  with  no  shocks  {L2), 
resulted  in  a  trend  for  creatine  phosphokinase  (CPK)  to 
increase  from  12  to  36-hours  post-exposure.  This  trend  was 
clearly  evident  in  two  subjects  (subject  1  and  4) ,  and  could 
be  identified  in  a  third  subject  (subject  3) ,  In  two 
subjects  (subject  1  and  3)  the  peak  concentration  of  CPK  was 
elevated  above  the  normal  range  of  200  U/L.  Iso-enzyme 
analysis  indicated  that  only  the  iso-enzyme  for  skeletal 
muscle  (CPK-MM)  was  elevated.  This  suggests  that  skeletal 
and  not  cardiac  muscle  was  damaged  as  a  result  of  exposure 
to  the  experimental  condition.  Exposure  to  2g  z-axis  shocks 
at  32 /-min  (L6)  caused  the  group  mean  value  of  CPK  to 
increase  at  36-hours  post  exposure,  but  this  increase  was 
not  statistically  significant. 

In  subject  2  (L6,  2g  z-axis  shocks  at  32 /-min)  CPK  was 
elevated,  although  not  above  the  normal  range,  at  the  onset 
of  the  experiment.  This  was  followed  by  a  gradual  reduction 
in  CPK  at  24  and  36  hours. 

Other  measures  which  could  be  indicative  of  muscle  damage 
(blood  urea  nitrogen,  lactate  dehydrogenase,  serum 
myoglobin,  uric  acid,  urinary  myoglobin)  did  not  change 
significantly  as  a  result  of  the  experimental  exposures. 

Electrolyte  shift 

Electrolyte  shifts  during  exposures  were  assessed  by 
changes  in  plasma  calcium,  magnesium,  potassium,  and  sodium. 
A  trend  for  calcixim  to  decrease  in  two  experimental 
conditions  (L2,  rms  vibration  with  no  shocks;  L7,  Ig  shocks 
in  X,  y,  and  z  at  32 /min)  was  observed,  but  the  observed 
change  did  not  exceed  the  normal  value  for  calcium.  Neither 
magnesium  nor  sodium  changed  significantly  as  a  result  of 
experimental  exposure.  Plasma  potassium  tended  to  increase 
after  exposure,  to  vibration  with  no  shocks  (L2),  although 
the  mean  value  did  not  exceed  normal  values .  Examination  of 
individual  data  for  potassium  indicated  that  subject  2  had  a 
definite  increase  in  potassium  at  the  end  of  exposure 
(peak) ,  compared  to  the  pre-  or  mid-exposure  value.  However, 
given  the  amount  of  hemolysis  reported  by  the  laboratory, 
this  value  is  suspect  since  hemolysis  will  elevate  potassiiom 
concentration . 

Fluid  shift 

Change  in  plasma  volume  at  each  sampling  point  was 
calculated  from  a  difference  in  hematocrit  and  hemoglobin 
compared  to  the  pre-exposure  value.  There  was  a  high 
variability  in  this  measure  (-3.0%  to  +13.5%).  A  trend 
towards  an  increase  in  plasma  volume  at  12,  24,  and  36-hours 
in  all  experimental  conditions  was  observed,  but  the 
increase  was  not  statistically  significant. 


Capillary  and  red  blood  cell  damage 

No  consistent  trend  was  observed  in  von  Willebrand's 
factor  antigen.  Hence,  there  was  no  indication  of  capillary 
endothelial  damage.  Damage  to  red  blood  cells  was  not 
observed,  based  on  total  bilirubin  or  haptoglobin  data. 

Free  hemoglobin  data  could  not  be  considered,  since  many  of 
the  blood  samples  were  hemolyzed,  and  the  variability  in  the 
data  was  unacceptable . 

Hypog ly c  emi a 

Blood  glucose  tended  to  decrease  in  the  two  hours  of 
experimental  exposure  in  all  experimental  conditions, 
including  control.  The  change  in  blood  glucose  was  not 
precipitous,  and  did  not 'fall  below  a  normal  value  for 
plasma  glucose. 

Fatigue 

The  group  mean  data  for  blood  ammonia  and  lactate  did  not 
indicate  cumulative  fatigue  from  muscular  activity  as  a 
result  of  experimental  exposure.  The  decrease  in  lactate 
concentration  from  pre-exposure  to  the  mid-  and  peak- 
exposure  value  in  all  experimental  conditions  was 
unexpected.  The  group  mean  ammonia  concentration  tended  to 
increase  in  response  to  2g  z-axis  shocks  at  32-/min  {L6)  ,  but 
this  change  was  not  significant. 

Inflammation 

No  large  change  in  white  blood  cell  count  was  observed  in 
response  to  any  of  the  experimental  conditions.  Subjects  1 
and  3  had  a  pre-exposure  white  blood  cell  count  above  normal 
(i.e.  >10.5  X  1000  cells-mm^)  in  several  experimental 
conditions,  suggesting  a  mild  infection  independent  of  the 
current  experiment.  This  confounded  the  interpretation  of 
the  group  mean  data. 

Joint  damage  or  bone  remodeling 

Hydroxyproline  was  measured  from  a  single  urine  sample 
(not  a  24  hour  urine  collection) .  Neither  alkaline 
phosphatase  nor  urinary  hydroxyproline  showed  a  consistent 
change  which  might  be  linked  to  bone  or  joint  cartilage 
stress . 

Kidney,  bladder,  or  urinary  tract  dysfunction 

Urine  samples  were  tested  (pre-  and  peak-exposure)  for 
many  variables  which  could  suggest  dysfunction  of  the  kidney, 
bladder,  or  urinary  tract  (i.e.,  the  presence  of  blood, 
protein,  white  blood  cells,  red  blood  cells,  epithelial 
cells,  casts,  crystals,  mucous) .  Frequency  analysis  did  not 


show  any  consistent  change  in  these  variables  as  a  result  of 
exposure  to  vibration  and  shocks . 


Discussion 


Response  to  individual  shocks 

Although  the  control  signal  generated  by  B.C.  Research  for 
the  MARS  contained  shock  waveforms  consisting  of  a  damped 
sinusoid  of  a  single  frequency,  the  resultant  shock  output 
measured  at  the  seatpad  contained  higher  frequency  components . 
Resonances  of  >100  Hz  observed  at  the  seat  were  removed  by  the 
low  pass  filtering,  and  with  a  few  exceptions  were  not  seen  in 
the  spinal  acceleration  records.  However,  the  filtered  seat 
acceleration  record  still  contained  significant  frequency 
components  above  the  nominal  shock  frequency  of  the  control 
signal,  which  distorted  the  shock  waveform.  As  the  primary 
objective  was  to  examine  the  relationship  of  response  ratio  to 
shock  frequency,  these  higher  frequency  components 
superimposed  on  the  shock  waveform  must  also  lead  to  a 
distortion,  or  error  in  the  response  data  for  a  given 
frequency.  It  was  found  that  filtering  at  lower  frequency 
cut-offs  progressively  attenuated  the  magnitude  of  the  shock. 
In  addition,  as  the  shapes  and  hence  frequency  components  of 
the  spinal  acceleration  waveforms  differed  from  those  of  the 
seat,  a  reduction  in  the  magnitudes  of  the  seat  and  spinal 
acceleration  peaks  by  filtering  could  also  lead  to  erroneous 
results.  To  determine  the  sensitivity  of  the  results 
(response  ratio)  to  the  cut-off  frequency  of  a  low  pass 
filter,  examples  of  seat  and  spinal  acceleration  data  were 
analyzed  in  an  unfiltered  condition,  and  with  a  low  pass 
filter  set  at  twice  the  nominal  shock  frequency  of  the 
individual  shocks  (8,  10,  12,  16  and  22  Hz  cut-off) .  A 
comparison  of  these  data  processing  techniques  is  shown  in 
Figure  102  for  3,  2  and  1  g  shocks  in  the  x-axis.  Although 
there  are  some  differences  in  magnitude  of  the  response  ratios 
between  the  two  methods,  they  generally  are  not  large,  and  the 
relationship  of  response  ratio  to  shock  frequency  at  each 
shock  amplitude  remains  essentially  unchanged.  Hence,  it  is 
considered  unlikely  that  the  study  substantially  altered  the 
conclusions  regarding  the  relationship  of  spinal  acceleration 
response  ratio  to  shock  frequency. 

Vertebra-skin  transfer  function 

The  calculation  of  vertebra-skin  transfer  functions  for 
the  accelerometer-tissue  sub-system  was  particularly 
difficult.  It  was  noted  that  there  was  often  an 
initial  very  high  frequency  response  on  release  of  the 
skin,  followed  by  lower  frequency  oscillations.  The 
values  of  05^  and  ^  calculated  for  any  given 


accelerometer  showed  a  large  variance  between 
successive  trials,  and  within  a  response,  depending  on 
the  particular  peaks  digitized.  As  the  natural 
frequencies  calculated  (30  to  80  Hz)  were  well  above 
the  nominal  shock  frequencies,  vertebra-skin  transfer 
function  had  no  significant  effect  on  the  spinal 
accelerations  measured  in  the  response  to  the  x-  and  y- 
axis  shocks .  In  the  case  of  the  negative  z-axis 
shocks,  the  frequency  of  the  spinal  acceleration 
approached  that  of  the  transfer  function,  particularly 
for  -3  g  shocks  at  4  Hz.  Thus,  the  transformed 
acceleration  peak  and  the  resultant  response  ratio  was 
sensitive  to  the  value  of  05^.  In  some  cases 
application  of  the  inverse  transfer  function  attenuated 
the  filtered  acceleration  peak.  In  a  few  cases 
however,  the  acceleration  peak  was  amplified  after 
application  of  the  inverse  transfer  function.  By  first 
filtering  the  raw  acceleration  data  at  a  low-pass  cut¬ 
off  frequency  of  60  Hz,  most  of  the  high  frequency 
components  were  removed,  and  the  ability  of  the  inverse 
transfer  function  to  amplify  the  acceleration  was  thus 
limited.  However,  the  results  obtained  in  this  study 
suggest  that  the  ass\imption  of  a  simple  Kelvin  element 
for  the  accelerometer-tissue  sub-system  may  be 
inappropriate  for  the  analysis  of  spinal  acceleration 
in  response  to  shocks .  It  has  been  shown  that  the  skin 
displays  non-linear  properties  and  that  it  behaves 
differently  at  high  frequencies  (Hinz  et  al.,  1988). 

Evidence  of  this  was  found  in  the  acceleration  data 
obtained  when  the  skin  was  pulled  and  released.  There 
was  also  evidence  of  high  frequency  vibration  at  the 
accelerometer  in  some  subjects  when  they  struck  the 
seat  following  a  negative  z-axis  shock.  These 
oscillations  may  not  have  been  transmitted  from  the 
vertebra  to  the  skin  surface,  but  rather  have  resulted 
from  a  secondary  resonance  at  the  skin  surface. 

Without  filtering  of  the  data  to  remove  high  frequency 
resonances ,  the  inverse  transfer  function  derived  from 
a  simple  Kelvin  Element  would  severely  distort  the 
accelerations  calculated  at  the  vertebra.  From  the 
limited  information  available  from  this  study  it  is 
concluded  that  without  a  more  complex  model,  the 
correction  of  accelerations  for  skin  movement  probably 
does  not  improve  accuracy,  particularly  for  shocks 
where  the  subject  leaves  and  then  strikes  the  seat. 

Spine  acceleration,  internal  pressure  and  Respitrace 
responses  to  shocks 

Despite  the  range  of  responses  between  individuals,  the 
pattern  of  mean  responses  for  the  various  measures  of  spinal 
acceleration  (in  addition  to  internal  pressure  and 
Respitrace)  are  remarkably  similar.  There  were,  however. 


differences  in  magnitude  and  shape  of  acceleration  responses 
among  the  axes  (x,  y  and  z)  and,  in  some  cases,  among  the 
shock  amplitudes.  The  largest  spinal  responses  were  to  -z- 
axis  shocks.  Responses  in  the  -z-axis  were  not,  however, 
due  to  the  shock  input,  but  rather  to  the  subject  impacting 
the  seat.  Acceleration  responses  to  +y-axis  shocks  were 
higher  in  magnitude  than  to  +x-axis  shocks. 

The  shape  of  the  spinal  acceleration  response  in  the  +x- 
axis  across  various  frequencies  of  shocks  showed  a  non¬ 
linearity  with  magnitude  of  shocks.  At  3  g,  there  was  a 
clear  peak  in  response  at  5  Hz,  both  in  the  thoracic  and 
lumbar  spine,  and  both  in  the  x-and  z-axis  response  to  the 
x-axis  shock  input.  There  was  also  a  smaller  peak 
noticeable  at  8  Hz .  At  2  g  and  1  g,  however,  the  peak 
response  shifted  to  4  Hz,  and  there  was  no  second  peak  at  8 
Hz.  This  pattern  is  clearly  different  from  that  suggested 
for  sinusoidal  vibration  by  the  ISO  2631  (1982) .  The  ISO 
curve  suggests  the  dominant  frequency  is  1-2  Hz,  with  a 
tapering  of  response  beyond  2  Hz.  There  was  also  a  non¬ 
linearity  of  response  between  the  various  amplitudes  of 
shocks .  Transmission  of  shock  peaks  was  highest  for  3  g 
shocks,  followed  by  2  g  and  finally  1  g  shocks. 

Interestingly,  the  z-axis  spinal  response  to  x-axis  shock 
inputs  was  higher  than  the  x-axis  response.  This  is  likely 
due  to  rotation  of  the  vertebrae  and  whiplash-like  effects. 
As  a  subject  flexes  in  response  to  a  shock,  there  are  x-and 
z-axis  components  due  to  curvature  of  the  spine. 

In  the  y-axis,  the  shape  of  the  spinal  acceleration 
response  curves  were  similar  for  all  amplitudes  of  shocks, 
for  lumbar  and  thoracic  locations,  and  for  both  y-and  z-axis 
responses .  The  non-linearity  in  shape  and  magnitude  found 
for  x-axis  shocks  was  not  apparent  for  y-axis  shocks.  The 
dominant  frequency  response  was  at  the  lowest  measurable 
frequency  (4  Hz  for  3  g  and  2  g  and  2  Hz  for  1  g) .  It  is 
impossible  to  know  whether  the  truly  dominant  frequency  is 
lower  than  these,  since  lower  shock  frequencies  could  not  be 
measured.  The  ISO  2631  suggests  that  for  sinusoidal 
vibration  the  dominant  frequency  for  the  y-axis  is  1-2  Hz. 

We  cannot  tell  from  this  data  whether  the  same  is  true  for 
response  to  shocks.  The  z-axis  spinal  component  in  response 
to  a  y-axis  shock  was  very  small. 

In  the  -z-axis,  the  shape  of  the  spinal  acceleration 
response  curves  for  3  g  and  2  g  shocks  resembled  that  of  the 
y-axis,  with  dominance  at  4  Hz  and  a  decrease  in  response 
with  increasing  frequency.  Again,  it  is  impossible  to  know 
whether  a  lower  frequency  shock  will  produce  a  higher 
response.  For  1  g  shocks,  however,  the  dominant  frequency 
was  4  Hz,  and  not  2  Hz.  The  ISO  2631  suggests  that  for 
sinusoidal  vibration  the  most  sensitive  frequencies  are  from 
4-8  Hz.  The  British  Standards  (BS  6841,  1987),  which  are 
meant  to  accommodate  shocks,  suggest  that  frequencies  as 


high  as  10  Hz  are  important  in  the  z-axis.  The  DRI, 
designed  for  single  large  +z-axis  impacts,  suggest  the 
dominant  frequency  is  11  Hz .  Results  from  this  study  with  - 
z-axis  shocks,  suggest  the  dominant  frequency  is  lower  than 
that  in  current  standards .  A  non-linearity  in  response  to 
shock  magnitude  was  apparent  in  the  -z-axis,  with  highest 
responses  recorded  for  higher  amplitude  shocks. 

The  x-axis  response  at  the  spine  to  z-axis  shocks  at  the 
seat  were  of  similar  magnitude  and  shape  to  the  spinal  z- 
axis  responses.  As  with  x-axis  shocks,  the  z-axis  shocks 
caused  a  whiplash-like  effect  in  the  body  and  forward 
rotation  accompanied  the  shocks.  At  the  lowest  frequencies, 
the  x-and  z-axis  response  at  the  lumbar  spine  was,  on 
average,  2.5-times  the  seat  acceleration.  This  large 
amplification  of  input  shock  at  the  lumbar  spine  is  reduced 
in  magnitude  by  one-half  in  the  upper  thoracic  spine. 

Clearly,  the  body  is  a  complex  system  that  does  not  respond 
as  a  single  mass  to  shock  inputs. 

When  the  subject  impacts  the  seat  following  a  -z-axis 
shock,  the  seat  accelerometer  registers  a  very  high- 
frequency  component  (approximately  150  Hz) .  A  high- 
frequency  response  (30-90  Hz)  is  transmitted  to  the 
accelerometers  at  the  spinal  vertebrae.  Data  in  the 
literature  suggests  that  the  body  damps  these  high  frequency 
components  (Fairley  and  Griffin,  1989) .  Our  results  showed 
higher  spinal  responses  to  2  and  4  Hz  shocks,  compared  with 
11  Hz  shocks.  However,  the  much  higher  frequencies  of  30-90 
Hz  were  not  being  damped  by  the  spine.  It  might  be 
suggested  that  these  high  frequency  components  are  skin 
movement.  However,  the  vertebra-skin  transfer  function  did 
not  remove  these.  In  addition,  examination  of  the  internal 
pressure  responses  also  showed  similar  high-frequency 
components .  It  may  be  that  for  a  large  enough  impact  the 
body  is  responding  as  a  non-linear  system  with  high 
frequency  components  being  transmitted. 

Results  from  measures  of  internal  pressure  in  response  to 
various  frequencies,  amplitudes  and  directions  of  shocks  are 
remarkably  similar  to  spinal  acceleration  responses . 

Highest  internal  pressure  responses  were  measured  in  the  z- 
axis .  Some  internal  pressure  magnitudes  produced  by  3  g 
shocks  exceeded  pressures  that  could  be  produced  by  subjects 
through  forceful  instantaneous  maneuvers.  Lowest  magnitudes 
were  for  y-axis  shocks,  however  often  double  peaks  were 
observed  in  the  y-axis , 

Due  to  time  constraints  in  the  pilot  study,  responses  to 
shocks  in  the  three  axes  were  measured  in  one  direction  for 
each  axes:  +x,  +y  and  -z.  Since  the  body  is  not 
symmetrical  about  the  x  and  z-axes,  there  will  likely  be 
differences  in  response  depending  upon  direction  of  the 
shock.  For  example,  the  body  responds  to  a  -z-axis  shock  by 


leaving  the  seat.  The  body's  acceleration  is  limited  to  1 
g,  but  a  high-frequency  event  that  is  different  from  the 
seat  shock  wave  occurs  as  the  subject  hits  the  seat.  In 
response  to  a  +z-axis,  the  subject  would  not  likely  leave 
the  seat,  but  the  spine  would  undergo  compression.  Cursory 
examination  of  spinal  responses  to  +z-axis  shocks  from  the 
long-term  experiments  showed  differences  that  should  be 
explored  in  future  experiments.  To  complete  the  data  set, 
short-term  experiments  should  be  conducted  in  the  -x  and  +z 
directions . 

Limitations  in  the  excursion  capability  of  the  MARS 
shaker  table  meant  that  very  low  frequency,  high  amplitude 
shocks  could  not  be  investigated.  At  3  g,  the  lowest 
frequency  of  shock  that  could  be  reproduced  was  4  Hz,  and  at 
1  g,  a  2  Hz  shock  was  the  lowest  frequency  possible.  With 
the  results  obtained,  in  most  cases  it  is  not  possible  to 
determine  whether  4  Hz  at  3  and  2  g  or  2  Hz  at  1  g 
represents  the  true  dominant  frequency  response,  since  lower 
frequencies  could  not  be  tested. 

Results  from  these  experiments  highlight  limitations  of 
current  standards  for  vibration  and  shock.  New  standards 
for  shocks  should  consider  the  non-linearity  of  the  body, 
which  may  result  in  curves  that  are  different  for  not  only 
different  axes  (x,  y  and  z) ,  but  different  directions  in 
each  axes  (+  and  -)  and  for  different  magnitudes  of  shocks. 
The  curve  may  be  a  totally  different  shape  for  very  high 
magnitude  shock,  compared  with  less  severe  shocks. 

Assessment  of  ECG  parameters  as  valid  measures 

Heart  rate  did  show  differences  between  experimental 
conditions.  However,  it  did  not  show  any  strong  trends  over 
time  in  response  to  shocks.  It  is  affected  by  many 
different  influences  including  the  psychological  and 
physiological  state  of  the  subject.  Thus  it  is  not  a 
reliable  index  for  investigating  fatigue  with  these  levels 
of  stressor. 

T-wave  amplitude  did  not  show  the  expected  decrease 
during  the  recovery  phase  compared  to  the  control  or  resting 
condition.  It  is  possible  that  there  was  an  inadequate 
recovery  period  to  fully  estimate  the  recovery  response  of 
the  subjects.  Conversely,  exposure  to  these  magnitudes  of 
shocks  may  not  have  been  long  enough,  or  the  magnitude  and 
shocks  per  minute  great  enough  to  induce  fatigue.  The  time 
periods  used  here  were  significantly  shorter  than  many 
operational  missions. 

It  was  anticipated  that  the  shock  exposure  would  produce 
a  narrowing  in  variance  of  the  R-R  spectral  components 
compared  with  the  control  and  resting  conditions .  This  was 
not  apparent  in  the  results.  There  were  large,  but 


seemingly  random  variations  in  the  two  frequency  bands  under 
consideration.  The  levels  of  variance  are  also  higher  than 
normally  reported  in  the  literature  as  the  data  points  are  a 
sum  over  150  seconds  (rest  and  recovery)  or  270  seconds 
(epochs)  and  in  some  cases  600  seconds  in  the  recovery 
phase.  A  more  standard  sampling  period  would  be  30  seconds 
and  the  natural  logarithm  of  the  variance  is  usually 
expressed. 

The  effect  of  undertaking  the  performance  task  on 
spectral  analysis  of  the  R-R  interval  produced  a  decrease  in 
variance,  as  expected.  This  effect  is  associated  with 
increased  effort  and  psychological  stress.  It  is  clear  that 
undertaking  cognitive  tasks  has  the  potential  to  affect  the 
physiological  response  which  may  compound  any  effect  induced 
by  the  shocks.  Since  TGV  crews  will  likely  be  involved  to  a 
greater  or  lesser  extent  in  cognitive  tasks,  a  simulation  of 
this  stress  should  be  included  in  any  future  experiments . 

EGG  response  to  shocks 

There  is  some  indication  from  the  results  that  2  g  z-axis 
shocks  (32/min)  are  more  stressful  than  3  g  z-axis  (1/2.5 
min).  However,  Subject  I's  results  are  in  doubt  due  to 
illness  during  the  control  experiment  and  possible  carry 
over  to  other  days,  and  data  was  only  available  for  two 
conditions  for  Subject  2.  There  are,  therefore,  only  two 
subjects  (3  &  4)  on  which  to  base  these  observations.  It  is 
clear  that  with  such  small  ntimbers,  no  conclusive  statements 
can  be  made.  Further  analysis  and  manipulation  of  the  data 
as  indicated  in  the  results  section  may  well  clarify  some  of 
the  outstanding  questions  and  issues. 

Electromyography 

Muscle  response  to  shocks 

The  main  purpose  of  the  muscle  response  to  an  impact 
acceleration  is  stabilization  of  the  spine  and  upper  body  to 
minimize  motion  and  preserve  the  seated  posture.  This  will 
maintain  the  individual's  ability  to  continue  performing  the 
task  at  hand,  as  well  as  minimize  the  possibility  of 
secondary  impacts  (through  loss  of  balance  or  striking 
surfaces)  that  may  prove  harmful.  The  pattern  of  muscle 
response  observed  as  bursts  and  silent  periods  in  the  EMG 
can  help  to  elucidate  the  strategy  used  to  compensate  for 
the  imposed  motion  of  an  impact.  Evidence  of  two  phases  in 
the  muscle  response  can  be  seen  in  the  EMG  response  to  z- 
axis  and  y-axis  impacts.  The  first  phase  is  one  of 
stabilization  in  which  there  is  co-contraction  of 
antagonistic  muscle  groups  (erector  spinae  and  rectus 
abdominus  for  z-axis,  and  both  right  and  left  erector  spinae 
for  y-axis) .  The  second  phase  is  one  of  positional 
correction  in  which  there  is  reciprocal  activity  of 


antagonistic  muscle  groups  to  correct  for  the  imposed  motion 
of  the  upper  trunk  and  to  restore  neutral  posture. 

In  contrast  to  the  early  findings  of  Cursiter  and  Harding 
(1974)  of  alternating  reciprocal  activity  of  paraspinal  and 
abdominal  muscle  in  response  to  z-axis  sinusoidal  vibration, 
a  z-axis  impact  acceleration  initially  produced  coincidental 
bursting  of  lumbar  erector  spinae  and  rectus  abdominus  which 
occurred  at  the  same  time  as  the  peak  internal  pressure  (see 
Figure  57) .  A  similar  pattern  was  also  seen  in  the  response 
of  right  and  left  paraspinal  muscles  to  a  y-axis  impact,  in 
which  there  was  often  a  brief  burst  of  activity  prior  to  the 
silent  phase  from  the  muscle  that  initially  undergoes 
shortening.  This  phenomenon  can  be  seen  in  the  right  lumbar 
muscle  EMG  in  response  to  a  positive  y-axis  impact  in  Figure 
50,  and  in  the  left  lumbar  muscle  response  to  a  negative 
impact  in  Figure  51.  Co-contraction  of  antagonistic  muscles 
will  stiffen  the  trunk  and  increase  the  resistance  to  motion 
in  any  direction.  It  will  also  contribute  to  an  increase  in 
pressure  within  the  abdominal  cavity.  This  may  be  an 
initial  response  that  minimizes  the  subsequent  correction 
required  for  positional  deviation  due  to  the  imposed  motion. 
The  increase  in  abdominal  pressure  will  distribute  the  force 
of  the  impact  across  a  larger  area  of  tissue,  since  it  will 
increase  the  likelihood  of  transmission  within  the  gut. 

This  may  decrease  stress  on  the  spinal  structures  but  result 
in  greater  stress  on  the  viscera  as  they  experience  a  larger 
pressure  transient. 

In  x-axis  and  z-axis  impacts,  the  paraspinal  muscle 
response  is  bilaterally  symmetrical  in  both  timing  and 
magnitude.  This  is  not  surprising  given  that  the  primary 
resultant  motions  will  be  fore-aft  and  spinal  flexion- 
extension.  The  relative  pattern  of  response  to  positive  and 
negative  x-axis  impacts  (Figures  48  and  49)  demonstrates 
that  a  paraspinal  EMG  burst  occurs  to  compensate  for  forward 
flexion,  while  a  silent  period  during  backward  extension 
minimizes  opposition  to  abdominal  muscle  activity  which 
resists  the  extension.  A  positive  x-axis  shock  would  result 
in  an  initial  backwards  displacement  of  the  seat,  which 
causes  flexion  of  the  trunk  due  to  the  inertia  of  the  upper 
body.  The  compensatory  response  is  therefore  an  initial 
increase  in  activity  of  the  back  muscles  to  restore  upright 
posture.  A  negative  x-axis  shock  produces  the  opposite 
displacement  and  compensatory  pattern  of  muscle  response. 

The  restoring  force  would  be  provided  by  the  abdominal 
muscles  while  the  activity  of  the  back  muscles  would 
decrease.  The  pattern  of  response  to  a  positive  z-axis 
impact  suggests  that  there  is  an  initial  back  extension  as 
the  seat  moves  downwards,  followed  by  flexion  as  the  seat 
returns  to  the  neutral  position.  This  results  in  an  EMG 
pattern  where  a  silent  period  precedes  a  burst  of  erector 
spinae  activity.  A  negative  z-axis  impact  results  in  an 


initial  spinal  flexion  followed  by  extension,  resulting  in 
an  EMG  pattern  in  which  a  burst  precedes  a  silent  period. 

In  response  to  y-axis  impacts,  there  is  a  reciprocal 
stretching  and  shortening  of  right  and  left  paraspinal 
muscle  groups  since  the  motion  is  in  a  side-to-side 
orientation.  The  pattern  of  muscle  response  to  y-axis 
impacts  is  therefore  dependent  on  the  side  of  the  body  that 
is  monitored  relative  to  the  orientation  of  the  impact.  The 
actual  motion  of  the  MARS  produced  a  small  displacement  in 
the  direction  opposite  that  required  for  the  intended  impact 
acceleration.  The  pattern  of  muscle  response  to  y-axis 
shocks  suggests  that  this  small  displacement  was  sufficient 
to  elicit  a  response  from  the  paraspinal  muscles.  Failure 
to  see  a  similar  response  in  either  of  the  x-  or  z-axes 
suggests  that  there  is  an  greater  sensitivity  to  motion  in 
the  y-axis,  yet  the  magnitude  of  muscle  response  is  lower  in 
the  y-axis  than  x-  or  z-axes. 

It  was  found  that  the  magnitude  of  muscle  response  was 
dependent  on  shock  frequency  for  each  axis  (see  Figure  54) . 
These  data  support  the  pattern  of  frequency-dependence  found 
in  shock  transmission  data.  It  further  suggests  that 
frequencies  below  4  Hz  may  elicit  larger  muscular  responses, 
particularly  in  the  z-axis,  which  should  be  examined  in 
future  studies . 

EMG-force  calibration 

The  force  calibration  procedure  that  was  employed  in 
these  experiments  may  be  improved  by  eliminating  measurement 
of  resting  EMG  and  by  increasing  the  number  of  measurements 
in  the  range  of  2  to  10  percent  MVC.  Resting  EMG  data 
proved  to  be  larger  than  that  measured  during  a  5  percent 
MVC  contraction  in  most  cases.  This  was  most  likely  due  to 
the  method  of  defining  rest  as  zero  force  measured  at  the 
load  cell.  This  condition  may  be  achieved  by  leaning 
forward  slightly,  which  increases  the  muscle  force  required 
to  support  the  mass  of  the  upper  torso. 

Fatigue 

The  mean  power  frequency  (MF)  and  lEMG  parameters 
measured  during  sustained  isometric  contractions  before  and 
after  the  long  term  exposures  showed  no  consistent  trends 
indicative  of  localized  muscle  fatigue.  This  is  probably 
because  most  of  the  isometric  contractions  failed  to  produce 
an  EMG  response  that  typifies  muscle  fatigue  within  the  30 
second  fatigue  trial  itself.  Many  of  the  sustained 
contractions  produced  an  increasing  mean  frequency  rather 
than  the  expected  decrease  that  signifies  localized  muscle 
fatigue.  This  phenomenon  has  been  reported  elsewhere  by 
Hagg  and  Suurkula  (1991) .  They  showed  that  an  increasing  MF 
or  zero-crossing  rate  was  positively  correlated  with  work- 
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related  myalgia  of  infraspinatus  and  descending  trapezius 
muscles.  Similarly,  lEMG  magnitude  tended  to  decrease 
during  some  sustained  contractions,  rather  than  increase  as 
expected  if  the  contraction  induced  localized  muscle  fatigue 
and  the  associated  recruitment  and  increased  firing  rate  of 
active  motor  units.  This  has  been  reported  elsewhere 
(Seidel,  Beyer  and  Brauer  1987,  Chapman  and  Troup  1970)  for 
sustained  contractions  of  liimbar  back  muscles  at  20  to  40 
percent  of  MVC.  However,  lEMG  was  usually  found  to  increase 
for  contractions  at  50  percent  MVC  by  Troup  and  Chapman 
(1972)  .  Observation  of  the  trends  in  lEMG  and  mean 
frequency  in  the  present  study  showed  that  the  relative 
activity  of  lumbar  and  thoracic  paraspinal  muscles  was  not 
complementary  or  consistent;  however,  there  was  a  similar 
trend  bilaterally  at  each  level. 

The  failure  to  produce  consistent  or  measurable  muscle 
fatigue  during  test  contractions  in  the  present  study  may  be 
attributed  to  the  complexity  of  the  musculature  of  the  back, 
inadequate  control  over  subjects'  posture  and  muscle 
temperature,  or  to  the  selection  of  a  sustained  force  that 
failed  to  produce  significant  fatigue  within  30  seconds. 

Most  subjects,  however,  reported  some  sensation  of  back 
muscle  fatigue  following  the  30  second  pulls,  and  a  50% 
maximal  contraction  required  a  visible  effort  to  sustain. 
Although  sitting  posture  was  reproduced  as  closely  as 
possible  in  both  trials,  subjects  were  allowed  to  brace 
against  the  floor  using  their  legs,  which  may  have  allowed 
for  changes  in  posture  and  muscle  utilization  during  the 
static  pulls.  Because  of  the  very  short  moment  arm  of  back 
muscles,  it  is  possible  to  produce  large  changes  in  joint 
torque  with  relatively  small  changes  in  the  length  of  the  . 
moment  arm.  It  is  difficult  to  control  the  angle  of  pelvic 
tilt  and  spinal  curvature,  which  could  significantly 
influence  the  length  and  activity  of  muscles  at  the  lumbar 
level.  It  is  also  difficult  to  control  the  posture  of  the 
upper  trunk,  which  will  influence  the  muscle  moment  arm,  and 
length  and  activity  of  muscles  in  the  thoracic  region.  The 
cargo  straps  that  were  used  to  secure  a  subject's  legs 
provided  a  degree  of  control  over  the  pelvis,  but  still 
allowed  for  considerable  alteration  in  pelvic  tilt.  A  means 
of  directly  securing  the  pelvis  would  greatly  improve 
postural  control. 

Additionally,  fatigue  may  have  been  induced  in  muscles 
other  than  those  monitored  by  EMG.  Monitoring  EMG  from  a 
greater  number  of  muscle  sites  during  a  sustained 
contraction  would  help  to  clarify  the  relative  recruitment 
of  muscles  involved  in  force  generation  during  this 
maneuver.  The  time  constant  for  mean  frequency  decay  was 
calculated  to  be  about  10  to  16  minutes  for  cases  that 
demonstrated  a  decrease  in  mean  frequency  during  sustained 
contractions  of  50  percent  MVC.  This  indicates  that  the 
duration  of  the  sustained  trials  would  need  to  be  increased 


significantly  to  provide  a  reliable  measure  of  fatigue  if  a 
contraction  level  of  50  percent  MVC  is  used. 

Another  consideration  in  assessing  muscle  fatigue  is  the 
fact  that  the  muscle  response  to  impacts  was  typically  less 
than  ten  percent  of  a  maximal  voluntary  contraction  and 
lasted  less  than  200  milliseconds.  At  this  level  of 
intermittent  dynamic  contraction,  one  would  not  expect  an 
appreciable  amount  of  localized  muscle  fatigue  unless 
exposure  to  the  motion  environment  was  continued  for  a  very 
long  period  of  time.  At  low  levels  of  contraction,  the 
motor  units  most  likely  to  experience  fatigue  are  the 
smaller  oxidative  fibres.  A  sustained  isometric  fatigue 
trial  at  a  higher  activity  level,  such  as  50%  MVC,  is  likely 
to  recruit  larger  fibers  and  a  greater  number  of  motor  units 
than  the  dynamic  burst  response  to  an  impact.  Localized 
fatigue  of  the  smaller  fibers  induced  by  the  motion 
environment  may  therefore  be  masked.  Fatigue  trials  in 
future  studies  would  better  represent  the  muscle  fibers  of 
interest  and  be  more  sensitive  to  muscle  fatigue  if  the 
level  of  contraction  in  the  fatigue  trial  approximates  the 
level  of  contraction  observed  in  response  to  the  motion 
environment .  Mean  power  frequency  could  then  be  used  to 
provide  an  estimate  of  fatigue  (Hagg  1991,  Hagg  and  Suurkula 
1991) ,  however,  the  lEMG  trend  at  low  levels  of  contraction 
is  not  reliable  (Seidel,  Beyer,  and  Brauer  1987,  Chapman  and 
Troup  1970) .  Furthermore,  sustained  contractions  to  derive 
a  measure  of  endurance  or  fatigue  would  require  lengthy 
durations  to  induce  measurable  muscle  fatigue.  An 
functional  approach  to  determine  the  fatiguing  effect  of 
repeated  impacts  on  the  paraspinal  muscle  response  may  be  to 
look  at  the  dynamic  response  to  individual  impacts  and  how 
this  response  changes  with  time. 

There  was  an  increase  in  the  magnitude  of  the  lumbar 
muscle  response  to  2  g  z-axis  impacts  after  2  hours  of 
exposure  at  32  impacts  per  minute.  This  result  is 
consistent  with  an  increased  recruitment  of  motor  units  to 
produce  equivalent  force  in  a  fatigued  muscle.  It  is  also 
consistent  with  observations  of  increased  stretch  reflex 
gain  found  to  result  in  a  larger  amplitude  response  to 
muscle  stretch  during  fatiguing  contractions  (Darling  and 
Hayes  1983)  .  An  altered  gain  of  the  stretch  reflex  could  be 
interpreted  as  a  compensatory  mechanism  that  maintains  an 
appropriate  level  of  force  production  despite  muscle 
fatigue.  This  would  therefore  become  deleterious  only  when 
a  maximal  muscle  response  is  required.  Given  the  low  level 
of  muscle  response  typically  observed,  this  is  unlikely  to 
occur  for  shocks  of  3  g  or  less.  An  increased  EMG  amplitude 
may  also  result  from  increased  motor  unit  synchronization; 
however,  this  was  not  visually  evident  in  our  signals. 
Spectral  analysis  can  assist  in  determining  the  contribution 
of  synchronization  to  an  observed  increase  in  EMG  amplitude. 
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Unfortunately,  this  is  not  a  practical  approach  for  an  EMG 
burst  lasting  less  than  100  milliseconds. 

In  addition  to  a  change  in  magnitude,  there  was  also  a 
change  in  the  timing  of  the  EMG  burst  relative  to  the  peak 
impact  acceleration.  However  the  direction  of  this  change 
^  was  not  consistent  in  the  data  of  the  two  subjects  analyzed. 

These  changes  may  be  associated  with  alterations  in  posture, 
or  to  the  preliminary  displacement  of  the  MARS  in  the 
opposite  direction  of  the  intended  impact.  Subjects  may 
^  have  anticipated  the  impact  which  could  alter  the  timing  of 

the  muscle  response.  Under  these  circumstances,  the  timing 
of  the  response  may  not  be  a  very  good  indicator  of  fatigue- 
dependent  changes . 

An  increase  in  background  EMG  was  seen  after  2  hours  in 
subject  3  but  not  in  subject  4.  Although  this  increase  may 
be  a  fatigue  effect,  it  may  also  be  due  to  changes  in 
posture  during  the  test.  Optotrak  data  could  be  analysed  to 
determine  systematic  postural  changes.  A  detailed  postural 
analysis  has  not  yet  been  undertaken  and  would  require  some 
programming  to  process  Optotrak  output  files. 

Qptotrak  displacement  and  acceleration 

A  preliminary  analysis  of  Optotrak  displacement  data 
suggests  that  the  system  would  be  beneficial  in  determining 
spinal  motion  in  response  to  shocks.  A  series  of  IRED's 
attached  to  the  spinous  processes  can  provide  the  coordinate 
data  necessary  for  detailed  biomechanical  analysis  of  the 
spine . 

The  accelerations  calculated  from  Optotrak  displacement 
data  compared  favourably  with  the  corresponding 
accelerometer  data.  The  Optotrak  has  the  advantage  of 
measuring  accelerations  in  three  dimensions  from  a  single 
sensor.  Thus  the  acceleration  data  derived  from  five  IRED's 
would  require  fifteen  accelerometers  attached  to  individual 
spinous  processes.  The  IRED’s  are  relatively  small  in  size, 
and  their  mass  can  be  further  reduced  by  grinding  down  the 
housing  material. 

The  results  obtained  from  Optotrak  analysis  indicate  that 
the  system  is  capable  of  providing  reliable  acceleration 
responses  to  shocks  of  1  to  3  g  magnitude  and  2  to  11  Hz 
^  frequency  in  the  x  and  y  axes.  Although  the  system  can  also 

measure  the  responses  to  similar  shocks  in  the  z-axis  with 
reasonable  accuracy,  the  sampling  rates  used  in  the  study 
did  not  capture  the  high  frequency  components  of 
acceleration  measured  by  the  accelerometers,  and 
superimposed  upon  the  underlying  shock  waveform.  These 
accelerations  only  occurred  in  the  higher  amplitude,  low 
frequency  shocks  and  were  associated  with  the  impact  as  the 
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subject  struck  the  seat,  rather  than  from  the  input  shock 
waveform. 

The  presence  of  missing  data  complicated  the  data 
analysis.  As  drop  outs  coincided  with  the  shocks,  they  may 
have  resulted  from  high  frequency  oscillations,  or  velocity 
affecting  the  resolution  of  the  system.  Further  ^ 

investigation  is  required  to  determine  the  cause  of  drop¬ 
outs,  and  to  determine  the  optimum  parameters  for  IRED's  to 
obtain  higher  sampling  rates,  without  incurring  missing 
data.  Unless  a  higher  sampling  frequency  (approximately  200  \ 

Hz)  can  be  achieved  from  the  Optotrak  system,  without 
missing  data,  it  is  not  considered  practical  to  replace 
accelerometer  data  with  Optotrak  data  in  the  Phase  4 
experiments . 

Biochemistry 

Biochemical  measures  are  subject  to  wide  inter-  and 
intra-individual  variation  at  rest  and  in  response  to 
physiological  stress.  The  small  subject  cohort  may  not 
provide  sufficient  analytical  power  to  identify  strong 
trends  in  the  data.  As  a  result  of  individual  variability, 
a  subject  identified  as  a  "responder"  may  be  masked  by  the 
group  data.  An  obvious  example  of  this  can  be  seen  when  the 
CPK  data  are  plotted  for  each  individual  (Appendix  E,  Figure 
99) .  Subject  1  and  4  have  a  clear  elevation  of  CPK  at  12 -, 

24-  and  36-hours  in  at  least  one  exposure  condition. 

Subject  3  has  a  more  variable  response,  and  subject  2  is 
essentially  a  "non-responder".  In  subject  2,  the  unusually 
high  baseline  value  for  CPK  in  L6  (2g  z  axis  shocks  at 
32 /-min)  is  not  typical.  Taking  into  consideration  the 
other  data  for  this  subject  and  the  lack  of  change  in  CPK  in 
response  to  other  experimental  conditions,  it  is  likely  that 
subject  2  experienced  some  muscular  trauma  or  extreme 
muscular  exertion  one  or  two  days  before  the  study  began. 

Overall,  the  group  mean  data  suggest  that  L2  (rms  only) 
produced  the  greatest  elevation  in  CPK,  which  would  indicate 
the  greatest  involvement  of  muscular  activity. 

Evidence  from  exercise  physiology  literature  suggests 
that  muscles  may  adapt  to  repeated  exercise  sessions,  and  -» 

that  the  release  of  CPK  into  the  blood  stream  will  be 
attenuated  in  repeated  exposures.  Subject  1  clearly  shows 
this  phenomenon.  When  the  data  are  categorized  by  date,  the 
greatest  elevation  in  CPK  is  observed  following  the  first  ^ 

two-hour  exposure  to  experimental  conditions  (L2,  rms 
vibration  with  no  shocks) .  Subsequent  experiments  resulted 
in  a  smaller  elevation  in  CPK  in  this  subject.  In  subject 
4,  CPK  was  elevated  in  only  one  condition.  It  is 
interesting  that  it  was  also  the  first  exposure  to  vibration 
for  this  subject  (L2,  rms  vibration  with  no  shocks) .  In 
subject  2  and  3  a  pattern  of  diminishing  response  for 
subsequent  exposures  could  not  be  clearly  identified. 
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Anderson  et  al . ,  (1977)  reported  a  small  elevation  in 

serum  CPK  immediately  following  nap-of-the-earth  helicopter 
flights.  It  is  difficult  to  compare  this  result  to  the 
present  study,  since  the  greatest  elevation  in  CPK  was 
observed  from  12  to  36  hours  after  the  experimental 
exposure . 

The  degree  of  elevation  of  CPK  in  the  current  study, 
while  clinically  significant,  is  small  compared  to  that 
expected  after  severe  physical  exercise.  Apple  and  Rhodes 
(1988)  reported  a  CPK  of  2,250  U/L  in  one  individual 
following  a  marathon  race,  which  is  roughly  ten-fold  higher 
than  that  observed  in  the  pilot  study.  This  suggests  that 
while  the  amount  of  muscular  involvement  in  response  to  the 
experimental  conditions  is  sufficient  to  produce  some 
leakage  of  CPK  from  skeletal  muscle  to  blood,  overall  it 
represents  a  moderate  physical  stress.  It  should  also  be 
noted  that  the  muscle  mass  under  stress  in  the  current 
experiment  is  much  more  localized  than  during  severe 
physical  exercise.  Therefore  a  lower  CPK  value  in  the 
present  study  could  still  be  associated  with  a  significant 
level  of  stress. 

The  elevation  in  CPK  suggests  that  some  of  the  subjects 
may  have  experienced  some  skeletal  muscle  damage  from  the 
vibration  and  shock  exposures  in  this  study.  In  a  practical 
sense,  a  "responder"  in  a  variable  which  suggests  damage  to 
skeletal  muscle  may  be  at  greater  risk  of  injury  if  other 
physical  activity  is  attempted  concurrently  or  immediately 
following  exposure  to  repeated  shocks  and  impacts. 

While  other  biochemical  results  did  not  produce  any 
indisputable  evidence  indicative  of  fatigue  or  tissue  damage 
in  the  present  experiments,  the  data  are  consistent  with 
other  measures  in  the  study.  Neither  EMG  nor  ECG  data 
reveal  obvious  fatigue,  which  is  supported  by  the  lack  of  a 
clear  biochemical  marker  (elevated  lactate,  ammonia,  or 
potassium) .  A  small  reduction  in  blood  glucose  in  all 
experiments  including  the  control  condition  underline  the 
need  for  regular  nutrient  intake  to  maintain  blood  glucose 
concentration.  However  the  duration  of  the  present 
experiment  was  not  sufficiently  long  to  result  in 
hypoglycemia  in  any  condition. 

There  were  other  problems  which  may  affect  the 
interpretation  of  some  of  the  biochemical  data  including: 
poor  quality  control  in  some  biochemical  analysis;  high 
resting  values  in  some  biochemical  variables;  general 
stress  associated  with  blood  sampling;  a  small  subject 
number  with  a  high  variability  in  subject  characteristics; 
and  two  subjects  who  appeared  to  have  had  a  sub- symptomatic 
illness  based  on  their  high  white  blood  cell  count.  It  is 
also  quite  likely  that  the  exposure  was  not  long  enough  or 
severe  enough  to  adversely  affect  the  subjects. 
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There  is  a  definite  need  for  biochemical  measures  to  be 
included  in  the  development  of  a  health  hazard  assessment. 

To  reduce  the  problems  in  the  interpretation  of  biochemical 
data,  future  experiment  should: 

•  have  longer  duration,  higher  exposure  experiments 
to  increase  the  likelihood  of  an  effect  on  the 
subject. 

•  include  kidney  function  tests  including  a  measure 
of  glomerular  filtration  rate  before  and  after  the 
experiment.  This  would  involve  a  24-Hour  urine 
collection  for  measurement  of  creatinine 
clearance . 

•  have  a  larger  number  of  subjects  in  the  long-term 
experiments . 

•  take  blood  samples  to  screen  the  subjects  for  any 
sub-symptomatic  illness  which  confounds  the  data. 

In  addition,  biochemical  data  obtained  following  "real" 
exposure  to  vehicle  stress,  for  example  in  prolonged  field 
studies,  would  provide  further  insight  to  the  physiological 
stress  associated  with  whole  body  vibration  and  repeated 
shock. 
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Conclusions 


*4 


1.  Responses  measured  in  spinal  acceleration,  internal 
pressure,  chest  and  abdominal  displacement,  and  EMG 
showed  similar  patterns  of  frequency  response. 

2 .  Responses  measured  for  shock  frequencies  of  2  to  11  Hz 
do  not  agree  with  transmission  (weighting)  curves  in 
current  standards  (ISO  2631,  1982;  BSI  6841,  1987; 
ASCC,  1982) . 

3 .  There  is  evidence  of  non-linearity  in  response  to 
shocks  as  reflected  in: 

•  Changes  in  frequency  of ‘peak  (transmission) 
response  in  the  x-and  z-axis  with  different  shock 
magnitudes . 

•  Changes  in  transmission  ratio  with  different 
magnitudes  of  shock. 

•  Shape  of  the  spinal  acceleration  response  to 
individual  shocks  in  the  negative  z-axis. 

3 .  The  magnitude  of  response  to  negative  z-axis  shocks  was 
higher  than  the  response  to  x-and  y-axis  shocks . 

4 .  The  dominant  spinal  acceleration  and  internal  pressure 
responses  to  negative  z-axis  shocks  are  associated  with 
the  subject  hitting  the  seat.  This  response  contained 
very  high  frequency  components  (>20  Hz) . 

5 .  The  pilot  experiments  did  not  show  conclusive  evidence 
of  fatigue  induced  by  2  hour  exposures  to  shock  and 
vibration  in  either  biochemical  indicators,  EMG 
response  or  ECG  parameters. 

6 .  There  was  biochemical  evidence  of  muscle  damage  in  some 
subjects  following  2 -hour  exposures  to  shocks  and 
vibration. 

7 .  The  magnitude  of  muscle  response  to  shocks  is  typically 
less  than  10  percent  of  maximal  voluntary  contraction. 

8.  The  pattern  of  muscle  response  to  shocks  involves  two 
phases:  stabilization  of  the  upper  torso  and  re¬ 
establishment  of  a  neutral  posture. 

9 .  Performance  measures  induced  changes  in  some  ECG 
parameters . 
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Recommendations 


Standards  developed  for  exposure  to  vibration  and 
repeated  shocks  should  account  for: 

•  non-linearity  of  response 

•  differing  responses  to  x,  y,  and  z-axis  inputs 

•  differing  responses  to  positive  and  negative 

directions  of  shocks  in  the  x-and  z-axes 

Further  investigations  of  individual  shock  responses 
are  required,  including: 

•  shocks  in  the  negative  x-axis  and  positive  z-axis 
directions 

•  shocks  at  low  frequencies  (for  example,  1  to  4  Hz) 

•  higher  frequencies  of  shocks  (for  example,  >20  Hz) 

•  larger  magnitudes  of  shocks 

Further  investigation  is  required  of  cumulative 
exposures  that  are  of  longer  duration  and  increased 
severity  to  more  accurately  simulate  a  typical  military 
mission. 

More  frequent  recovery  measures  should  be  taken  over  a 
longer  recovery  period  to  observe  possible  fatigue. 

The  fourth  phase  of  the  project  should  include  the 
following  measures:  acceleration  at  the  spine; 
displacement  of  the  spine  (measured  by  Optotrak) ; 
internal  pressure;  EMG  (including  more  muscles  and 
sustained  contractions  at  levels  similar  to  those 
induced  by  shocks) ;  ECG;  biochemical  markers  (including 
hydroxyproline,  lactate,  K  ,  CPK  and  glomerular 
filtration  rates);  and  performance  measures. 
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Appendix  B 

Subjective  Assessment  Technicmes 


The  following  questions  will  be  asked  verbally  of  the 
subject  during  the  short-time  experiments  and  long-time 
experiments  and  recorded  by  an  experimenter.  In  the  short- 
tenn  experiments  (5.5  minutes)  the  questions , will  be  asked 
during  the  30  second  standard  exposure,  at  the  beginning  of 
the  exposure  and  just  prior  to  completion.  In  the  long-term 
experiments  (1-hour  and  2  hour) ,  questions  will  be  asked 
during  the  standard  exposure,  at  the  beginning  of  the 
exposures  and  during  the  20  minute  periods  when  subjects  are 
not  performing  the  multi-task,  as  well  as  just  prior  to 
completion  of  the  experiment. 

All  responses  will  be  made  to  a  seven  point  scale.  The 
scale  is  shown: 


not  at  all  extremely 

Questions : 


1.  How  uncomfortable  is  the  vibration? 

2.  How  uncomfortable  are  the  shocks? 

3.  How  tired  do  you  feel? 

4.  How  has  the  vibration  and  shocks  affected  your 
performance  of  tasks? 


the  method  used  in  each  biochemical  analysis,  the  normal  range  expected 
e,  and  the  coefficient  of  variation  (CV)  of  each  method. 
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Transmission 

Frequency _ Mean _ Delay  Mean 

4  0.737  0.159 

5  0.482  0.106 
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Chest  respitrace  response  ratio  and  response  delay  for  -2  g  z-axis  shocks 
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Table  65 


EC6  Responses 

to  Control 

Condi  t: 

Lon 

(2 

hour  exposure) 

Subject 

Exp.  Phase 

Heart  rate 
(bpm) 

Mean  ±  SD 

T-wave 

amplitude 

(volts) 

Mean  ±  SD 

R-R 

interval 

(sec) 

Mean  ±  SD 

Spectral 
Variance 
(0.06-0.1  Hz) 

Spectral 
Variance 
(0.12- 0.4  Hz) 

1 

Rest 

95.2 

± 

1.41 

0.12 

± 

0.02 

0.63 

± 

0.03 

10.48 

203.84 

Epoch  1 

96.22 

± 

1.89 

0.12 

± 

0.03 

0.62 

± 

0.03 

5.95 

204.58 

2 

95.33 

± 

1.63 

0.62 

± 

0.02 

2.85 

170.97 

*3 

96.22 

± 

2.22 

0.62 

± 

0.03 

6.14 

206.26 

6 

93.33 

± 

1.26 

0.14 

± 

0.04 

0.64 

± 

0.05 

22.04 

787.94 

10 

88.89 

± 

1.5 

0.67 

± 

0.05 

10.87 

476.44 

*11 

90.67 

± 

1.71 

0.66 

,  ± 

0.03 

7.34 

415.74 

12 

95.78 

± 

1.71 

0.11 

± 

0.04 

0.63 

± 

0.06 

31.56 

635.04 

(1 50  sec) 

Recovery 

90.4 

± 

0.71 

0.19 

± 

0.02 

0.66 

± 

0.05 

62.08 

777.42 

2 

Rest 

72.40 

± 

1.41 

0.30 

± 

0.04 

0.83 

± 

0.03 

4.51 

222.81 

Epoch  1 

73.78 

± 

2.50 

0.28 

± 

0.06 

0.81 

± 

0.04 

16.6 

246.81 

2 

72.44 

± 

2.65 

0.83 

± 

0.06 

24.14 

490.86 

‘3 

71.11 

± 

1.50 

0.36 

± 

0.03 

0.84 

± 

0.02 

4.06 

185.59 

6 

68.44 

± 

0.96 

0.88 

± 

0.03 

13.74 

350.05 

10 

64.00 

± 

2.50 

0.94 

± 

0.07 

25.8 

867.03 

*11 

67.78 

± 

2.08 

0.88 

± 

0.03 

5.75 

286.87 

12 

70.44 

± 

2.52 

0.34 

± 

0.1 

0.86 

± 

0.08 

13.73 

885.43 

(600  sec) 

Recovery 

64.30 

± 

1.62 

0.40 

0.04 

0.93 

± 

0.05 

24.88 

586.49 

3 

Rest 

59.20 

± 

1.41 

0.32 

± 

0.03 

1.01 

± 

0.03 

13.31 

431.77 

Epoch  1 

62.89 

± 

2.75 

0.29 

± 

0.05 

0.95 

± 

0.07 

27.86 

897.28 

2 

60.67 

± 

1.71 

0.98 

± 

0.07 

26.82 

959.36 

*3 

62.89 

± 

3.30 

0.95 

± 

0.06 

12.18 

512.48 

6 

59.33 

± 

2.87 

0.33 

± 

0.05 

1.01 

± 

0.08 

41.68 

1728.90 

10 

63.11 

± 

2.99 

0.95 

± 

0.07 

23.61 

725.04 

*11 

71.78 

± 

2.22 

0.83 

± 

0.05 

19.23 

640.04 

12 

69.78 

± 

1.91 

0.27 

± 

0.07 

0.86 

± 

0.10 

53.64 

2179.00 

4 

Rest 

78.00 

± 

1.41 

0.23 

± 

0.03 

0.76 

± 

0.05 

35.86 

365.03 

Epoch  1 

76.67 

± 

2.06 

0.24 

± 

0.02 

0.78 

± 

0.05 

17.73 

493.39 

2 

75.56 

± 

2.22 

0.79 

± 

0.06 

27.18 

963.4 

*3 

72.44 

± 

2.83 

0.83 

± 

0.06 

19.11 

1953.6 

6 

71.78 

± 

0.82 

0.24 

± 

0.04 

0.84 

± 

0.08 

39.78 

1787.4 

10 

70.44 

± 

3.77 

0.85 

± 

0.09 

35.52 

2313.4 

*11 

67.78 

± 

1.26 

0.88 

± 

0.08 

32.51 

2407.5 

12 

71.33 

± 

3.11 

0.25 

± 

0.04 

0.84 

± 

0.08 

16.43 

2366.2 

(150  sec) 

Recovery 

84.00 

± 

0.00 

0.26 

± 

0.03 

0.94 

± 

0.11 

186.58 

3222.4 

Synchronised  Work  Battery  performed  during  this  epoch 


Table  66 

ECG  Response  to  3g  z-ajcis  1  8hock/2.5  min  Condition 
_ _ (2  hour  exposure) _ _ 


Subject 

Exp.  Phase 

Heart  rate 
(bpm) 

Mean  ±  SD 

T-wave 

amplitude 

(volts) 

Mean  ±  SD 

R-R 

interval 

(msec) 

Mean  ±  SD 

Spectral 
Variance 
(0.06- 
0.1  Hz) 

Spectral 
Variance 
(0. 12- 
0.4  Hz) 

1 

Rest 

79.60 

± 

0.71 

± 

0.02 

0.77 

± 

0.06 

20.54 

1475.4 

Epoch  1 

83.33 

± 

2.22 

0.18 

± 

0.03 

0.74 

± 

0.05 

6.44 

946.22 

2 

78.67 

± 

2.16 

0.77 

± 

0.04 

12.34 

664.56 

*3 

83.11 

± 

0.96 

0.73 

± 

0.04 

7.13 

339.10 

6 

74.89 

± 

2.75 

0.22 

± 

0.03 

0.81 

± 

0.07 

44.23 

1018.6 

10 

79.33 

± 

3.69 

0.76 

± 

0.10 

93.52 

1422.00 

•11 

81.33 

± 

1.73 

0.74 

± 

0.04 

10.97 

505.4 

12 

78.89 

± 

6.19 

± 

0.05 

0.78 

± 

0.08 

41.34 

1071.2 

(150  sec) 

Recovery 

77.60 

± 

2.12 

0.19 

± 

0.02 

0.77 

± 

0.59 

64.13 

1065.8 

2 

Rest 

81.20 

± 

2.12 

0.25 

± 

0.04 

0.74 

± 

0.02 

4.09 

42.87 

Epoch  1 

77.78 

± 

1.5 

0.29 

± 

0.04 

0.77 

± 

0.03 

4.74 

81.00 

2 

79.33 

± 

1.5 

0.76 

± 

0.04 

17.99 

269.97 

•3 

^7.33 

± 

0.5 

0.77 

± 

0.03 

9.79 

87.89 

6 

74.00 

± 

2.16 

0.32 

± 

0.03 

0.81 

± 

0.03 

9.45 

199.33 

10 

73.56 

± 

0.5 

0.81 

± 

0.04 

11.82 

233.35 

•11 

75.33 

± 

1.26 

0.79 

± 

0.02 

3.98 

83.29 

12 

73.78 

± 

1.89 

0.31 

± 

0.04 

0.81 

± 

0.04 

10.07 

209.84 

(600  sec) 

Recovery 

74.30 

± 

2.57 

0.33 

± 

0.06 

0.81 

± 

0.03 

12.64 

147.18 

3 

Rest 

56.00 

± 

0.00 

0.28 

± 

0.03 

1.07 

± 

0.04 

36.25 

498.71 

Epoch  1 

58.89 

± 

1.71 

0.27 

± 

0.04 

1.02 

± 

0.07 

25.5 

961.03 

2 

59.78 

± 

2.89 

1.00 

± 

0.09 

43.26 

1468.5 

•3 

64.22 

± 

3.4 

0.93 

± 

0.07 

24.60 

997.06 

6 

58.89 

± 

3.27 

0.29 

± 

0.05 

1.02 

± 

0.09 

50.90 

1089.50 

10 

58.67 

± 

3.87 

1.02 

± 

0.1 

56.45 

1353.4 

•11 

65.11 

± 

1.29 

0.92 

± 

0.07 

19.49 

1638.5 

12 

59.11 

± 

2.5 

0.32 

± 

0.05 

1.02 

± 

0.07 

38.47 

989.19 

4 

Rest 

74.89 

± 

0.96 

0.19 

± 

0.02 

0.8 

± 

0.08 

34.28 

1332.4 

Epoch  1 

75.1 

± 

1.29 

0.17 

± 

0.03 

0.8 

± 

0.07 

22.50 

2212.5 

2 

75.56 

± 

2.16 

0.79 

± 

0.06 

17.84 

1985.9 

•3 

72.00 

± 

1.15 

0.83 

± 

0.08 

29.12 

2597.3 

6 

73.11 

± 

2.75 

0.21 

± 

0.04 

0.82 

± 

0.08 

26.21 

2799.2 

10 

74.00 

± 

2.22 

0.81 

± 

0.09 

59.60 

1948.2 

•11 

77.78 

± 

2.5 

0.77 

± 

0.06 

9.26 

1162.8 

12 

72.22 

± 

2.94 

0.22 

± 

0.04 

0.83 

± 

0.09 

27.85 

2367.8 

(150  sec) 

Recovery 

± 

1.41 

0.22 

± 

0.03 

0.82 

± 

0.07 

31.66 

962.48 

•Synchronised  Work  Battery  performed  during  this  epoch 
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Table  67 

EC6  Response  to  2g  z-axis  32 /min  Condition 
_ _ (2  hour  exposure) _ p__ 


Subject 

Exp.  Phase 

Heart  rate 
(bpm) 

Mean  ±  SD 

T-wave 

amplitude 

(volts) 

Mean  ±  SD 

R-R 

interval 

(msec) 

Mean  ±  SD 

Spectral 
Variance 
(0.06- 
0.1  Hz) 

Spectral 
Variance 
(0.12- 
0.4  Hz) 

1 

Rest 

88.00 

± 

2.83 

0.12 

±  0.02 

0.76 

± 

0.07 

72.44 

1549.60 

Epoch  1 

86.22 

± 

3.40 

0.14 

±  0.03 

0.75 

± 

0.06 

15.37 

1135.60 

2 

75.33 

± 

1.41 

0.79 

± 

0.04 

18.40 

460.58 

*3 

83.11 

± 

0.96 

0.73 

± 

0.03 

3.85 

515.47 

6 

71.33 

± 

2.50 

0.2± 

0.04 

0.84 

± 

0.05 

14.35 

455.12 

10 

73.11 

± 

4.57 

0.82 

± 

0.06 

14.14 

468.68 

*11 

80.67 

± 

1.29 

0.74 

± 

0.03 

7.37 

363.17 

12 

81.11 

± 

4.55 

0.15 

±  0.06 

0.76 

± 

0.05 

16.87 

673.47 

(600  sec) 

Recovery 

82.20 

± 

2.95 

0.12 

±  0.05 

0.75 

± 

0.06 

38.20 

1058.70 

3 

Rest 

54.80 

± 

0.71 

0.34 

±  0.03 

1.09 

0.04 

38.59 

1884.80 

Epoch  1 

64.89 

± 

2.22 

0.39 

±  0.14 

0.96 

± 

0.08 

30.18 

1172.40 

2 

72.44 

± 

8.38 

0.90 

± 

0.10 

35.04 

2875.30 

•3 

65.56 

± 

2.45 

0.92 

± 

0.05 

12.29 

475.59 

6 

72.00 

± 

5.16 

0.38 

±  0.17 

0.91 

0.09 

25.49 

2411.90 

10 

64.67 

± 

2.89 

0.94 

± 

0.06 

29.66 

543.47 

*11 

75.56 

± 

2.65 

0.81 

± 

0.06 

18.44 

344.70 

12 

65.78 

± 

1.71 

0.39 

±  0.09 

0.91 

0.07 

23.3 

578.14 

(150  sec) 

Recovery 

58.40 

± 

0.71 

0.36 

±  0.02 

1.03 

± 

0.03 

8.48 

316.94 

4 

Rest 

78.40 

± 

4.95 

0.19 

±  0.06 

0.76 

± 

0.07 

23.57 

776.90 

Epoch  1 

76.44 

± 

1.0 

0.22 

±  0.04 

0.78 

± 

0.07 

13.97 

611.64 

2 

78.89 

± 

2.06 

0.76 

± 

0.04 

9.95 

619.36 

*3 

79.33 

± 

2.36 

15.81 

1198.7 

6 

78.00 

± 

0.71 

0.22 

±  0.05 

0.77 

± 

0.06 

20.14 

1322.9 

10 

76.67 

± 

2.06 

0.78 

± 

0.06 

10.80 

1909.8 

*11 

83.56 

± 

1.29 

0.75 

± 

0.06 

34.50 

1287.2 

12 

34.89 

± 

6.95 

0.22 

±  0.07 

0.78 

± 

0.07 

26.81 

1250.80 

•Synchronised  Work  Battery  performed  during  this  epoch 

«  ’ 
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Table  68 

Background  and  impact  response  XEM6  for  the  first 
emd  last  sampling  intervals  of  a  2  hour  exposure 
to  2  g  z-eucis  impacts  at  32/min. 


Background  lEMG  (mV) 

Impact  Response  lEMG  (mV) 

First  Interval 

Last  Interval 

First  Interval 

Last  Interval 

Subject 

Mean 

St.  Dev. 

Mean 

St.  Dev. 

Mean 

St.  Dev. 

Mean 

St.  Dev. 

3 

9 

4 

61 

20 

59 

29 

75 

31 

4 

7 

i 

2 

9 

4 

24 

7 

35 

11 
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Table  69 

Meem,  nuucimum,  and  standard  deviation  of  the  percent 
HVC  muscle  response  to  3  g  inqpact  accelerations. 


Frequency 

(Hz) 

x-axis  muscle  response 
(%MVC) 

y-axis  muscle  response 
(%MVC) 

z-axis  muscle  response 
(%MVC) 

max 

mean 

St.  dev. 

max 

mean 

st.  dev. 

max 

mean 

st.  dev. 

11 

13.2 

3.6 

4.7 

14.1 

1.7 

3.7 

7.6 

2.6 

2.0 

8 

32.5 

8.4 

9.0 

19.7 

2.9 

5.5 

31.9 

5.5 

7.9 

6 

18.7 

6.6 

4.9 

12.6 

2.4 

3.4 

34.8 

6.4 

8.6 

5 

15.7 

6.8 

4.6 

24.1 

3.0 

5.9 

34.6 

7.8 

8.4 

4 

13.2 

6.3 

3.9 

23.1 

3.6 

6.1 

65.4 

14.3 

18.0 

Table  70 

The  freguency  response  (n=4)  of  urinazy  biochemical 
variables  measured  pre-  and  post-  exposure  to  at  two-hour 
exposure  in  four  experimental  conditions. 

(Exp.  L1  -  control; 

Exp.  L2  -  RMS; 

Exp.  L6  -  2  g  z-axis  at  32/min: 

Exp.  L7  - 1  g  X,  y,  z-axis  at  32/min) 


A.  Urinary  Appearance 


Pre 

Post 

Frequency 

Percent 

Frequency 

Percent 

L1  Clear 

3 

75 

2 

50 

Slightly  Cloudy 

1 

25 

0 

0 

Cloudy 

0 

0 

1 

25 

Hazy 

0 

0 

0 

0 

Turbid 

0 

0 

0 

0 

Missing 

0 

0 

1 

25 

L2  Clear 

1 

25 

3 

75 

Slightly  Cloudy 

0 

0 

0 

0 

Cloudy 

0 

0 

0 

0 

Hazy 

3 

75 

1 

25 

Turbid 

0 

0 

0 

0 

Missing 

0 

0 

0 

0 

L6  Clear 

2 

50 

3 

75 

Slightly  Cloudy 

0 

0 

0 

0 

Cloudy 

1 

25 

0 

0 

Hazy 

0 

0 

0 

0 

Turbid 

1 

25 

1 

25 

Missing 

0 

0 

0 

0 

L7  Clear 

2 

50 

2 

50 

Slightly  Cloudy 

0 

0 

0 

0 

Cloudy 

2 

50 

0 

0 

Hazy 

0 

0 

0 

0 

Turbid 

0 

0 

0 

0 

Missing 

0 

0 

2 

50 
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Table  70  (cont'd) 


B.  Urinary  Bacteria 


Pre 

Post 

Result 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

4 

100 

4 

100 

1+ 

0 

0 

0 

0 

L2  Negative 

3 

75 

3 

75 

1+ 

1 

25 

1 

25 

L6  Negative 

3 

75 

3 

75 

1+ 

1 

25 

1 

25 

L7  Negative 

4 

100 

4 

100 

1+ 

0 

0 

0 

0 

C.  Urinary  Bilirubin 

Pre 

Post 

Result 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

4 

100 

3 

75 

Missing 

0 

0 

1 

25 

12  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L6  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L7  Negative 

4 

100 

2 

50 

Missing 

0 

0 

2 

50 

D.  Urinary  Blood 

Pre 

Post 

Result 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

4 

100 

3 

75 

Missing 

0 

0 

1 

25 

L2  Negative 

3 

75 

4 

100 

Missing 

1 

25 

0 

0 

L6  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L7  Negative 

4 

100 

2 

50 

Missing 

0 

0 

2 

50 
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Table  70  (cont'd) 


E.  Urinary  Casts 


Pre 

Post 

Result 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

3 

75 

4 

100 

1+ 

1 

25 

0 

0 

2+ 

0 

0 

0 

0 

3+ 

0 

0 

0 

0 

L2  Negative 

4 

100 

4 

100 

1  + 

0 

0 

0 

0 

2+ 

0 

0 

0 

0 

3+ 

0 

0 

0 

0 

L6  Negative 

3 

75 

4 

100 

1+ 

0 

0 

0 

0 

2+ 

0 

0 

0 

0 

3+ 

1 

25 

0 

0 

L7  Negative 

4 

100 

4 

100 

1+ 

0 

0 

0 

0 

2+ 

0 

0 

0 

0 

3+ 

0 

0 

0 

0 

F. 

Urinary  Crystals  (Calcium  Oxalate) 

Pre 

Post 

Result 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

3 

75 

3 

75 

Few 

1 

25 

1 

25 

L2  Negative 

2 

50 

2 

50 

Few 

2 

50 

2 

50 

L6  Negative 

2 

50 

4 

100 

Few 

2 

50 

0 

0 

L7  Negative 

3 

75 

4 

100 

Few 

1 

25 

0 

0 

Table  70  (cont'd) 


G.  Urinary  Color 


Pre 

Post 

Frequency 

Percent 

Frequency 

Percent 

L1  Yellow 

4 

100 

3 

75 

Straw 

0 

0 

0 

0 

Pinkish-Yellow 

0 

0 

0 

0 

Missing 

0 

0 

1 

25 

L2  Yellow 

3 

75 

4 

100 

Straw 

1 

25 

0 

0 

Pinkish-Yellow 

0 

0 

0 

0 

L6  Yellow 

2 

50 

1 

25 

Straw 

2 

50 

3 

75 

Pinkish-Yellow 

0 

0 

0 

0 

L7  Yellow 

3 

75 

2 

50 

Straw 

0 

0 

0 

0 

Pinkish-Yellow 

1 

25 

0 

0 

Missing 

0 

0 

2 

50 

H.  Urinary  Epithelial  Cells 

Pre 

Post 

Result 

Frequency 

Percent 

Frequency  i 

Percent 

LI  Negative 

0 

0 

2 

50 

^+ 

3 

75 

2 

50 

2+ 

1 

25 

0 

0 

3+ 

0 

0 

0 

0 

L2  Negative 

1 

25 

2 

50 

1  + 

2 

50 

2 

50 

2+ 

1 

25 

0 

0 

3+ 

0 

0 

0 

0 

L6  Negative 

0 

0 

1 

25 

^+ 

2 

50 

2 

50 

2+ 

1 

25 

0 

0 

3+ 

1 

25 

1 

25 

L7  Negative 

1 

25 

2 

50 

1-1- 

2 

50 

2 

50 

2+ 

1 

25 

0 

0 

3+ 

0 

0 

0 

0 
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Table  70  (cont'd) 


1.  Urinary  Glucose 

Pre 

Post 

Result 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

4 

100 

3 

75 

Missing 

0 

0 

1 

25 

12  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L6  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L7  Negative 

4 

100 

2 

50 

Missing 

0 

0 

2 

50 

J.  Urinary  Ketones 

Pre 

Post 

Result 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

4 

100 

3 

75 

Missing 

0 

0 

1 

25 

L2  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L6  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L7  Negative 

4 

100 

2 

50 

Missing 

0 

0 

2 

50 

K.  Urinary  Leukocytes 


Result 

Pre 

Frequency 

Percent 

Post 

Frequency 

Percent 

LI  Negative 

4 

100 

3 

75 

Missing 

0 

0 

1 

25 

L2  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L6  Negative 

4 

100 

3 

75 

Missing 

0 

0 

1 

25 

L7  Negative 

4 

100 

2 

50 

Missing 

0 

0 

2 

50 

4 
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Table  70  (cont'd) 


L  Urinary  Mucous 


Pre 

Post 

Result 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

1 

25 

2 

50 

Trace 

0 

0 

1 

25 

1+ 

1 

25 

1 

25 

2+ 

2 

50 

0 

0 

3+, 

0 

0 

0 

0 

4+ 

0 

0 

0 

0 

L2  Negative 

1 

25 

2 

50 

Trace 

2 

50 

0 

0 

1+ 

1 

25 

1 

25 

2+ 

0 

0 

0 

0 

3+ 

0 

0 

0 

0 

4+ 

0 

0 

1 

25 

L6  Negative 

2 

50 

3 

75 

Trace 

0 

0 

0 

0 

1  + 

0 

0 

1 

25 

2+ 

1 

25 

0 

0 

3+ 

4+ 

1 

25 

0 

0 

L7  Negative 

4 

100 

2 

50 

Trace 

0 

0 

1 

25 

1  + 

0 

0 

1 

25 

2+ 

0 

0 

0 

0 

3+ 

0 

0 

0 

0 

4+ 

0 

0 

0 

0 

M.  Urinary  Myoglobin 


Result 

Pre 

Post 

Frequency 

Percent 

Frequency 

Percent 

LI 

Negative 

4 

100 

4 

100 

L2 

Negative 

3 

75 

4 

100 

Missing 

1 

25 

L6 

Negative 

4 

100 

4 

100 

L7 

Negative 

4 

100 

4 

100 
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Table  70  (cont'd) 


N.  Urinary  Nitrites 


Result 

Pre 

Post 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

4 

100 

3 

75 

Missing 

0 

0 

1 

25 

L2  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L6  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L7  Negative 

4 

100 

2 

50 

Missing 

0 

0 

2 

50 

O.  Urinary  Protein 


Result 

Pre 

Frequency 

Percent 

Post 

Frequency 

Percent 

LI  Negative 

4 

100 

3 

75 

Missing 

0 

0 

1 

25 

L2  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L6  Negative 

4 

100 

4 

100 

Missing 

0 

0 

0 

0 

L7  Negative 

4 

100 

2 

50 

Missing 

0 

0 

2 

50 

Table  70  (cont'd) 


P.  Urinary  Red  Blood  Cells 


Result 

Pre 

Post 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

2 

50 

2 

50 

1+ 

1 

25 

1 

25 

2+ 

0 

0 

0 

0 

3+ 

0 

0 

1 

25 

4+ 

1 

25 

0 

0 

L2  Negative 

4 

100 

2 

50 

1+ 

0 

0 

2 

50 

2+ 

0 

0 

0 

0 

3+ 

0 

0 

0 

0 

4+ 

0 

0 

0 

0 

L6  Negative 

3 

75 

2 

50 

1+ 

1 

25 

1 

25 

2+ 

0 

0 

1 

25 

3+ 

0 

0 

0 

0 

4+ 

0 

0 

0 

0 

L7  Negative 

2 

50 

3 

75 

1+ 

1 

25 

0 

0 

2+ 

0 

0 

0 

0 

3+ 

0 

0 

0 

0 

4+ 

1 

25 

1 

25 

Q.  Urinary  White  Blood  Ceils 

Pre 

Post 

Result 

Frequency 

Percent 

Frequency 

Percent 

LI  Negative 

0 

0 

1 

25 

1+ 

3 

75 

3 

75 

2+ 

1 

25 

0 

0 

L2  Negative 

1 

25 

1 

25 

1+ 

2 

50 

2 

50 

2+ 

1 

25 

1 

25 

L6  Negative 

3 

0 

2 

50 

1+ 

1 

25 

2 

50 

2+ 

0 

0 

0 

0 

L7  Negative 

2 

50 

3 

75 

1+ 

1 

25 

0 

0 

2+ 

1 

25 

1 

25 
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accel-LX  (m/sGc''2) 


123.25 
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123.38  123.5 


123.63  123.75  123.88  124 

Time  (Seconds) 


Acceleration  at  the  spine  (L4  z)  for 
axis  shock.  Dotted  line:  unfiltered  d 
band  pass  filtered  at  0.5  Hz  to  60  Hz 
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60.0 

124 


Igur 


-25.0 


123.24  123.3  123.36  123.42  123.48  123.54  123.6  123.66 

Time  (Seconds) 


Figure  8.  Acceleration  at  the  spine  {L4  z)  for  a  -3  g,  4  Hz  z- 
axis  shock.  Dotted  line:  unfiltered  data;  broken 
line:  filtered  data;  full  line:  filtered  data 
multiplied  by  inverse  transfer  function. 


el-SX  (m/sec''2  ) 


accel-SX  (m/sec''2) 


accel-SY  (m/sec''2)  „  accel-SY  {m/sec''2) 


Igure  16. 


Acceleration  measured  at  the  seat,  lumbar  and 
thoracic  spine  for  a  3  g,  4  Hz  y-axis  shock.  Dotted 
line:  seat  Sy;  broken  line:  lumbar  L3  y;  full  line: 
thoracic  T2  y. 


Figure  17.  Acceleration  measured  at  the  seat,  lumbar  and 

thoracic  spine  for  a  3  g,  11  Hz  y-axis  shock.  Dotted 
line:  seat  Sy;  broken  line:  lumbar  L3  y;  full  line: 
thoracic  T2  y. 
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Figure  26 


z-axis  acceleration  measured  at  the  seat  and  x-axis 
acceleration  measured  at  the  lumbar  and  thoracic  spine 
for  a  3  g,  4  Hz  z-axis  shock.  Dotted  line:  seat  Sz; 
broken  line:  lumbar  L2  x;  full  line:  thoracic  Tl  x. 
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Figure  29. 


Internal  pressure  measured  for  a  3  g,  _4  Hz  x-^is 
shock.  Dotted  line:  seat  Sx  acceleration;  full  line 
internal  pressure. 
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Figure  30:  (a,b,c)  Internal  pressure  response  to  seat  x 

acceleration  for  3,  2,  &  1  g  shocks 
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Figure  31.  Internal  pressure  measured  for  a  3  g,  4  Hz  y-axis 

shock.  Dotted  line:  seat  Sy  acceleration;  full  line 
internal  pressure. 
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Figure  35.  Abdominal  respi trace  displacement  for  a  3  g,  4  Hz  x 
axis  shock.  Dotted  line:  seat  Sx  acceleration;  full 
line:  abdominal  displacement. 
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Figure  36:  (a,b,c)  Abdominal  respi trace  response  to  seat  x 

acceleration  for  3 ,  2 ,  &  1  g  shocks 


E-29 


acceleration  for  3,  2,  &  1  g  shocks 


E-32 


acceleration  for  3,  2,  &  1  g  shocks 


E-35 


(uidq)  nau  V«*H 


is  s  e  s  8  s  s 

(uriq)  nm  umH 


ggsesssss 

ttHlO  too)  WWWA  M=»0s 


P 

O) 

•H 

b 


E-38 


E-40 


40.125  40.25  40.375  40.5  40.625  40.75  40.875  41.0  41.125 

Time  (  Seconds ) 


Figure. 49  A  typical  response  of  lumbar  muscle  (volts)  to  a 
negative  1  g  x-axis  impact  acceleration  at  a 
frequency  of  6  Hz . 
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Figure  50  .^The-  typical  response  of  both  right  and  left  lumbar 
muscle  (volts)  to  a  positive  1  g  y-axis  impact 
acceleration  at  a  frequency  of  6  Hz. 
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Figure  51  The  typical  response  of  both  right  and  left  lumbar 


muscle  (volts)  to  a  negative  1  g  y-axis  impact 
acceleration  at  a  frequency  of  6  Hz . 
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Figure  52  A  typical  response  of  lumbar  muscle  (volts)  to  a 
positive  1  g  z-axis  impact  acceleration  at  a 
frequency  of  6  Hz . 
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Figure  54  The  mean  response  (n=20)  of  lumbar  muscle  to  negative 
3  g  impact  accelerations  at  frequencies  of  4  to  11  Hz 
in  X,  y,  and  z  axes.  Lximbar  muscle  response  is 
normalized  to  percent  of  maximal  voluntary 
contraction  (%MVC)  for  each  subject. 
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Figure  55  The  internal  pressure  (mm  Hg)  and  rectus  abdominus  . 

(RA)  muscle  activity  during  a  maximal  voluntary 
contraction  ("bearing  down")  in  which  the  subject  was 
instructed  to  produce  as  much  internal  pressure  as 
possible . 
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Figure  56  The  internal  pressure  (itimHg)  and  activity  of  the  _ 

rectus  abdominus  (RA)  and  lumbar  (LL)  muscles  during 
first  coughing  (2  to  6  seconds)  and  then  laughing  (7 
to  14  Seconds) . 


Figure  57:  The  internal  pressure  (mm  Hg)  and  activity  of  the 
rectus  abdominus  (RA)  and  lumbar  (L3)  muscles  in 
response  to  a  positive  z-axis  3  g  impact 
acceleration,  shown  as  measured  over  the  fourth 
lumbar  spinous  process  (L4) . 
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Figure  58.  Seat  displacement  (full  line)  and  calculated 

acceleration  (dotted  line)  measured  by  Optotrak 
for  3  g,  6  Hz  x-axis  shock. 
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Figure  59.  Seat  displacement  (full  line)  and  acceleration 

(dotted  line)  measured  by  Optotrak  for  3  g,  5  Hz  y- 
axis  shock. 
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Figure  60.  Seat  displacement  (full  line)  and  acceleration 

(dotted  line)  measured  by  Optotrak  for  -3  g,'4  Hz  z- 
axis  shock. 
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Figure  62.  3  g,  5  Hz  y-axis  shock  at  the  seat  measured  by 
accelerometer  (top)  and  derived  from  Optotrak 
(bottom) 
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Figure  63.  -3  g,  4  Hz  z-axis  shock  at  the  seat  measured,  by 
accelerometer  (top)  and  derived  from  Optotrak 
(bottom)  . 
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64.  Acceleration  waveforms  of  seat  accelerometer  Sx  (full 
line)  and  Optotrak  marker  (dotted  line)  superimposed: 
3  g,  6  Hz  x-axis  shock. 


Figure  65.  Acceleration  waveforms  of  seat  accelerometer  Sz  (full 
line)  and  Optotrak  marker  (dotted  line)  superimposed: 
-3  g,  4  Hz  z-axis  shock. 
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Figure  68.  The  mean  response  (n=4)  of  alkaline  phosphatase  to  a 
two  hour  exposure  in  four  experimental  conditions 
(Exp.  LI  -  Control;  Exp.  L2  -  RMS;  Exp.  L6  -  2g  z- 
axis  shocks  at  32.min-l). 
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Figure  69.  The  mean  response  (n=4)  of  blood  ammonia  to  a  two 
hour  exposure  in  four  experimental  conditions. 
Experimental  conditions  are  described  in  Figure  68, 
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Figure  70.  The  mean  response  (n=4)  of  total  bilirubin  to  a  two 
hour  exposure  in  four  experimental  conditions . 
Experimental  conditions  are  described  in  Figure  68. 
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Figure  71.  The  mean  response  (n=4)  of  the  change  in  blood  voliime 
to  a  two  hour  exposure  in  four  experimental 
conditions.  Experimental  conditions  are  described  in 
Figure  68. 
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Figure  72.  The  mean  response  {n=4)  of  blood  urea  nitrogen  (BUN) 
to  a  two  hour  exposure  in  four  experimental 
conditions.  Experimental  conditions  are  described  in 
Figure  68. 
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Figure  73.  The  mean  response  (n=4)  of  calcium  to  a  two  hour 
exposure  in  four  experimental  conditions. 
Experimental  conditions  are  described  in  Figure  68. 
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ire  74. 


The  mean  response  (n=4)  of  cortisol  to  a  two  hour 
exposure  in  four  experimental  conditions. 
Experimental  conditions  are  described  in  Figure  68 
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Figure  75.  The  mean  response  (n=4)  of  creatine  phosphokinase 
(CPK)  to  a  two  hour  exposure  in  four  experimental 
conditions.  Experimental  conditions  are  described  in 
Figure  68. 
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ire  76.  The  mean  response  (n=4)  of  creatinine  to  a  two  hour 
exposure  in  four  experimental  conditions. 
Experimental  conditions  are  described  in  Figure  68. 
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Figure  77.  The  mean  response  (n=4)  of  glucose  to  a  two  hour 
exposure  in  four  experimental  conditions . 
Experimental  conditions  are  described  in  Figure  68 
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Figure  78.  The  mean  response  (n=4)  of  haptoglobin  to  a  two  hour 
exposure  in  four  experimental  conditions . 
Experimental  conditions  are  described  in  Figure  68. 
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Figure  79.  The  mean  response  (n=4)  of  hematocrit  to  a  two  hour 
exposure  in  four  experimental  conditions. 
Experimental  conditions  are  described  in  Figure  68. 
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ire  80.  The  mean  response  (n=4)  of  hemoglobin  to  a  two  hour 
exposure  in  four  experimental  conditions . 
Experimental  conditions  are  described  in  Figure  68. 
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Figure  81.  The  mean  response  (n=4)  of  lactate  to  a  two  hour 
exposure  in  four  experimental  conditions . 
Experimental  conditions  are  described  in  Figure  68, 
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are  82.  The  mean  response  (n=4)  of  lactate  dehydrogenase 
(LDH)  to  a  two  hour  exposure  in  four  experimental 
conditions.  Experimental  conditions  are  described  in 
Figure  68. 
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Figure  83.  The  mean  response  (n=4)  of  magnesiim  to  a  two  hour 
exposure  in  four  experimental  conditions . 
Experimental  conditions  are  described  in  Figure  68. 
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Figure  84.  The  mean  response  (n=4)  of  myoglobin  to  a  two  hour 
exposure  in  four  experimental  conditions. 
Experimental  conditions  are  described  in  Figure  68 
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Figure  85.  The  mean  response  (n=4)  of  norepinephrine  to  a  two 
hour  exposure  in  four  experimental  conditions . 
Experimental  conditions  are  described  in  Figure  68. 
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Figure  88.  The  mean  response  (n=4)  of  plasma  volume  to  a  two 
hour  exposure  in  four  experimental  conditions. 
Experimental  conditions  are  described  in  Figure  68 
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Figure  89.  The  mean  response  (n=4)  of  platelets  to  a  two  hour 
exposure  in  four  experimental  conditions . 
Experimental  conditions  are  described  in  Figure  68 


(SodfumJ  mEq/l.  (v  [Potassium]  mEq/L 


Potassium 


- ExpL1 

- ExpL2 

-  -  -  -  ExpL6 

—  -  -  ExpL7 


[WBCJ  1000^3 .  <1)  lyygq  1000/tnm3 


White  Blood  Cell  Count  (WCB) 


Time 


94.  The  mean  response  {n=4)  of  white  blood  cell  count 
(WBC)  to  a  two  hour  exposure  in  four  experimental 
conditions.  Experimental  conditions  are  described  in 
Figure  68. 
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Figure  95.  The  mean  response  (n=4)  of  urinary  hydroxyproline  to 
a  two  hour  exposure  in  four  experimental  conditions . 
Experimental  conditions  are  described  in  Figure  68. 
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Figure  102  Effect  of  low-pass  filter  cut-off  frequency  on  the 
seat  to  spine  transmission  ratios  for  3,  2,  St  1  g 
x-axis  shocks . 
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Appendix  F 
Glossary 


fraction  of  critical  damping 
logarithmic  decrement 

time  constant  for  exponential  decay  of  EMG  mean 
power  frequency 

angular  natural  frequency 

positive  x-axis  vibration  or  shock  according  to 
biodynamic  convention:  forward  (ISO  2631,1982) 

positive  y-axis  vibration  or  shock  according  to 
biodynamic  convention:  to  left  (ISO  2631,1982) 

positive  z-axis  vibration  or  shock  according  to 
biodynamic  convention:  upward  (ISO  2631,1982) 

negative  x-axis  vibration  or  shock  according  to 
biodynamic  convention:  backward  (ISO  2631,1982) 

negative  y-axis  vibration  or  shock  according  to 
biodynamic  convention:  to  right  (ISO  2631,1982) 

negative  z-axis  vibration  or  shock  according  to 
biodynamic  convention:  downward  (ISO  2631,1982) 

silver/silver  chloride 

heart  beats  per  minute 

British  Standards 

blood  urea  nitrogen 

cervical  vertebra  #7 

calcium  ion 

centimetre  squared 

creatine  phosphokinase 


DC 

ECG 

EMG 

FM 

g 

g/dl 

GEDAP 

GFR 

GI 

h 

Hz 

lo 


direct  current 

electrocardiogram  or  electrocardiography  'J' 

electromyogram  or  electromyography 

frequency  modulation 

gravitational  acceleration:  9.8  m/s 

grams  per  deciliter 

Generalized  Data  Aquisition/Analysis  Programs 
glomerular  filtration  rate 
gastrointestinal 
hour ( s ) 

Hertz  (cycles /second) 

integrated  EMG  extrapolated  to  time  zero  during 
sustained  contractions 


lEMG 

IRED 

ISO 

ITFC 

IV 

K 

LI 

L2  X 

L3 

L3  y 


rectified  and  integrated  magnitude  of 
electromyogram 

infra-red  emitting  diodes 

International  Standards  Organisation 

iterative  transfer  function  compensation 

intra-venous 

potassium  ion 

lumbar  vertebra  #1 

lumbar  vertebra  #2  with  accelerometer  mounted  in 
x-axis 

Iximbar  vertebra  #3 

lumbar  vertebra  #3  with  accelerometer  mounted  in 
y-axis 


F-2 


L4  z 

L5 

lbs 

LDH 

In 

m 

MANOVA 

MARS 

mEq/dl 

MF 

MFo 


lumbar  vertebra  #4  with  accelerometer  mounted 
z-axis 

liambar  vertebra  #5 
pounds 

lactate  dehydrogenase 
natural  logarithm 
slope  of  lEMG  time  series 
multiple  analysis  of  variance 
Multi-axis  ride  simulator 
millieguivalent  per  deciliter 
mean  power  frequency 

mean  frequency  extrapolated  to  time  zero  of 
sustained  contraction 


mg/dl 

Mg 

min 

ml . 

mmHg 

mmol/L 

msec 

mV 

mV/sec 

MVC 

N 

n 


milligrams  per  deciliter 

magnesium  ion 

minute 

milliliter 

millimetres  of  Mercury 
millimoles  per  liter 
milliseconds 
millivolts 

millivolts  per  second 
maximum  voluntary  contraction 
Newtons 
number 


Na 

ng/ml 

ODAU 

PC 

pg/ml 

RBC 

rms 

RT-AB 

RT-CH 

SAE 

sec 

SST 

SWB 

Sx 

Sy 

Sz 

T1  X 

T2  y 

T3  z 

T4 
T8 
T9  ■ 
TGV 


sodium  ion 

nanograms  per  milliliter 

Optotrak  Data  Aquisition  Unit  (timing  pulse) 

personal  computer 

picograms  per  milliliter 

red  blood  cells  (erythrocytes) 

root  mean  square 

abdominal  circumference  (Respi trace) 
thoracic  circumference  (Respitrace) 

Society  of  Automotive  Engineers 
seconds 

serum  separator  tubes 
Synthetic  Work  Battery 

accelerometer  on  seat  mounted  in  x-axis 

accelerometer  on  seat  mounted  in  y-axis 

accelerometer  on  seat  mounted  in  z-axis 

thoracic  vertebra  #1  with  accelerometer  mounted  in 
x-axis 

thoracic  vertebra  #2  with  accelerometer  mounted  in 
y-axis 

thoracic  vertebra  #3  with  accelerometer  mounted  in 
z-axis 

thoracic  vertebra  #4 
thoracic  vertebra  #8 
thoracic  vertebra  #9 
tactical  ground  vehicle 


TTS 


temporary  threshold  shift 
U/L  international  units  per  liter 

ug/dl  microgram  per  deciliter 

ug/ml  microgram  per  milliliter 

umol/24  hours  micromoles  excreted  in  24  hours 

USAARL  United  States  Army  Aeromedical  Research  Laboratory 

V  Volts 

VDV 
WBC 


WBV 


vibration  dose  value 

white  blood  cells  (leukocytes) 

whole  body  vibration 


